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ABSTRACT 


The problem of determining the vertical axis control, 
requirements for landing a VTOL aircraft on a moving ship deck in 
various sea states i3 examined. Both a fixed-base piloted 
simulation and a non-piloted simulation were used to determine the 
landing performance as influenced by thrust-to-weight ratio, 
vertical damping, and engine lags. 

The piloted simulation wa3 run using a fixed-base simulator at 
N.A.S.A. Ames Research Center. Simplified versions of an existing 
AV-8A Harrier model and an existing head-up display format were 
used. The ship model used was that of a DD963 class destroyer. 

Simplified linear models of the pilot, aircraft, ship motion, 
and ship air-wake turbulence were developed for the non-piloted 
simulation. A unique aspect of the non-piloted simulation was the 
development of a model of the piloting strategy used for shipboard 
landing. This model was refined during the piloted simulation until 
it provided a reasonably good representation of observed pilot 
behavior. Further refinement could lead to a model suitable for 
prediction of landing performance of VTOL aircraft on ships and as 
the basis of control logic for automatic landing. 

A surprising result of this simulation was that, with a good 
station keeping control system and with statistical ship motion 
displayed on the head-up display, pilots could- consistently perform 
safe landings in sea state 6, with handling qualities that were 
adequate at thrust-to-weight ratios greater than 1 .03 and even 



marginally adequate down to thrust- to-weight ratios of 1.01. These 
results should hold quite generally provided that a thrust-to-weight 
ratio of 1 + 4 is interpreted a3 meaning that the pilot always has 
the capability of accelerating the aircraft at Ag upward even in the 
presence of ground effect and hot gas reingestion. 



TABLE OF CONTENTS 


Page 

ABSTRACT 

LIST OF FIGURES ...vi 

LIST OF TABLES x 

LIST OF SYMBOLS xi 

I. INTRODUCTION 

II. ANALYSIS OF HOVER HEIGHT CONTROL 5 

A. CONTROL TASK 8 

B. THE NON-PILOTED SIMULATION.... 8 

C. ANALYTICAL MODELS 12 

1. Aircraft Model 12 

2. Pilot Model 13 

3. Ship Model 15 

4. Turbulence Model 20 

D. RESULTS FOR THE NON-PILOTED SIMULATION 22 

III. SIMULATION EXPERIMENT 25 

A. TEST MATRIX .,.25 

B. AIRCRAFT MODEL... 25 

1. Basic Aircraft ..25 

2. Modifications to the Simulation Model 27 

3. Control Augmentation System 28 

4. Head-up Display........ ......30 

C. SHIP MODEL 32 

D. SIMULATION FACILITIES 32 



TABLE OF CONTENTS , continued 


Page 

1. Simulator . 32 

2. Pilots . 34 

E. OPERATIONAL TASK ,.36 

F. RESULTS FOR THE PILOTED SIMULATION 38 

1, Task Performance and Pilot Ratings ...38 

2. Pilot Comments.... 55 

IV. REVISED ANALYTICAL PREDICTIONS 58 

A. MODIFIED ANALYTICAL MODEL. 58 

B. COMPARISON OF ANALYTICAL AND SIMULATOR RESULTS 63 

V. CONCLUSIONS 75 

VI. REFERENCES 78 

VII. APPENDIX A A1 

A. APPROXIMATE INVERSE LAPLACE TRANSFORM METHOD. A2 

B. A.I.L. IMPLEMENTATION A3 

C. DEVELOPMENT OF THE SHIP MOTION MODEL.... A1 0 

D. TURBULENCE MODEL A1 3 

E. HOVER HEIGHT INPUT A1 4 

F. DATA PLOTTING A14 

G. LISTING OF NON-PILOTED SIMULATION COMPUTER 

VARIABLE DEFINITIONS AND PROGRAM LISTINGS 

VIII. APPENDIX B B1 

A. TEST VARIABLES ... B2 

B. HUD DYNAMICS B3 

iv 



TABLE OF CONTENTS, continued 


Page 


C. DATA OUTPUT OF PILOTED SIMULATION » B4 

I 

D. PILOTS B5 

E. LISTING OF PILOTED SIMULATION RUNS AND ! 

RESULTS B9 j 


i 

i 



V 




LIST OF FIGURES 




Figure No. 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 


10 . 

11 . 

12 . 

13. 

14. 

15. 


16 . 


Title 


Page 


Non-Piloted Simulation Computer Program 

Flow Path .....9 


Block Diagram of the Original Non-Piloted 
Simulation Model... 11 


Block Diagram and Example Time History of 

the Original Flight Path Command Logic ....16 


Wind Waves at Sea 17 

U. S. Non-Aviation Ships 19 


Axis Definitions for the Non-Piloted 

Simulation Model 21 

Non-Piloted Simulation Touchdown Velocity 

Results 23 

Test Matrix for the Piloted Simulation .26 

Block Diagram of Vertical Velocity Damping 

Through Thrust as Implemented in the 

Simulation 27 

Throttle Control Sensitivity 29 

Head-up Display Format and Symboloqy 

Definitions 31 

The Chair 6 Fixed-Base Simulator 

Handling-Qualities Rating Scale 

Mean Flight Times for the Piloted Simulation ....... 38 

Influence of Thrust-to-Weight Ratio, Engine 


and Airframe Dynamics on Touchdown Sink Rate — 

Baseline HUD with Velocity Command Control 

System 


Influence of Thrust-to-Weight Ratio, Engine 
and Airframe Dynamics on Pilot Ratings — 

Baseline HUD with Velocity Command Control 

System 


1 

j 

i 


t 



i 



1 



vi 



LIST OF FIGURES, continued 

Figure No. Title Page 


17. Comparison of Touchdown Sink Rate and Pilot 

Rating as Influenced by Thrust-to-Weight Ratio 
for a Selected Condition — Baseline HUD with 
Velocity and Attitude Control Systems. 42 

18. Influence of Thrust-to-Weight Ratio, Engine and 
Airframe Dynamics on Touchdown Sink Rate — 

Augmented HUD with Velocity Command Control 
System. 44 

19. Influence of Thrust-to-Weight Ratio, Engine and 
Airframe Dynamics on Pilot Ratings — Augmented 

HUD with Velocity Control System ..45 

20. Comparison of Touchdown Sink Rate and Pilot 

Rating as Influenced by Sea State for a 
Selected Test Condition — Augmented HUD with 
Velocity Command Control System.. ....47 

21 . Influence of Thrust-to-Weight. Ratio on the 
Positioning Side Task — Augmented HUD with 

Attitude Command System. ; 48 

22. Distribution of Touchdown Sink Rate and Ship 

Deck Velocity — Baseline and Augmented HUD 
Breakdown for all Landings Completed for the 
Simulation 48 

23. Influence of Thrust-to-Weight Ratio on Thrust- 
to-Weight Ratio Use Histograms — Baseline HUD 

with Velocity Command Control System .50 

24. Influence of Thrust-to-Weight Ratio on Thrust- 
to-Weight Ratio Use Histograms — Augmented HUD 

with Velocity Command Control System 

25. Influence of Engine and Airframe Dynamics, and 
Thrust-to-Weight Ratio on Thrust-to-Weight 
Ratio Use Histograms — Augmented HUD with 

Velocity Command Control System 

26. Influence of HUD Format and Control System on 

Thrust-to-Weight Use for a Selected Test 
Condition 53 


vii 





LIST OF FIGURES , continued 


Figure No. 


Title 


Influence of Sea State on Thrust- to- Weight 
Ratio U3e— Augmented HUD with Velocity Command 
System. 54 

Block Diagram for the Modified Non-Piloted 
Simulation Model 59 


Block Diagram of the Modified Flight Path 
Command Logic 


Example Time History for the Modified Flight 
Path Command Logic..... 


Aircraft and Ship Motion Time History as 
Output from the Non-Piloted Simulation. . . 


Aircraft and Ship Motion Time History as 
Output from the Non-Piloted Simulation... 


Comparison of the Influence of Thrust-to- 
Weight Ratio on Touchdown Sink Rate — 
Augmented HUD with Velocity Command System 
of Piloted Simulation and Non-Piloted 
Simulation. 


Comparison of the Influence of Thrust-to- 
Weight Ratio on Flight Time — Augmented HUD 
with Velocity and Attitude Command Systems 
of Piloted Simulation and the Original and 
Modified Flight Path Command Logic of the 
Non-Piloted Simulation. 


Comparison of the Influence of Thrust-to- 
Weight Ratio and Engine Dynamics on Touchdown 
Sink Rate — Augmented HUD with Velocity 
Command System and Modified Flight Path 
Command System 


Comparison of Thrust- to-Weight Use Histograms 
as Obtained from the Piloted and Non-Piloted 
Simulations 


Variation of Thrust- to-Weight Ratio Use 
Histograms Due to a Statistically Small Number of 
Run3 Used Per Data Point 73 


viii 



LIST OF FIGURES, continued 


Figure Ho. Title Page 


A1 . Block Diagram of the A.I.L. Implementation in 

the Non-Piloted Simulation...... A9 

A2. Examle Output from SHPREF.FOR with Plot of 

Ship Motion Specturm for the DD962 in Heave All 

A3. Example Output from SHPREF.FOR with Plot of 

Ship Motion Spectrum for the DD962 in Pitch A12 

A4. Strip Chart Recording of Vertical Acceleration 

Due to Turbulence as Obtained from the AV-8A 
Fixed-Base Simulation Facilities and Used to 
Determine Magnitude of Turbulence Values for 
the Turbulence Modeling in the Non-Piloted 
Simulation.. A15 

A5. Example Plot of Thrust-to-Weight Ratio Use as 

Output from the Non-Piloted Simulation A17 

A6. Example Plot Showing the Variation in Variance 

and Means for Flight Time and Touchdown Velocity 
as Influenced by the Number of Runs Made in the 
Non-Piloted Simulation A18 

A7. Example, Plot Showing the Variation in Variance 

and Means for Ship Motion as Influenced by the 
Number of Runs Made, Where Number of Runs i3 
Seconds in This Case A19 

A8. Example Output of Run Conditions and Results 

for the Non-Piloted Simulation.... A20 

B1 • HUD Hover Point Dynamics .B4 

B2. Example Data Output from the Piloted 

Simulation B6 

b 3. Piloted Simulation Test Matrix BIO 


ix 



LIST OF TABLES 


Table No. Title Page 


1. Example Lull Criteria....... 3 

2. Simulation Environmental Conditions for the 

D0963 Class Destroyer 33 

3. Touchdown Velocity Statistics Indicating 
Number of Landings with Sink Rates Greater 

Than 12 ft/3 for Selected Test Conditions..... 49 

4. Influence of Thrust-to-Weight Ratio and 

Engine and Airframe Dynamics on Flight Path 
Command Logic Sequence' Use... 62 

A1 . Values for Variables Which Must be Edited into 

Program with Change in Sea State A16 

B1 . Fixed-Base Simulation Computer Variables arid 

Related Variables and Values B2 

B2. Variable Values Use in Correcting for the Mean 

Aerodynamic Forces B3 


x 



list of syh bols 


Definition 

Aircraft transfer function 
Amplitude 

Acceleration magnitude' 

Cosine 

Command 

Degree 

Engine transfer function 
Engine lag time constant 
Random phase angle 
Foot, feet 

Gravitational constant 
Altitude 

Significant wave height 
Inches 

Ship motion component 
Pilot gain 

Nautical miles per hour 

Meters 

Millisecond 

Random number 

Laplace variable 

Second 

Second 


xi 



Greek 

Symbol 

A 

(i) 

a 

x 

A 


LIST OF SYMBOLS, continued 
Definition 

Time (seconds) 

Thrust 

Time increment 
Lag time constant 
Lead time constant 
Modal period 
Velocity 
Weight 

Position along the X axis 

True longitudinal acceleration 

Estimated true acceleration 

Estimate of high frequency acceleration 

Pilot transfer function 

Vertical velocity damping coefficient 


Increment 

Frequency 

Statistical variance 
Time constant 

Component of the ship motion spectrum 



xii 


LIST OF SYMBOLS, continued 


Greek 

Symbol Definition 



Roll angle (deg) 

e 

Pitch angle (deg) 

4» 

Yaw angle (deg) 

<P 

Ship heading (deg) 

c 

Laplace variable 

Subscript 

c 

Command 

eg 

Center of gravity 

en 

Encounter 

eng 

Engi ne 

* P 

Landing pad 

long 

Longitudinal 

v - 

Velocity 

w 

Wind 

wod 

Wind over deck 

X 

Position along x axis 

y 

Position along y axis 

Z 

Position along x axis 

<{> 

Roll axis 

0 

Pitch axis 

* 

Yaw axis 


xiil 





xiv 



I. INTRODUCTION 


The problem of landing V/STOL aircraft aboard destroyer class 
ships has been investigated in the past (References 1-3). Several 
methods have been used to determine the feasibility of, and the 
control/display systems' needed, to accomplish this task. 

Many of the researchers in this area began with the premise 
that for the successful completion of this task, it would be 
necessary for the pilot/aircraft system to have the capability of 
in-phase chasing of the ship deck. The vertical task then was to 
3tart at a specified altitude, descend at a reasonable rate and 
begin to match the vertical motion of the ship deck. If the ship 
deck motion can be matched in both phase and amplitude, then it is 
only a matter of establishing a small relative descent rate and a 
reasonable landing can- be made. There are, of course, problems with 
this technique. In high sea states (5-6) the frequency of the ship 
motion is near the maximum piloting frequency (~4 rad/sec) . In 
attempting to match both the phase and amplitude of the ship motion, 
the pilot is forced to operate at close to his break frequency and 
at fairly high gain. The aircraft must incorporate a high thrust- 
to-weight ratio to achieve the maximum amplitudes in the time 
required. The combination of high piloting frequency and gain with 
high thrust-to-weight ratio can cause lags and a tendency to 
overshoot. If in addition there is a large system lag due to engine 
spool time, display lags, and pilot delay times, the phase lag can 
become excessive to the point of producing large touchdown 
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velocities and an unstable 3ystera. This has been demonstrated in 
computer simulations (Reference 1). As a result it is generally 
concluded that deck chasing is not a reliable method for landing an 
aircraft with reasonable accuracies or consistently low touchdown 
velocities. This conclusion is confirmed in helicopter operations 
onto small ship3. 

A second approach was based on the idea that the pilot could 
loiter until a lull in the ship motion occurred. Some of the 
research indicates that adequate lulls are not frequent occurrences 
or are too short in duration to-be useful. For example Reference 1, 
which investigated results for the two lull criteria given in Table 
1, determined that for the more conservative criterion no lulls 
occurred in the DD963 ship motion model over a period of 1800 
seconds. For the less stringent criterion, 52 lulls occurred in 
this time period, or the average of 1 every 33 seconds. This 
indicates lulls of very short duration. 

Reference 1 concludes that looking for lull conditions under 
high sea states in order to make a landing is not very feasible. 
Other research conducted in the area of lull prediction (Reference 
2) indicates that lull conditions (defined as the time from which 
there are 2 successive peaks under the mean value of the positive 
peak amplitude envelope until 2 successive peaks over that value) 
occurs at the rate of 1 every 70 seconds and are of 10 to 60 seconds 
duration. Another consideration is found in Reference 3 where a 
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Table 1 s Example Lull Criteria 


Reference i LULL CRITERIA 


Motion 

Component 

Criteria 

Criteria 

Liait 

1 

2 

Longitudinal Vel . 

2 ft/sec 

3 ft/sec 

Lateral Velocity 

2 ft/sec 

3 ft/eec 

Vertical Velocity 

2 ft/sec 

3 ft/sec 

Pitch 

1° 

1.5° 

Roll 

2° 

3° 

Pitch Rate 

2 °/sec 

8 °/sec 

Roll Rate 

2 °/sec 

8 °/sec 

No. of Occurrances 

0 

52 

Based on 1800 sec . 
of DD963 ship model 

motion for 

Sea State 5. 


description of sea trials performed using a SH-2F helicopter 
indicated that pilots were often unable to determine visually wherr a 
lull was occurring. The above information indicates that in general 
it is not practical for fixed wing VTOL aircraft to loiter for the 
required time periods in the high fuel use state of hovering while 
waiting for optimum landing conditions. 

Another area of research involved the use of ship notion 
prediction schemes. Research in this area has shown some promising 
results. Computer studies have shown that ship motion can be 
predicted with reasonable accuracies for 10-15 seconds in advance 
(Reference 4). Given this capability, it has been demonstrated in 
computer studies using optimal control modeling techniques, that 
autopilot landings can be made with touchdown velocities on the 
order of 1 ft/sec in sea state 5 conditions (Reference 5). The 
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advantage to using a system in which the ship position i 3 predicted 
in advance and updated as the approach progresses comes from being 
able to adopt a control strategy in which the aircraft is chasing a 
slower moving prediction point rather than the real time deck 
motion. This adds lead, thereby requiring less effort on the part 
of the system a3 demonstrated by adequate performance at thrust-to- 
weight ratios of 1.05. The major problem here is that motion 
prediction for destroyer class ships in high seas has not been 
demonstrated . 

The research conducted for this report is directed at the 
question of what thrus.t-to-weight ratios, and vertical velocity 
damping are required to allow a pilot to make an acceptable landing 
given adequate situational information. 

It is clear that it is desirable to land with low touchdown 
velocities using low thrust-to-weight ratios, and without wasting 
time in the high fuel consumption state of a hovering loiter waiting 
for lulls in the ship motion. 
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1 XI. ANALYSIS OF HOVER HEIGHT CONTROL 


l 


To provide some insight into the problem and the models being 
used in this research, the following facts are presented: 1 ) The 

current AV-8A Harrier aircraft has a maximum touchdown velocity 
limit of 12 ft/sec. 2) The DD963 ship motion model has a maximum 
ship deck heave velocity (3 sigma value) of 17.5 ft/sec. 3) In sea 
state 6, the 3hip has a heave velocity equal to or greater than 
12 ft/sec for less than 0.6% of the time. 4) The maximum heave 
velocity values occur in a region close to the ship motion mean 
position. 5) The lo.wer velocities occur in the peak and trough 
regions of the motions. The above facts in conjunction with the 
ship motion histogram for sea state 6 would indicate that if the 
pilot did nothing but maintain a descent velocity of less than 
2 ft/sec, he would- touch down within the AV-8A gear limits 
approximately 97% of the time. The ship motion statistics show that 
a maximum heave velocity of 12.5 ft/sec i3 encountered at the 2 
sigma heave amplitude for sea 3tate 6 conditions. These same 
statistics show that the 2 sigma heave amplitude values are reached 
at reasonable time intervals of approximately 1 per minute. This 
frequency increases to 1 every 40 seconds for amplitudes of 1 foot 
below the 2 3igma value. These facts suggest a landing strategy 
differing from both the "deck chasing" and "lull waiting." This 
strategy has the pilot descend to the 2 sigma height above the mean 
deck position. At this height the pilot waits and watches the ship 
motion. If it appears that a 3 sigma amplitude (high velocity) deck 
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motion is imminent, the pilot has enough height, and therefore lead 
time to begin an ascent, landing at the more desirable higher 
altitude, (lower ship and relative velocity) as the deck catches up 
to him. If it appears that the deck position is going to peak 
somewhat below his present altitude, there is again enough lead time 
to begin a slow descent and land near the crest (low ship velocity 
position) of the deck motion. The 2 sigma height above mean deck 
position meets the desired criteria for the strategy. It offers an 
easily obtained position without continous deck chasing, provides 
the buffer needed to escape the high velocity portions of the deck 
motion, and presents the pilot with numerous landing opportunities 
without a lengthy loiter period. To accomplish this, the pilot must 
be presented with a suitable indication of mean deck position and 
the 2 and 3 sigma values of deck position relative to the. aircraft's 
landing gear. This information requires measurement of the ship 
motion for several minutes prior to the arrival of the aircraft and 
the transmission of this information to the aircraft for use in the 
head-up display. Real time information is also required to show the 
pilot where the deck is currently positioned within the bounds of 
the probable travel. The pilot can then monitor the ship motion, 
obtain an accurate deck position relative to the mean, and make his 
prediction as to how fast it is moving and where its position will 
be in a couple of seconds. 

An investigation of the proposed landing strategy was conducted 
along two lines: namely, a non-piloted simulation and a piloted 
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simulation. The non-piloted simulation incorporates a linear pilot 
model and vertical axis aircraft model. In addition there is a 
flight’ path command logic section, a command flight path subroutine, 
a ship motion subroutine, and a turbulence modeling subroutine. 

Input variables are- pilot gain, maximum thrust-to-weight ratio, sea 
state, vertical velocity damping coefficient, pilot time delay, and 
engine lag. 

The piloted simulation used the chair 6 simulator at the 
N.A.S.A. Ames Research Center. This fixed-base simulator consists 
of a cab containing the normal, cockpit controls and a single forward 
looking window through which a visual image of the outside 
environment can be obtained. The visual image is provided by a 
camera- terrain board imaging system. A head-up display can be 
superimposed on the outside scene to provide flight situation 
information to the pilot. The cockpit also contains a unique 
throttle/nozzle control quadrant used for V/STOL simulation. 

Slightly modified versions of an existing math model of the AV-8A 
(Reference 6) and an existing head-up display format were used. The 
math model was linearized in the vertical axis and the display 
format simplified to represent only information needed to fly the 
vertical axi3. 




A. CONTROL TASK 


The task to be flown begins with the aircraft at an initial 
altitude of 45 ft above the mean deck position in a stabilized hover 
directly over the bull's-eye on the ship deck landing pad. The 
pilot uses the throttle to control altitude and vertical velocity to 
descend to an initial hover altitude. This initial hover altitude 
is the 2 sigma value of ship deck position above the mean, as 
designated by a line on the head-up display (HUD). The pilot then 
lands at his discretion, based on the ship deck motion information 
presented on the HUD-. 

There are two variations of this task. In the first variation, 
the ship deck motion boundaries and reference lines are not 
displayed on the HUD and the pilot makes the landing without the 2 
sigma reference. In the second variation, an attitude command 
control system replaces the translational velocity command control 
system which creates an effective sidetask in that the pilot must 
actively maintain the aircraft position over the bull's-eye, using 
the control stick, while ; performing the vertical task. 

B. THE NON-PILOTED SIMULATION 

The non-piloted simulation was run through a control program 
which links the main program with flight path, ship motion, and 
turbulence subroutines, data files, and subroutines for output of 
statistical and plotted data. Figure 1 3hows the flow path' diagram 








































between the various programs and files. The main program contains 
the flight path command logic, pilot neuromuscular model, and 
aircraft models. A block diagram showing the model transfer 
functions is shown in Figure 2. 

The rationale for developing a non-piloted simulation was that 
it provided a relatively low cost test of the feasibility of the 
proposed landing strategy and its ability to minimize the required 
T/W. In addition it was conjectured that comparison of the results 
of such a simulation with the piloted simulation would provide a 
test of the validity of the intuitive notions underlying the assumed 
way the pilot would implement the strategy. Perhaps not too 
surprisingly, the task turned out to be much more invloved than 
originally thought, and refinements to incorporate additional 
features suggested by the piloted simulation was a continous 
process. Nonetheless, the insight gained by the exercise was 
invaluable. 

The Approximate Inverse Laplace Transform (A.I.L.) method (see 
Appendix A), was used to solve the differential equations describing 
the pilot and aircraft transfer functions. From a review of the 
literature and looking at the task to be flown, it appeared from the 
beginning that the pilot transfer function would be the major 
problem. For the type of problem being looked at, a generic 
aircraft transfer function could be used, but because of the 
amplitudes and frequencies involved in the ship motion, it was 
^ apparent that the time constants used in the pilot transfer function 
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would be more critical. It was also desired to keep the pilot 
transfer function separated frbm the aircraft transfer function as 
much as possible, so that either the pilot or aircraft portion of 
the program could be moved as a block, either for use in other 
programs or to facilitate looking at other systems in the current 
program. This need to separate the transfer functions complicated 
the A.I.L. equation set-up. 
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__ (g/ELC) 

[S 3 + (1/ELC + Z )s 2 + (Z /ELC)s] 
w w 


A limiter in put on the pilot output of T/W to prevent any input 
above the maximum T/W ratio or below 0; i.e., the pilot can not 
command a negative thrust. 


The engine lag constant, ELC, is one of the parameters on which 


the control requirements depend. 


2. Pilot Model 

The pilot transfer function was first set up using a first 
order lead- lag multiplied by the Pad£ approximation (Figure 2). 

This transfer function was looked at using the Linear Systems 
Analysis Program (LSAP). Bode plots, root locus plots, and time 
history plots were looked at to determine values for lead and lag 
time constants which gave the best results for both a step input and 
a sine-wave input. The values which subjectively produced outputs 
similiar to those expected from a piloted simulation were then used 
in the non-piloted simulation. Because the computer program is 
nonlinear, several runs were made to determine which values for lead 
and lag still gave a good combination of rise time for step inputs 
and small phase and position errors in following a sine-wave. 

Values for the lag time constant, TG, of 0.1 secs and for the lead 
time constant, TL, of 0.5 secs were finally selected. These values 
are within the range of values usually quoted for pilot models 
^ (Reference 7). 
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Because a pure time delay could be programmed into the non- 
pi loted simulation, the Pad€ approximation, which was needed in the 
analysis using LSAP, was replaced. 

The pilot transfer function is 


y = K (TL s + 1) 
p p (TG s + 1 ) 


The pilot's logic is represented by a logic section in the 
simulation. The logic section consists of 4 basic sections. The 
first section is a series of logic statements which determine which 


of the other 3 sections will be used to provide the commanded flight 


path. These sections will be referred to subsequently as ABORT TO 


HOVER HEIGHT, CHASE, and RUN FROM. 

If the aircraft is more than 6 ft above the ship deck, or has 


followed the ship deck below a specified abort chase altitude, or 
has exceeded a specified vertical velocity (a function of the 
maximum available T/W), the flight path command logic enters the 
ABORT TO HOVER HEIGHT section. The hover height is the 2 sigma 
value of the ship deck heave above the ship deck mean position. The 
CHASE sequence i3 entered if the ship deck is within 6 ft of the 
aircraft and the ship deck velocity is less than 2 ft/sec 
(approaching the aircraft) and decreasing. When the CHASE sequence 
is entered, the aircraft flight path is commanded to match the ship 
motion. The RUN FROM sequence is entered if the ship deck velocity 
is greater than 2 ft/sec or if the velocity is increasing; i.e., the 
ship is accelerating toward the aircraft. In the RUN FROM sequence, 
the commanded altitude is the ship position plus an exponential 
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smoothing function (a function of time, based on when the sequence 
was entered). Because the exponential function dies out with time, 
the aircraft is prevented from climbing to excessive altitudes 
before one of the other sequences is initiated. A block diagram and 
example time history of the flight path command' logic prior to the 
piloted simulation are presented in Figure 3. 

After results were obtained from the piloted simulation, models 
representing the time the pilot spends flying a particular portion 
of the task and delays in perception were added. Since the pilot 
scans the situation and instruments and corrects errors in a 
sequence rather than in parallel, the time delay was divided into a 
combination of the pure delay and a gap where' the input was 
maintained at a given value for the time a pilot could be considered 
flying another axis or scanning other instruments and therefore not 
activily flying the vertical task. The length of the time delay, 
the length of the gap, and how often the gap occurs are variables 
that can be initialized at the beginning of a run. In addition 
there are inputs for pilot preception error noise and pilot internal 
noise. These are discussed in more detail later. 

3. Ship Model 

An understanding of sea state can be gained by referring to the 
chart in Figure 4. sea state is shown with the associated wind, 
wave heights, lengths, and periods. Sea state 6 is considered 
significant in that it is estimated that operational capability 
under these conditions would provide for use of aircraft 67% of the 
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SA1, SA2, SA3 are generated by 
the hover altitude subroutine, 
IN4.F0R. SSL1, SSL2, SSL3 are 
generated by the ship motion 
generating subroutine, CENSM.FOR. 


SI Is the command Input to the 
pilot transfer function. 


Aircraft follows flight path 
command to hover altitude. 



51 - SSL1 + 

EXP(1.9*(l-TE/2))-l 

52 - SSL2 

53 - SSL3 



Aircraft within 6 ft. W C 
of descending ship deckA 
enter 'chase' mode. \ 

D) Aircraft descent rate \ 
^ exceeds 1.5*ZHTDT, abort 
1 to hover altitude. . 




Aircraft descends 
below abort chase 
altitude, abort to 
hover altitude. 


J Rapidly rising 
ship deck, nir- 
enters 'run from* 
mode until. In 
this case, touch- 
down occurs. 
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Thla table applies to vavea generated by the local vlnd and dowa not apply to awell orglnatlng elsewhere. 

WARNING NOTE: Presence of swell makes accurate wave observations exceedingly difficult. 

a. The height of waves Is arbitrarily chooscn as the height of the highest 1/3 of the 
waves. Occasional waves by Interference between waves or between waves and swell 
may be considerably larger. 

b. The above values are only approximate due to lack of precise data and to the diffi- 
culty In expressing It In a single easy way. 

c. Below the surface the wave motion decreases by 1/2 for every 1/9 of a wave length 
of depth Increase. 


Figure 4: Wind waves at Sea 
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time, in the North Atlantic daring January. If operations can be 
conducted only under sea state 3 conditions, operations would be 
limited to 31% of the time. 

It should be noted that ship motion may be considerably greater 
than indicated by the wave amplitude values. For example in sea 
state 6 the wave height maximum is 20 ft from crest to trough. The 
heave motion for the DD963 class destroyer can approach 40 ft from 
crest to trough under sea state 6 conditions. 

The DD963 Spruance class destroyer model was used for both the 
piloted and non-piloted simulation. This ship' was chosen because a 
motion program for use with the piloted simulation was already 
available. The DD963 is considered to be typical of the- type of ship 
from which VTOL operations could be conducted. The non-piloted 
simulation was set up so that data for other ships, as contained in 
Reference 8, could also be used. To provide some perspective of the 
ship used, Figure 5 shows a listing of ships according to type and 
class, information on the number of helicopters (assumed replaceable 
by VTOL aircraft' that can be carried, dimensional information, and 
a chart showing displacement. 

A single degree of freedom ship model was programmed using the 
method outlined in Reference 9. The data for the DD963 class 
destroyer were obtained from Reference 8. The model consists of the 
superposition of twelve sine waves, six representing the heave 
motion at the ships center of gravity and six representing the pitch 
motion. The pitch motion is multiplied by the appropriate moment 
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SHIP 

TYPE 

CLASS 

HELI- 

COPTER 

DIMENSIONS 

RELIANCE 

C. G. CUTTER 

WMEC 615 

1 

210.5 

34 

10.5 

BEAR 

C. G. CUTTER 

WMEC 901 

1 1 

270 

38 

13.5 

GARCIA 

FRIGATE 

FF 1040 * 

1 1 

414 

44 

24 

WEST WIND 

ICE BREAKER 

WAGB 83 

2 

269 

63.5 

29 

O.H. PERRY 

FRIGATE 

FFC 7 * 

2 

; 

445 

45 

24.5 

KNOX 

FRIGATE 

FF 1052 * 

1 

430 

46.8 

24.8 

MACKINAW 

ICE BREAKER 

WAGB 83 

1 

290 

74 

19 

SPRUAHCE 

DESTROYER 

DD 963 ** 

2 

563.2 

55.1 

29 

BELKNAP 

G. MISSLE CRUISER 

CG 26 

■ 1 

547 

54.8 

28.8 

GLACIER 

ICE BREAKER 

WAGB 

2 

309.6 

74 

29 

AUSTIN 

AKPHIB. TRANS. DOCK 

LPD 4 

6 

570 

100 

23 

MARS 

COMBAT STORES 

AFS 1 

2 

581 

79 

24 

IWO JIMA 

AMPHIBIOUS ASSUALT 

LPH 2 *** 

11-20 

602 

84 

26 

KILAUEA 

AMMUNITION 

T-AE 26 

2 

564 

81 

28 

WICHITA 

REPLENISHMENT OILER 

AOR 1 

2 

659 

96 

33.3 

BLUE RIDGE 

AMPHIBIOUS COMMAND 

LCC 19 

1 

620 

82 

29 

SACRAMENTO 

FAST COMBAT SUPPORT 

AOE 1 

2 

793 

107 

39.3 

SAIPAN 

AMPHIBIOUS ASSUALT 

LHA 2 

19-26 

820 

106 

26 

— 

AMPHIBIOUS ASSUALT 

LHD 1 

19-26' 

840 

106 

26 


DISPLACEMENT 


* OTHER SHIPS CONTAINED ItJ REF. 4 DATA BASE 

** SHIP MODELED IN THIS SIMULATION 

*** THIS SHIP IS ALSO MODELED AT N.A.S.A. AMES 


TONS X 1000 


Figure 5: U. S. Non- Aviation Ships 



arm and the vertical component is extracted and added to the heave 
motion to obtain the heave motion at the landing pad. Reference 9 
reports that this method gives ship motion accuracies to within 5% 
of a model containing 30 superimposed sine waves per axis. 

This ship motion program uses the Bretschneider wave spectrum, 
transformed based on ship velocity and heading relative to the 
waves. The ship motion spectrum is then obtained by combining the 
wave spectrum with the ship response functions and phase 
differences, and the six amplitudes, Ajj, and frequencies, a> en , 

are extracted. The appendix contains more information on how this 
is accomplished. A component of the ship motion is then represented 
by 

6 

i ( t ) = I A cos (ut-$. . +E) (5) 

. n en ii n 
n=1 

The phase differences' are directly available from the ship motion 
data base information contained in Reference 8. E n is a random 
phase angle which is obtained from a random number generator with an 
output scaled to give valued between 0 an 6.242 radians. Figure 6 
shows a graphical definition of the axis system used in the 
simulation. 

4. Turbulence Model 

The turbulence model (Reference 9) consists of white noise, 
shaped by the following filters: 
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NON-PILOTED SIMULATION 
AXIS DEFINITION 


NON-PILOTED SIMULATION SHIP 
OREIt.TATION TO WAVES 


Figure b: Axis Definitions for the Non-Piloted Simulation Model 
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• 2 1/2 
a r = WA R + cr v (2w n ) n(1/TC ' ) 


A = A + TC A 
R R R 


( 6 ) 


Where Aj^ is turbulence induced accelerations on the aircraft in 
g's. The bandwidth, was obtained from Reference 9. Values of 
the variance, a v , were obtained from a strip-chart recording of the 
turbulence induced vertical acceleration for the AV-8A during a 
fixed-based simulation. The term containing the simulation time 
increment, TC, is a correction to the power spectrum to allow for 
digitization. The quantity, n, is a random number from a Guasian 
distributed sequence with zero mean and unity variance. 


O. RESULTS FOR NON-PILOTED SIMULATION 


Touchdown velocities obtained from the non-piloted simulation 
incorporating the flight path command logic devised prior to the 
piloted simulation are shown in Figure 7. This figure show3 the 
average over many runs during which damping constants, and various 
parameters in the flight path command logic, as well as pilot gain, 
were being manipulated to roughly determine the range of values that 
could be expected for touchdown velocities. These runs were also 
monitored as they were occurring to determine (subjectively) which 
values provided realistic piloting responses. As such, the results 
should be interpreted as a rough indication of the trends. 


22 





Figure 7 shows that in sea state 6 conditions, the orginal 
flight path command logic suggested the very surprising result that 
the landing task could be accomplished without ever exceeding a 
touchdown velocity of 12 ft/sec, for T/W values down to 1.01 and 
with engine lags as high as 0.7 sec. Moreover, the mean value of 
touchdown velocity was only mildly dependent on T/W and engine lag, 
ranging between 5.5 ft/sec and 6.5 ft/sec. Flight times from the 
initial altitude of 20 ft above the mean deck height were again only 
mildly dependent on T/W and engine lag, ranging from 20 to 30 sec. 

The unexpected nature of the preliminary non-piloted simulation 
results strongly indicated the need for experimental verification. 
This experiment was performed on a fixed-base simulator at N.A.S.A. 






Ames Research Center using an existing model of the AV-8A. It was 
recognized from the outset that an important requirement of this 
experiment was the need to supply the pilot with exactly the same 
information assumed in the construction of the command logic in the 
non-piloted simulation. Fortunately, an existing head-up display 
format termed SUPAR HUD (Reference 10) was available and was 
modified by removing all augmented information and adding the basic 
situation information of mean deck position, and the 2 sigma and 3 
sigma deck positions relative to the mean. 
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III. SIMULATION EXPERIMENT 

t • 

1 A description of the test matrix, aircraft and ship models, 

j simulation cab, and results obtained for the piloted simulation are 

1 

» 

presented in the following sections. 


A, TEST MATRIX 

Based on the information gained from the non-piloted 
simulation, the test matrix shown in Figure 8 was set up. Sea 
state, thrust- to-weight ratio, vertical velocity damping through 
thrust, pilot, and HUD format were the variables. The task was then 
simulated using the fixed-base simulation facilties as further 
descibed below. 

B. AIRCRAFT MODEL 
1. Basic Aircraft 

The AV-8A model used in the piloted simulation is described in 
Reference 6. It is a real time digital computer program developed 
to simulate the take-off and landing of V/STOL aircraft aboard 
ship. The unmodified aircraft model includes nonlinear 
aerodynamics, engine and reaction control systems, stability 
augmentation and actuator dynamics, and a landing gear model. 
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Figure 8: Test Matrix for the Piloted Simulation 



















































2 . Modifications to the Simulation Hod e 1 


Only the turbulence portion of the airwake simulation was 
used. In order to create a simulation environment as close as 
possible to that assumed in the non-piloted simulation, the mean 
wind and its variation with height above deck were not used. In 
addition, the ground effects inherent in the AV-8A were deleted from 
the simulation to match the non-piloted simulation, and to more 
effectively isolate the parameters being studied. 

A block diagram of the vertical velocity damping through thrust 
is presented in Figure 9, along with a definition of the variables 
used. 



CHDOT - POWER LEVER ANCLE <deg)/(ft/sec) 

A/C CAIN - 0.45 <ft/sec J ) /POWER LEVER ANCLE (deg) 


Figure 9: Block Diagram of Vertical Velocity Damping 

through Thrust as Implemented in the Simulations 
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3. control Augmentation System 

la running the simulation, the throttle quadrant was used with 
a fixed range of travel throughout the test matrix. A3 a result, 
the throttle sensitivity (vertical acceleration per degree of 
throttle) varied with thrust- to-weight ratio. A plot of commanded g 
per degree of throttle travel may be found in Figure 10. 

Two control systems were used for the fixed-based simulation. 
For the portion of the simulation in which only the vertical axis 
was being flown, a translational velocity command system was used. 
With this system the stick controls velocity in the longitudinal and 
lateral directions. When the pilot leaves the stick centered, the 
aircraft maintains position (zero velocity). The pilot only had to 
control the throttle to perform the task. 

The second control system was an attitude command system. It 
was used to provide a positioning sidetask for the pilot. Using 
this system, the pilot flew the vertical task using throttle and had 
to actively maintain position over the ship deck by commanding 
attitude through the stick. 
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4. Head-Op Display 

Two HUD displays were used in the piloted simulation. These 
are referred to as HUD1 , and HUD 3, and are shown in Figure 11. HUD 3 
was used most extensively in the testing. HUD3 contains the ship 
position reference lines and aircraft hover height and "abort chase" 
altitude lines. HUD1 contains only the symbol indicating current 
ship deck position. HUD1 was used to determine how much (if any) 
advantage there was to displaying the extra information on ship 
motion boundaries. It should be noted that the HUD3 format 
collapses into the HUD1 format for calm seas (sea state 0). 

The head-up display superimposes vertical and horizontal 
situation information. The trident symbol is fixed in the center of 
the display and shows the aircraft's vertical position as the 
distance of the three 'pads' of the trident from the top of the ship 
deck reference symbol. In addition, the trident provides horizontal 
situation indications to the pilot when actively flying the sidetask 
with the attitude command control system. When the sidetask is 
flown, a dagger is added to the display to indicate the undisturbed 
position of the deck bull's-eye. The distance between the trident 
and the dagger shows the error in lateral and longitudinal 
position. The HUD format also contains symbols showing a horizon 
line, pitch bars, pitch reference, and side slip. The lines added 
for the experiment show the 3 sigma, 2 sigma, abort chase height, 
and mean deck positions. 
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The equations describing the transfer function for moving the 
dagger relative to the trident are described in appendix B. To 
assist the pilot in performing the station keeping task the dagger 
motion is augmented with input of aircraft true velocity, the 
estimated acceleration, and stick position. The HUD symbology was 
driven by a Digital Corporation PDP-1VES computer at an update frame 
time of 110 msec. 

C. SHIP MODEL 

Ship dynamics are modeled as six degree of freedom sinusodial 
motion. The ship is assumed to have a fixed mean position about 
which it oscillates. Wind over the deck is composed of a steady 
induced wind equal to the ship speed plus a separate North or East 
component of natural wind which can be specified. At present a 
turbulence model developed for the DD9G3 class destroyer is used 
(Reference 8). This subroutine calculates the free air turbulence 
as well as ship wake turbulence. The latter varies with position 
relative to the landing pad. Table 2 gives an indication of the 
environmental conditions for the simulation. 

D. SIMULATION FACILITIES 
1 ■ Simulator 

The fixed-base chair (Ch. 06), is used primarily to develop 
controls and head-up displays for use in VTOL aircraft and 
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helicopters (Figure 12). Configured as a 3ingle-seat cockpit, it is 
equipped with a conventional stick and rudder pedals with adjustable 
trim. The chair is provided with a virtual image TV display, on a 
single forward looking window. The outside scene is provided by a 
camera- terrain board system. A head-up display can be superimposed 
on the outside scene to provide flight information, to the pilot. 

The cockpit also contains a unique throttle/nozzle control quadrant 
used for V/STOL simulation. Aircraft dynamics are provided by a 
Xerox Sigma 9 digital computer. The aircraft dynamics are updated 
at a frame time of 55 msec for the AV-8A model. 

2. Pilots 

Two pilots were used for the simulation. Both pilots were used 
in the sea state 6 portion of the test matrix. Pilot B filled in 
the remaining tested portions of the test matrix. 

Pilot A is a research engineer at NASA Ames Research Center. 

He has been flying both the fixed and moving base simulators for 15 
years. Most of this time has been spent using the AV-8A and YAV-8B 
simulation models. He was also heavily involved in simulations of 
advanced lift- fan transport (ALFT) aircraft and the Navy research 
and technology aircraft (RTA). 

Pilot B is also a research engineer at NASA Ames Research 
Center. He has been flying both the fixed and moving base 
simulators for over 4 years. Most of hi3 simulation time has also 
been with the AV-8A and YAV-8B models. He also has some helicopter 


simulation time 
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The pilots were asked for comments and a pilot rating following 
each series of the five runs that made up a data point. The pilot 
rating is based on the Cooper-Harper rating scale (Figure 13). More 
information on the Cooper-Harper scale may be found in Reference 12. 

E. OPERATIONAL TASK 





The task to be flown was the same as described in the non- 
piloted simulation. The aircraft was positioned at an initial 
altitude of 45 ft above the mean deck position, in a stabilized 
hover, directly over the bull's-eye on the ship deck landing pad. 

The pilot, through use of throttle, was then to- control altitude and 
vertical velocity in descending to an initial hover altitude. The 
intial hover altitude was the 2 sigma value of the ship deck 
position, above the ship deck position mean, as designated by a line 
on the HUD (HUD3). The pilot was then to land at his descretion, 
aided by the ship deck motion information presented on the HUD. 

There were two variations of the task. In the first, the HUD 
was reconfigured so that the ship deck motion boundaries and 
position reference lines were not displayed (HUD1 ) . The pilot was 
then to complete the task with only the deck position symbol as a 
reference, in the second variation, the control system was changed 
from the translational velocity command system to an attitude 
command system. Thi 3 was done to provide the pilot with the 
sidetask of actively using the stick to maintain the aircraft 
position over the deck bull's-eye while performing the vertical 

•jc 
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task 




Figure 13s Handling-Qualities Rating Scale 




F. RESULTS FOR THE PILOTED SIMULATION 


1 . Task Performance and Pilot Ratings 

A record of the test matrix and runs completed may be found in 
Appendix B. 

As might be expected from the random nature of the ship motion, 
the time, to complete the task wa3 the most variable parameter. The 
average time over each series of 5 runs as a function of sea state, 
HUD, and T/W (Figure 14) lies in a band from 19 to 36 seconds. 



Individual times varied from a low of 7.2 seconds for a particular 
sea state 4, HUD 3 run, to a hiqh of 79.0 seconds for a sea state 6, 
HUD1 run. It should be noted that the addition of the sidetask does 
not appear to change the amount of time the pilot takes to complete 
the task. 
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The non-piloted simulation showed that the mean time and 
standard deviation for flight time did not approach a steady value 
until after approximately 40 runs, indicating why such a large 
variation exits for the piloted simulation groups of 5 runs each. 

Average touchdown velocity as a function of T/W for HUD1 , sea 
state 6 conditions are presented in Figure 15 for two values of 
engine response. The following points should be noted: 1) In 
general there is a less than expected change in touchdown, velocity 
with change in T/W ratio (less than 4.0 ft/sec for the T/W ratio 
range tested), 2) The faster engine time response (0.3 secs) 
produce touchdown velocities of 0.5 to 1.5 ft/sec less than the 
slower responding engine (0.7 sees) and 3) The greater vertical 
velocity damping the lower the touchdown velocity (as much as 4 
ft/sec) . 

Pilot ratings as a function of T/W for HUD1, sea state 6 
conditions are presented in Figure 1 6 for two values of engine 

response. The following items should be noted: 1) The pilot 

rating is higher (worse) as T/w ratio goes down over the range 

tested, 2) Pilot rating is higher (worse) as engine lag is 

increased and as vertical velocity damping is decreased (as much as 
2.5 points difference) 3) With the HUD1 system and low vertical 
damping there are ratings in the Inadequate range even for T/W 
ratios of 1.1. 

Average touchdown velocity and pilot ratings as a function of 
the control system and HUD format are presented in Figure 17. 
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Figure 16: Influence of Thrust- to-Weight Ratio, Engine 

and Airframe Dynamics' on Pilot Ratings — 

Baseline HUD with Velocity Command Control System 
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As shown, the touchdown velocity and pilot ratings are both improved 
with the additional information of the HUD3 format. There is as 
much as a 2.5 ft/sec improvement in average touchdown velocity in 
going from HUD1 to HUD3. It should also be noted that HUD3 produced 
an Adequate pilot rating throughout the range of T/W ratio tested 
for the stated engine response and vertical damping (Figure 17). 

The addition of the sidetask did not seem to influence the result in 
any systematic manner. 

Average touchdown velocity as a function of T/W for HUD3, sea 
state 6 conditions are presented in Figure 18. Again note that 
there is little change in average touchdown velocity over the T/W 
ratio range tested (less than 3 ft/sec for the worse case). Average 
touchdown velocities are 0.5 to 3 ft/sec better with the HUD3 format 
than for the HUD1 format (Figures 15 and 18). The engine time 
constant and vertical velocity time constant have a much smaller and 
less consistent effec.t on average touchdown rate for the HUD3 format 
(Figures 15 and 18). 

Pilot rating as a function of T/W for HUD3, sea state 6 
conditions are presented in Figure 19. The most interesting thing 
to note here is that the average ratings are Adequate for all 
conditions tested throughout the T/W ratio range tested. The pilot 
ratings also, show more clearly that the faster responding engine and 
greater vertical velocity damping are important to the pilot as 
evidenced by the consistent effect they have on pilot rating. There 
is as much as one pilot rating improvement attributable to either a 
fast responding engine' or good damping. 
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Figure 39: Influence of Ihrust-to-Weight Ratio, Engine and 

Airframe Dynamics on Pilot Ratings — Augmented HUD 
with Velocity Control System 
















n: 


Figure 20 provides some perspective on the effect of sea state 
on average touchdown velocity and pilot rating, for a selected, 
condition using the HUD3- format. Note .the relatively small change 
in touchdown velocity over the range of sea state, the greatest 
change occurring for T/W ratio of 1.01 (an increase of 5.0 ft/sec 
from calm sea to sea state 6). The pilot ratings show Satisfactory 
ratings for sea state 0 and 4 and Adequate ratings at sea state 6 
for all T/W ratios for the selected condition. 

The miss distance when performing the sidetask associated with 
the attitude command system is presented in Figure 21 for various 
values of T/W. The miss distance is the horizontal distance from 
f the aircraft center of gravity to the center of the bull's-eye on 

the ship deck at the time of landing. It is interesting to note 
that the average miss distance changes less than a foot with the 
fast responding engine while changing 5.5 feet for the slow 
responding engine over the T/W ratio range tested. In both cases 
the miss distance decreases with decreasing T/W ratio. Possible 
reasons for this result would include; too much control sensitivity 
at higher T/W, or a change in piloting technique; i.e., the pilot 
may be spending a different percentage of time, or changing the 
frequency with which he samples the information at the different T/W 
ratios. 

A better understanding of the relative performance of HUD1 and 
HUD3 can be gained from Figure 22, which shows a histogram of the 
landing performance for HUD1 and HUD 3 with a histogram of the ship 
■ — deck motion for sea state 6 conditions. 
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Influence of Thrust-to-Helght Ratio on the 
Positioning Side Task — Augmented HUD with 
Attitude Command System 
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Figure 22: Distribution of Touchdown Sink Rate and Ship 

Deck Velocity — Baseline and Augmented HUD 
Breakdown for all Landings Completed for the 
Simulation 
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The velocity histogram for HUD1 peaks at 8.5 ft/sec, which is 
just above the 1 sigma value for the ship motion velocity for sea 
state 6. The touchdown velocity for HUD3 peak3 at 6.5 ft/sec or 2 " 
ft/sec less than HUD1. The histogram also shows the HUD3 data to 
have a smaller standard deviation. 

Table 3 indicates the precentage and number of landings over 12 
ft/sec which occurred overall and for several sets of test 
conditions. 


Table 3: Touchdown Velocity Statistics Indicating Number 

of Landings with Sink Rates Greater Than 12 ft/s 
for Selected Test Conditions. 


CONDITION 

NUMBER OF 
LANDINGS 

NUMBER OF 
landings 

> 12 fc/sec 

X 

LANDINGS 
> 12 ft/sec 

Total Simulation 

995 

28 

2.81 

S.S.O, S.S.4 Overall 

260 

0 

0.0 

S.S.6 Overall 

735 

28 

3.81 

S.S.6, HUD1 Vel. Cad. Sys. 

212 

18 

8.49 

S.S.6, HUD3 Vel. Cad. Sys. 

419 

7 

1.67 

S.S.6, 1NJD3 Atd. Cad. Sys. 

104 

3 

2.89 

S.S.6, HUD3 V.C. t -0.3 

216 

1 

0.46 

S.S.6, HUD3 A.C. T -0.3 
eng 

60 

1 

1.67 


A series of T/W ratio histograms are presented in Figures 23 through 
27. Most of these curves show a fairly sharp spiked peak. In 
observing the simulation and through pilot comments, it is concluded 
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that this shape of curve can be attributed to the pilots technique 
of setting up an initial descent rate and holding it either until 
near the ship deck for the HUD1 system, or until the 3 sigma line on 
the HUD3 system was approached, and then slowing or arresting the 
descent. The initial descent phase, which takes 90 to 95% of the 
total time produces the predominant peak in the histogram. 

The histograms for HUD1 , sea state 6 conditions are presented 
in Figure 23. 


HUD1 

Sea State 6 
Velocity Cad. Sys 
Z “ -0.4 sec -1 

t - 0.3 sec 
eng 


T/W 

max 

01.10 
□ 1.05 
& 1.01 


percent 

OF 
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THRUST-TO-l.’ZIGHT RATIO, T/U 


Figure 2! Influence of Thrust-to-Weight Ratio on 
Thrust- to-Weight Ratio Use Histograms — 
Baseline HUD with Velocity Command Control 
System 


The peak occurs at the T/W ratio used most often in the initial 
descent. The peak is somewhat broader for T/Wmax of 1.1 probably 


E 



indicating that the pilot is operating with lower gains. The base 
of the peak is fairly narrow and corresponds to the observed 
pilot behavior, with the HUD1 format, of chasing the deck less, 
possibly due to lack of positional cues. Landings, using the HUD1 
format generally occurred as the deck caught the aircraft. The area 
under the curves tends to bunch up at the high end of the available 
T/W as T/Wraax i3 decreased. This is due to the pilot tending to fly 
more conservatively (lower descent velocities). 

The histograms for HUD3, sea state 6 conditons are shown in 
Figure 24. 
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Figure 24: Influence of Thrust- to-Weiaht Ratio on 

Thrust- to-Weight Ratio Use Histograms — 
Augmented HUD with Velocity Command Control 
System 







Notice that the curves broaden somewhat in comparison to those 
for the HUD1 format. Because of the positional information 
available to the pilot with the HUD3 format, pilot gains are 
somewhat lower, and the pilot ha3 more opportunity to chase or run 
from the ship deck. This tends to broaden the peak and base of the 
curves. The area under the curve for T/Wmax = 1.01 is much more 
bunched up at the high end, than for the HUD1 case. This is due to, 
the pilot's tendency to fly more conservatively (lower descent 
velocities, higher thrust settings), since the HUD3 format makes him 
more aware of the limitations of the available control power. 

The effects of vertical velocity and engine time lag for HUD3, 
sea state 6 conditions are presented in Figure 25. 

With zero vertical damping, the pilot had a much more difficult 
time controlling vertical velocity, and this is compounded with a 
slow responding engine. The curve becomes much less peaked and is 
spread out over a larger range of available T/W. This is in 
contrast to the curve for higher damping and faster responding 
engine where the pilot is able to control vertical velocity with 
much less throttle movement (i.e., the throttle becomes a vertical 
velocity command control) . 

The effects of HUD format and type of control system on T/W 
histograms is presented in Figure 26. An interesting thing to note 

here is that the peak for the attitude command system occurs at a 

higher T/W ratio than for the velocity command system using either 

HUD. Also the curves for the attitude command system tend to be 
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Figure 25: Influence of Engine and Airframe Dynamics, 

and Thrust-to-Weight Ratio on Thrust-to-Weight 
Use — Augmented HUD with Velocity Command Control 
System 
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Figure 26-. influence of HUD Format and Control System on 

Thrust-to-Weight Use for a Selected Test Condition 
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much more symmetrical about the peak. The explaination for this is 
that the pilot was observed to spend at least as much time on the 
positioning sidetask as with the vertical task. As a result, he 
tended to fly more conservatively (at a higher thrust, slower 
descent rate) and the errors in desired vertical position/velocity 
are therefore more likely to occur randomly rather than only on a 
conservative side as when flying only a vertical task. 

The effect of sea state on T/W histograms for a selected 
condition is presented in Figure 27. 

As expected, at the low sea states the pilot spends more time 
in the initial descent and less time correcting for deck position. 
As a result, lower sea states produce a histogram with a sharper 
peak and with a smaller base. 



Figure 27: Influence of Sea State on Thrust-to-Weight Ratio Use- 

Augmented HUD with Velocity Command System 



2. Pilot Comments 


The pilots were given an explanation of the task and 
discr'iption of the HUD symbology. From this they developed some 
individual techniques. Pilot A, especially for the higher values of 
T/W with fast responding engine, often waited for the deck symbol to 
crest just below the 2 sigma line., and then smoothly rolled off the 
throttle providing a quick descent in which the aircraft would catch 
the deck on its downward motion with a usually low value of relative 
descent rate. Pilot B used this technique also although, not as 
often. 

Both pilots became more adept at picking out "lulls" in the 
ship motion as the simulation progressed. It was often possible to 
tell that the deck was in a lull condition when its motion was slow 
and position was a couple of feet above the mean line. The pilots 
would make a quick descent and attempt to catch the deck before the 
more extreme motions reoccurred. 

The pilots also understood that if they chased the deck below 
the mean line (or were unable to arrest a descent until below the 
mean line), then in most cases it was better to continue and attempt 
to catch the deck near its lowest position (and therefore low 
velocity) than attempt to climb back to the 2 sigma line. When the 
attempt to pull up was made, especially with low T/W, the deck 
tended to catch the aircraft near the mean deck position, often with 
a high velocity, and therefore, high relative velocity. 





The most general comment expressed about technique was, "always 
be in a position to gradually take off power, and don't get caught 
needing it." 

Pilot comments indicated the following: 

1 ) The greater the T/W, the more controllable and easier it is 
to perform the task (through the range tested). 

2) The higher the value of T/W, the less sensitive the pilot 
workload is to engine lag and vertical damping. 

3) The higher values of damping provided better control of 
vertical velocity, which in turn aided the initial descent to the 2 
sigma line. The lower values of damping provide quicker vertical 
response and therefore greater agility during the final phase of 
descent. 

4) The slow responding engine was considered unfavorable, even 
though a few cases occurred when the engine compensated for an 
overcontrol by the pilot. 

5) Both pilots commented on the importance of engine noise as 
a cue. 

6) HUD 3 was perferred over HUD1 . The pilots commented 
favorably on having the ship motion boundaries and 2 sigma lines as 
references in giving precise situation information. 

7) Comments concerning HUD1 generally focused on the feeling 
of not knowing either the position or vertical velocity of the 
aircraft relative to the mean deck position. 
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8) When using HUD1 , the effects of damping and engine lag were 
less evident, since errors were not as easily detected due to the 
lack of references. One of the pilots commented that the task 
workload was less using HUD1 because of the lack of references, but 
gave it a higher (worse) pilot rating because of the greater 
uncertainty involved. 



XV. REVISED ANALYTICAL PREDICTIONS 


After the piloted simulation was run, the non-piloted 
simulation was modified in an attempt to match the measured T/W 
histograms, as well as to better represent the observed piloting 
technique. 

A. MODIFIED ANALYTICAL MODEL 

The following modifications were made to the program: 

1 ) A feed back loop was added to provide vertical velocity 
damping through thrust, and the coefficient for vertical velocity 
damping through airframe was reduced to more closely match the 
fixed-based simulation model. 

2) The flight path command logic was rewritten to provide a 
better match to the observed landing strategy adopted by the 
pilots. Details of the logic are given in the next section. 

3) Two additional noise sources were added to account for 
pilot perception error and internal pilot noise. 

4) The pure time lag, used to represent the pilot’s 
information processing time interval, was divided into a pure lag 
and a secondary time in which the input to the pilot's neuromuscular 
dynamics was held constant for a specified time. The purpose here 
was to simulate the pilots concentration on a sidetask. A block 
diagram of the modified model set-up is shown in Figure 28. 

















The modified version of the flight path command logic consists 
of 6 basic sections. The first section determines which of five 
defined regions of relative position and velocity that the aircraft 
is currently in. This section then specifies which of the remaining 
5 sections is to be used to supply the commanded flight path 
information to the pilot transfer function. These sections are 
referred to according to their basic strategy: RUN FAST, RUN, 

CHASE, ABORT TO HOVEl HEIGHT, and CHOP THROTTLE. The RUN FAST 
sequence is initiated whenever the relative velocity of the aircraft 
and ship exceeds a given value, typically 5.5 ft/sec. The RUN FAST 
logic commands an altitude of 12 ft above the present altitude using 
a cosine smoothing function. This command causes the simulated 
pilot to apply full throttle. If the relative position is less than 
9 ft and the relative velocity exceeds 4.5 ft/sec the RUN sequence 
is initiated. The RUN logic commands an altitude of 6 ft above the 
present altitude through a cosine smoothing function, causing the 
simulated pilot to apply approximately 95% of full throttle. If the 
aircraft was in the RUN FAST sequence, and then switches to the RUN 
sequence, the commanded altitude is lowered 6 ft from the previously 
commanded altitude, thus causing the throttle to be reduced from 
full to approximately 90%. If the relative position of the aircraft 
and ship deck is less than 3 ft and the relative velocity hasn't 
exceeded 4; 5 ft/sec then the CHASE sequence is initiated. The 
CHASE logic commands a cosine function descent, modulated by the 
ship deck motion, starting at the previously commanded altitude and 
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ending at the ship deck. If the aircraft position descends below a 
designated abort height, the ABORT TO HOVER HEIGHT sequence is 
initiated. The abort height is the height above the ship deck mean 
position which, if a descent is continued, will not provide enough 
time to gain the necessary height to prevent a hard landing in the 
event the next segment of the motion is around the 3 sigma value. 
This height, of course, varies with sea state. Up to sea state 4 
the abort height is zero because the relative velocity can be 
maintained below the gear limits anywhere in the ship motion 
boundaries. The hover height is the 2 sigma height above the ship 
deck mean. The ABORT TO HOVER HEIGHT logic commands the hover 
height altitude through a cosine smoothing function from the 
altitude in the previous sequence. The hover height altitude is 
then maintained until conditions require use of another logic 
section. The CHOP THROTTLE sequence is designed to mimick a normal 
pilot landing technique. It is initiated when the relative aircraft 
to ship position is less than the Chop Throttle Now Height (CTNH), 
which is another adjustable variable. Unlike the other sequences, 
which can be abandoned for a more apporpriate one at any time during 
the sequence, once the CHOP THROTTLE sequence is initiated it 
continues until a landing occurs or a specified time has elasped 
The CHOP THROTTLE logic commands a 12 ft descent in altitude through 
a cosine smoothing function. This command effectively produces near 
zero thrust output. The sequence maintains this reduced altitude 
for a specified period of time. If during this tine period a 
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landing has not occurred, the sequence commands an ascent back to 
the hover altitude. This logic models the pilot behavior after he 
sees a landing opportunity; i.e. , the aircraft is 0.5 ft off the 
deck but with the deck beginning to descend. 

The modified analytical model produced a much more accurate 
representation of observed piloting technique. For example, it was 
observed in the piloted simulation that with the slow responding 
engine, more landings were made running from the ship than chasing 
it, whereas with the fast responding engine, approximately the same 
number of landings were made chasing the ship as running from it. A 
series of computer runs were made which duplicated this result 
(Table 4). 


Table 4: influence of Thrust-to-Weight Ratio and Engine and Air- 

frame Dynamics on Flight Path Command Logic Sequence Use 


TEST CONDITION 

T/W t Z 

eng w 

sec sec -1 

PERCENT OF LANDLNCS 

Flight Path Logic Section 
12 6 7 

1.07 0.7 0.4 

52.5 10.0 

35.0 2.5 

1.05 0.7 0.4 

55.0 10.0 

15.0 20.0 

1.03 0.7 0.4 

55.0 5.0 

27.5 12.5 

1.01 0.7 0.4 

50.0 15.0 

30.0 5.0 

Average 

63.1 

36.9 

1-01 0.3 0.4 

42.5 10.0 

27.5 20.0 

1.07 0.3 0.4 

42.5 5.0 

32.5 20.0 

1.05 0.3 0.4 

40.0 12.5 

30.0 17.5 

1.03 0.3 0.4 

45.0 2.5 

45.0 7.5 

Average 

50.0 

50.0 

Flight Path Logic Sect lo 

ns: 


1 - Run Fast Sequence 

2 - Run Sequence 

6 - First section of the Chop Throttle 

Sequence 

| 7 - Second section of the Chop Throttle Sequence | 



a 
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The original flight path command logic did not provide a chop 
throttle sequence, and the landing nearly always occurred while in 
the HUN FROM sequence. A schematic diagram of the modified flight 
path command logic .s shown in Figure 29. An example time history 
showing some of the flight path command logic aspects is shown in 
Figure 30. In addition, two representive time histories are 
presented in Figures 31 and 32. 

B. COMPARISON OF ANALYTICAL AND SIMULATION RESULTS 

A comparison of touchdown velocities for the non-piloted 
simulation and the piloted simulation are shown in Figure 33. The 
non-piloted simulation results, using either the original or 
modified flight path command logic, produced lower touchdown 
velocities than were achieved i:i the piloted simulation. The 
modified logic generally produced the lowest touchdown velocities. 
Further adjustment of the parameters in the non-piloted simulation 
to produce a closer match to the piloted simulation results could 
provide further insight into the pilot's capabilities. 

The results for the modified version of the flight path command 
logic when compared to the piloted simulation data show a good 
correspondance in average flight time for the lower T/w ratios, but 
a large gap is evident for T/W = 1.1 as shown in Figure 34. 
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Figure 29: Block Diagram of the Modified Flight Path Command Logic 


























Figure 30s Example Time History for the Modified Flight Path Command Logic 


INITIAL LET DOWN TO HOVER 
■ HEICHT 


MAINTAINING HOVER HEICHT 


Relative velocity >5.5 
ft/sec , RUN FAST SEQUENCE 

A/C <3.0 ft above deck 
CHASE SEQUENCE 

A/C hits abort height, 

ABORT TO HOVER HEIGHT SEQUENCE 


Relative velocity >5.5 

ft/sec. RUN EAST SEQUENCE 

A/C <3. ft of deck, CHASE SEQUENCE 

A/C >3.0 ft above deck, 

ABORT TO HOVER HEIGHT V. 

SEQUENCE 


A/C <3.0 ft above deck 
CHASE SEQUENCE 


A/C hits abort height, 

ABORT TO HOVER HEIGHT SEQUENCE 

A/C <9.0 ft above deck, velocity 
of deck >4.5 ft/sec, RUN SEQUENCE 

A/C < CTNH, CHOP THROTTLE . * * ' 
SEQUENCE ’ # / 

A/C >3.0 ft above deck, [ 

ABORT TO HOVER HEIGHT SEQUENCE V 

Relative velocity >5.5 ft/sec, N 
RUN FAST SEQUENCE 



Relative velocity <4.5 ft/sec, 
RUN SEQUENCE, till touch down 
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Figure 31 : Aircraft and Ship Motion Tine History as Output from the Non-Piloted Simulation 
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Figure 32: Aircraft and Ship Motion Time History as Output from the Non-Piloted Simulation 











THRUST-TO-VEICUT RATIO, T/W 

Figure 33: Comparison of the Influence of Thrust- to-Weight Ratio 

on Touchdown Sink Rate — Augmented HUD Velocity Command 
System of Piloted Simulation with Non-Pilated System 



THRUST-TO-WEICHT RATIO, T/U 

Figure 34: Comparison of the Influence of Thrust- to-Ueight Ratio 

on Flight Time — Augmented HUD with Velocity and 
Attitude Command Systems of Piloted Simulation and 
the Original and Modified Flight Path Command Logic 
of Hon-Piloted Simulation 








A comparison of average touchdown velocities for a selected 


condition is shown, in Figure 35. 
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Figure 35: Comparison of the Influence of Thrust- to-Weight Ratio 

and Engine Dynamics on Touchdown Sink Rate — 

Augmented HUD with Velocity Command System and Modified 
Flight Path Command System. 


As can be seen the trends are roughly correct for either 
command logic. However, there is still a fairly large bias that is 
unaccounted fcr. It should be noted that the data shown for the 
non-piloted simulation is a composite of data in which the pilot 
gains, aircraft vertical velocity damping, and variables mentioned 
in the T/W histogram section were being manipulated in an attempt to 
find optima (touchdown velocity and flight time) for each T/W 
ratio. Most of these runs were made shortly after the piloted 
simulations began, but before the data from the piloted simulation 
were analysed. A more precise comparison requires another series of 




non-piloted simulations to be run with variables sirailiar to those 

us”d in the piloted simulations. 

Initially, there was only interest in obtaining an indication 

from the non-piloted simulation as to whether or not the lower 

thrust-to-weight values were practical for landing in high sea 

states. In the early stages the program was used to provide 

estimates of touchdown velocity means and standard deviation and the 

time required to land. As the work progressed it became clear that 

* 

the fixed-base simulation output data used to construct thrust-to- 
weight ratio histograms, could then be used as another matching 
variable in determining the accuracy of the non-piloted simulation. 

There are 17 variables whose values determine in some way the 
shape of the histogram. These are: 

1. The initial starting altitude. 

2. The letdown time. 

3. The hover height altitude. 

4. The abort height altitude. 

5. The chop throttle height. 

6-8. Pilot gains. 

9. The pilot pure lag time. 

10. The pilot sidetask time. 

11. The ratio of 9. and 10. 

12. The pilot preception error noise amplitude. 

13. The pilot preception error noise frequency. 

14. The pilot internal noise amplitude. 



15. The pilot internal noise frequently. 

16. Pilot lead tine, time constant. 

17. Pilot lag time constant. 

Some of these non-piloted simulation variables are easily fixed. 
Initial height, hover height, and abort heights are all displayed 
for the pilot and so are the same in both piloted and non-piloted 
simulations. Pilot internal noise amplitude was based on having 
each pilot, • and any other observers present, guess the velocity at 
which the pilot touched down as viewed on the head-up display 
monitors. The error was then determined and the noise standard 
deviation was set at 1.5 ft/.sec after averaging over the number of' 
observers approximately 180 landings-. The res.t of the variables 
were set by making an initial guess and then making 5 runs and 
looking to determine hbw the variable had changed the histogram. 

This process was repeated by either changing the variable again, or 
going to the next one. In this way a good match to the histogram 
was made for one of the cases (Figure 36). It should be noted that 
there may be as much, as 20-30% change in any one point on the 
histogram for a given set of 5 runs (as can be seen in Figure 37). 
This is because 5 runs are not enough to get a good statistical 
representation, thus aggravating the difficulty of trying to obtain 
a good match. The non-piloted simulation showed that approximately 
40 runs (depending on the T/W ratio being used) is needed to p. ovide 
adequate statistics. When histograms representing only 5 runs are 
constructed based on 10 unknown variables, there is some difficulty 
in producing a good match. Another problem occurs when trying to 
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Figure 37: Variation of Thrust- to-Weight Ratio Use Histograms 

Due to a Statistically Small Number of Runs Used 
Per Data Point 


match the touchdown velocity and flight time averages for the same 
cases. It was found that once the major influence of variable 
changes on the histogram shape was obtained, it was fairly easy to 
get a rough match of histograms. The T/W ratio for which the 
histogram peaks can be obtained through selection of the initial 
letdown velocity of the aircraft. The base can be broadened by 
selecting higher values for the pilot gains. The sharpness of the 
fillets between the base and peak were found to change somewhat with 
the selection of pilot gains and the values used for frequency and 
magnitude in the pilot internal and preception noise models. The 
shape of the base; i.e., the number and size of the peaks appeared to 
be mainly dependant on the ratio of three pilot gains used. It is 
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also fairly easy to get a reasonable match of touchdown velocities 
and flight times. It i3 not easy to get. all three of these results 
to give a reasonable match simultaneously. This may be due to 
something inherent in either the piloting technique or aircraft 
model being, flown on the fixed-base simulation that is not being 
modeled accurately in the non-piloted simulation, or that the right 
combination of values for the variables has not been found. The 
latter should be explored further using parameter identification 
techniques . It would ultimately be hoped that all three could be 
matched for a couple of cases and then the values of the variables 
be determined analytically to produce matching values for other 
cases. This pilot model would then be a good tool for predicting 
the landing performance of any VTOL aircraft onto any type of ship. 
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V. COHCI.USIQNS, 

The problem of determining the vertical axis control 
requirements for landing a VTOL aircraft on a moving ship deck in 
various sea states is examined. Both a fixed-base piloted 
simulation and a non-piloted batch simulation were used to determine 
the landing performance as influenced by thrust- to-weight ratio, 
vertical damping, and engine lags-. 

The piloted simulation was run using a fixed-base simulator at 
N. A. S. A. Ames Research Center. Simplified versions of an existing 
AV-8A Harrier model and an existing head-up display format were 
used. The ship model used was that of a DD963 class destroyer. Two 
pilots were used to obtain data and to give pilot ratings based on 
the Cooper-Harper pilot rating scale. 

A surprising result of this simulation was that, with a good 
station keeping control system and with statistical ship motion 
displayed on the head-up display, pilots could consistently perform 
safe landings in sea state 6, with handling qualities that were 
adequate at thrust- to-weight ratios greater than 1.03 and even 
marginally adequate down to thrust-to-weight ratios of 1.01. These 
results should hold quite generally provided that a thrust-to-weight 
ratio of 1 + A is interpreted a3 meaning that the pilot always has 
the capability of accelerating the aircraft at Ag upward even in the 
presence of ground effect and hot gas reingestion. 

Preliminary work with a non-piloted simulation showed that with 
. a good strategy and the right information, a pilot should be able to 
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land a VTOL type aircraft vertically aboard a DD963 class destroyer 
under sea state 6 conditions, in an adequately controllable manner, 
with thrust- to-weight ratios as low as 1.01, engine lags as high a3 
0.7 sec, and vertical velocity damping of 0.2 secs, without 
exceeding a 12 ft/sec landing gear limit. This non-piloted 
simulation showed an overall average touchdown velocity of 5.8 
ft/sec and an average flight time of 32.5 seconds. Results were 
then obtained from a piloted fixed-base simulation in order to 
verify the non-piloted results. Similiar results were obtained, 
with an average touchdown velocity of 6.7 ft/sec and flight time 
average of 33 seconds. In addition, the pilot ratings indicate 
satisfactory (level 1) handling qualities for sea state 6 conditons 
and thrust-to-weight ratios as low as 1.03 and adequate (level 2) 
handling qualities for thrust-to-weight ratios as low as 1.01. 

Pilot ratings showed the expected results of being more 
favorable as T/W ratio increased, up to the maximum tested of 1.1, 
and with increasing vertical velocity damping, up to the maximum of 
-0.4 sec -1 tested, and with the faster responding engine, engine lag 
of 0.3 secs. The pilots also demonstrated lower touchdown 
velocities when presented with the ship motion boundaries, and 
aircraft hover and abort chase height lines on the head-up display, 
then when presented only with the ship deck position symbol. The 
simulation also showed that the pilot was capable of obtaining 
similiar results when flying either the translational velocity 
command system or attitude command system. Pilot ratings indicate 
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that the translational velocity command system was perferred 
however. 

Based on the current non-piloted simulation, it is believed 
that an extension can be made to determine piloted results for other 
T/W ratios, engine lags, vertical velocity damping, and ship 
classes, under various sea state conditions. This is based on the 
assumption that, using parameter iridentif ication techniques, touch- 
down velocities, flight times, and T/W use histograms can be made to 
match the current piloted simulation data 1 . 

Although the simulation indicates that aircraft can be landed 
vertically at much lower T/W ratios than previously suspected, even 
with the positioning sidetask, it remains to be seen how well low 
thrust- to-weight ratios will work when a full range of sidetasks 
inherent in flying actual aircraft are employed. In addition, more 
work should be done to determine the effects on minimum T/W ratio of 
non-linear elements which were not examined in this simulation. The 
effects of suckdown, fountains, hot gas ingestion, and height 
dependant mean winds can all significantly effect minimum T/W. 

As a practical point regarding the simulation, additional 
research would have to be conducted to ascertain how well 
instrumentation aboard ship can determine the mean ship position and 
the motion boundaries, how long in advance of the aircraft arrival 
the motion would have to be monitored to provide accurate results, 
and how well the aircraft/ship systems can determine real time deck 


position 
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The computer program for the non-piloted simulation was 
discussed in general in the main body of the report. Some more 
detailed information regarding certain aspects of the program are 
described below. 

A. APROXIMATE INVERSE LAPLACE TRANSFORM METHOD. 


The numerical iteration method- used for solving the 
differential equations used in the computer program is the 
Approximate Inverse Laplace (A.I.L.) method (Reference 13). It is 
assumed that the equations for the transfer functions are in the 
form: 


X(s) 


N ( S ) 
DCs) 


A s 11 + A .s" ' + ...+ A. 
n n -1 0 

m m -1 

Bs +B s +...+ B _ 
m m -1 0 


( A1 ) 


where m > n and B m is not equal to zero. 

Then using the following recurrence formula: 


m -1 

C. = — [A . - I B.C. .) 

i B n-i+1 . _ 3 i-m +3 

m j=u 


( A2) 


where A = B =0 for x < 0. and C =0 for x < 0. 
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X(s) can be rewritten as 


. n-m n-m -1 _ n-m -2 

X(s) = C s + C s + C s + .. 

1 2 3 


(A3) 


The inverse transform of equation (A3) is 
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X(t + dt) 1 - C 1 + Cjdt + C 3 


( A4) 


The derivatives can then be found frota differentiation: 

• dt^ 

X(t + dt) - C„ +• C,dt + C. ... 

2 3 4 21 


CAS) 


•• J • 

X(t + dt) « C, + C dt + C. — - ... CA6) 

3 4 5 21 

The A.I.L. method allows the use of time varying coefficients 
in the transfer functior. . The method is based on the assumption 
that a time interval, At, can be found during which all of the 
coefficients can be considered as constant. A set of points 
calculated in one interval of time is then used a3 tne initial 
conditions for the next calculation, and the process is repeated. 
Hore information on the A.I.L. method may be obtained from 
Reference 13. 


B. A.I..L. IMPLEMENTATION 

The specific equations used for the non-piloted simulation were 
derived as follows: 
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Then 
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Using Laplace transformation, 

Y - r 

Y “ aT - Y(0) 

Therefore, 
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dH 
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or r , d Il. 

s 2 + 1/TG 3 

Using' the A.I.L. method: 


Where dH is considered 
a constant for the time 
interval under 
consideration. (A9) 
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(dH will be reintroduced 
at a later point) 
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And from the reversion forumula: 
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C4 A “ -B1 * C3A C4B •• -Bt * C3B 
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Then 

YA » Cl + C2A * IT1 + C3A * IT2 + ... 

YB ** C2B * IT1 + C3B * IT2 + ... 

YA - C2A + C3A * IT1 + C4A * IT2 + ... 

YB = C2B + C3B * IT1 + C4B * XT2 + ... 

For the next iteration 
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Substituting (A15) for Hts 
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Therefore, (A17) can be rewitten as 
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Using Laplace transform: 
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where B is considered a constant for the tine 
interval under consideration. 


Or 
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Using the A.I.L. method: 

AAO = A 

AA1 “i^ z(0) + ( ife + Z w )z(0) + z <°> 
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And from the recursion formula: 
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Cl = A3 


r" 


C2 *» A2 - B3 * Cl 

C3 » A1 - B2 * Cl - B3 * C2 

C4 = AO - B1 * Cl - B2 * C2 - B3 * C3 

C5 = -B1 * C2 - B2 * C3 - B3 * C4 (A25) 


(A26) 


(A27) 

A block-type diagram showing how the A.I.L. equations are 
implemented in the- computer program for the pilot servo and aircraft 
transfer functions is shown in Figure A1 . 

A series of runs were made to determine the number of A.I.L. 
coefficients and size of time step required to obtain accurate 
results and a stable program. It was found that 6 terms and a time 
step of 0.01 seconds provided good results. 

AO 


Then, 




z = Cl + C2 * IT1 + C3 * IT2 + 
z = C2 + C3 * IT1 + C4 * IT2 + 
z = C3 + C4 * IT1 + C5 * IT2 +■ 

For the next iteration. 


z(0) n = z n-1 


z(0) = z , 

n n-1 


z(0) = z , 

n n-1 






















C. DEVELOPMENT OF THE SHIP MOTION MODEL 


To obtain the amplitudes for. the individual sinusoidal 
components, the following method wa3 employed. First, the two- 
parameter Bretschneider wave spectrum i3 used to obtain the wave 
spectrum: 


S w (u) 
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This is transformed using the following relationship: 
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To obtain the effective frequency based on ship velocity and heading 
relative to the waves. The ship motion spectrum can then be 
obtained as follows. 


4 > (w ) = S (tu ) -r — RAO . ( w ) 
ii e w e 6u e 
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(A30) 


A computer program, SHPREF . FOR, was developed to accept a table of 
RAO values corresponding to a set of values, and then store the 
information in a data file. This table i3 then read by the program, 
and the above calculations are made. A numerical table and a plot 
of 4> versus w e is then output as the finished result; see Figures A2 
and A3. The plot was then sectioned by hand into six components and 
the amplitudes for these six frequencies were determined by the 
following relationship: 
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Figure A2: Example Output from SHPREF.FOR with Plot of 

Ship Motion S.oecturm for the DD962 in Heave 
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Figure A3: Example Output from SHPREF. FOR with Plot of 

Ship Motion Spectrum for the DD962 in Pitch 
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A component of the ship motion could then be calculated as: 

6 

i(t) = Z A C03 ( w t - 4 + E ) (A32) 

n e ii n 

n=1 n 

where 4^ is directly available from the ship data base information 
contained in Reference 8. E n is a random phase angle which is 
calculated from a random number generator with the output scaled to 
give values- between 0 and 6.242 radians. 

In the computer-run simulation the ship motion is calculated in 
the GENeral Ship Motion, GENSM.FOR, subroutine. This subroutine 
use3 the- associated amplitudes and phase angles corresponding to a 
particular sea state as stored in the data file, GENSM.DAT. 

D. TURBULENCE MODEL. 

The turbulence model for the computer-run simulation is located 
in the subroutine TURB3.F0R. It consists of a random number 
generator and the following principle equations: 

. 1/2 
A^ =* w n A^ + cf^ {2w^) n (1/TC ) 

a r = a r + TC a r (A33) 

The bandwidth, w n , was obtained from page 10, Figure 15, of 
Reference 9. The sigma values were obtained from a strip-chart 
recording of the acceleration of the vertical axis for the AV8-B 
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during a fixed-base simulation, as presented in Figure A4. The 
(1/TC ** 0.5) term is a correction for the effects of using digital 
computation. The random number input is represented by n. Becuase 
this is a filter acting to shape white-noise, it is necessary to 
precycle it initially before inputlng values into the simulation. 
This is also accomplished as part of the subroutine. 

E. HOVER HEIGHT INPUT 

The subroutine IN4.F0R provides the values for the flight path 
command logic hover height. The sequence starts at an initial 
altitude and follows a shallow cosine path to a selected hover 
altitude as a function of time. The hover altitude is then 
maintained as a constant for the remainder of the run. It is also 
possible to configure the subroutine to change to a second hover 
height during the run, although this function was not used other 
than for initial testing of pilot lead and lag time constants. 
Variables which can be adjusted before running a simulation are, the 
initial altitude, AMPA + AMPB, the hover altitude, AMPA, and the 
rate of descent, through the frequency term of the cosine function, 
WNS. 


F. DATA PLOTTING 

The subroutine MPLT2 .FOR provides the instruction to the system 
library and DISSPLA software tb produce a plot of the desired data 
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Figure A4: 


Strip Chart Recording of Vertical Acceleration 
Due to Turbulence as Obtained from the AV-8A 
Fixed-Base Simulation Facilities and Used to 
Determine Magnitude of Turbulence Values for 
the Turbulence Modeling in the Kon-Piloted 
Simulation 
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at the end of program execution. Currently, plots of the T/W ratio 


\ 


used as a percentage of time, and statistical variations as a 
function of the number of runs made are output (for examples see 
Figures A5-A7). A plot of the aircraft flightpath and ship deck 
position time history can also be set up with small changes in 
program configuration (for examples see Figures 31 and 32). A 
tabular listing of the runs made for each te3t configuration is also 
output, an example is given in Figure A8. 

The following table is - a listing of variable names and their 
appropriate values which must be edited into the indicated programs 
when a change in sea state is made. 


Table A1 : Valued for Variable Which Must Be Edited 

into Program with Change in Sea State 


VARIABLES which are fvkctioh of sea state 

Subroutine 

Variable 

Mane 

0 

Sea State 
4 5 

6 

TRAHO.FOR 

VH 

2.70 

6. 15 

7.79 

7.79 


SIGMA 

0.001 

0.007 

0.01 

0.01 

IN4.F0R 

AMPA 

0.0 

5.0 

'7.0 

9.0 


AMP 3 

40.0 

35.0 

33.0 

31.0 

VP A1C 

ABRTHT 

0.0 

1.5 

3.0 

6.0 

CENSM.DAT 
(Version No. 

- 

- 

;4 

;5 

;6 

VASCON.DAT Supplies the input data for Che test conditions 
and is also changed as apporpriate. 
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Tec- 0.1500000 

ELO 

0.7000000 7U* 

1.010000 

ZU- 0.1000000 

SEP S.- 

6 KP- 

7.0000002E-02 




RLM 

1. 1. 

F . T . 

TCVEL 

TIVRD 

TPPEfih 

THi-flX 

THMIh 

! 

2569.7 

21.2 

6.7 

7.9 

0 . S 68 

1.010 

0.877 


121.4 

24.1 

11.4 

11 .6 

0 .986 

1.010 

0.950 

3 

2026.6 

24.2 

8.8 

C ! 

« • i 

0.991 

1 .010 

0.950 

■5 

384.2 

23.4 

2. *3 

e.s 

0.985 

1.010 

0.950 

5 

2524.8 

19.7 

8.8 

10.8 

0.987 

1.010 

0.950 

6 

774.4 

17.6 

6.8 

G .4 

0.967 

1 .010 

0.977 

7 

2825.2 

27.0 

* *5 

8.8 

0.988 

l.oto 

0.950 

S 

262.6 

22.2 

5.1 

9.8 

0.936 

1 .010 

0.950 

9 

2729.3 

22.5 

6.7 

6.7 

0.985 

1.010 

0.950 

10 

2491.3 

22.8 

6.9 

9.7 

0.985 

1 .010 

0.950 

11 

2404.7 

1 1 . a 

5.6 

10.7 

0.987 

1.010 

0.977 

12 

945. 4 

24.9 

8.7 

8.3 

0.987 

1 .010 

0.950 

13 

1816.2 

IS . 4 

0.4 

11.1 

0.989 

1.010 

0.950 

M 

672.8 

22.3 

3.4 

6.5 

0.990 

1 .010 

0 . S 5 Q 

15 

3490.0 

24.6 

6.2 

9.9 

0.986 

1.010 

0.950 

16 

962.4 

16.9 

3.2 

8.6 

0.987 

1 .010 

0.877 

‘7 

3467.7 

19.2 

4.2 

6.5 

0.938 

1.010 

0.977 

.8 

543.1 

18.9 

10.9 

11.8 

0.986 

1.010 

0.950 

19 

3075.1 

24.9 

4.6 

8.3 

0.991 

1.010 

0.977 

20 

3225.8 

20.7 

10.2 

10.2 

0.985 

1 .010 

0.950 

21 

2802.3 

21.1 

8.2 

8.7 

0.984 

1 .010 

0.950 

22 

1931.2 

12.0 

3.1 

S.5 

0.984 

1 .010 

0.950 

23 

1901.3 

27.4 

11.0 

11.0 

0.983 

1 .010 

0.950 

2-1 

1974.7 

24.2 

4.8 

11.7 

0.991 

1 .010 

0.950 

25 

2514.5 

22.6 

3.5 

6 . 3 

0.989 

1.010 

0.977 

26 

2255.3 

12.3 

6.3 

8.0 

0.966 

1 .010 

0.977 

27 

3118.8 

19.9 

5.5 

7.6 

0.988 

1 .010 

0.977 

28 

1216.4 

30.6 

2.3 

8.0 

0.993 

1 .010 

0.950 

29 

1298.2 

19.1 

0.6 

7.4 

0.968 

1.010 

0.950 

30 

2666.2 

17.4 

9.6 

10.6 

0.957 

1.010 

0.977 

31 

2999.0 

24.4 

3.7 

7.3 

0.986 

1 .010 

0.950 

32 

1157.0 

23.2 

3 . 3 

5.0 

0.939 

1.010 

0.877 

33 

165.6 

12.4 

3.5 

8.5 

0.987 

1.010 

0.977 

3*1 

2190.3 

8.3 

8.6 

14.6 

0.990 

1 .010 

0.978 

35 

866.7 

20.2 

5 . 0 ' 

9.8 

0.984 

1.010 

0.950 

36 

3417.2 

28.2 

2.3 

8.9 

0.282 

1.010 

0.950 

37 

2984.7 

17.2 

5.2 

8.5 

0.888 

1 .010 

0.950 

38 

615.9 

16.1 

6.5 

S.O 

0.986 

1.010 

0.977 

39 

521.1 

25.8 

2.6 

12.1 

0.990 

1 .010 

0.S50 

40 

2124.9 

33.2 

10.3 

10.3 

0.991 

1 .CIO 

0.950 

41 

1051.6 

20.5 

3.1 

5.9 

0.888 

1 .010 

0.950 

42 

165.8 

12.2 

2.8 

9.5 

0.287 

1.010 

0.977 

1 

3314.5 

33.4 

10.9 

11.5 

0.991 

1.010 

0.850 

i4 

2679.8 

24.7 

0.7 

10.4 

0.989 

1 .010 

0.950 

45 

2616.7 

27.7 

8.6 

8.6 

0.987 

1 .010 

0.950 

46 

2302.2 

25.9 

1.4 

10.6 

0.989 

1 .010 

0.950 

47 

2956.3 

21 .2 

7.7 

7.7 

0.984 

1 .010 

0.950 


Figure A8: Example Output of Run Conditions and Results 

for the Hon-Piloted Simulation 
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48 

■' 1083.0 

16.0 

3.4 

8.6 

0.986 

1 .010 

0.S77 

49 

1784.8 

27.9 

4.2 

6.8 

0.981 

1.010 

0.950 

50 

2304.6 

22.7 

1.1 

10.4 

0.988 

1.010 

0.950 


MEAN<SEC) : 

21.55800 

FT 

SIG(SEC) : 

5.464322 



FTMflX: 33.40000 FTMItS : 8.300000 

TCVEL MEAh<FT/S): 5.526001 TDVEL SIG(FT/S): 3.060933 

TDVELMAX: 11.40000 TOVELMIh: 0.4000000 

riTCUT MEAN: 0.9879200 MTCUT SIG: 2 . 3633221E-03 

MTCUTMAX: 0.9940000 MTGUTM1N: 0.9840000 

TIME: 11:11 :54 DATE: 4-JAN-85 
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Figure A8, continued: Example Output of Run Conditions and Results 

for the Non-Piloted Simulation 
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G. LISTING OF NON-PILOTED SIMULATION COMPUTER 
VARIABLE DEFINITIONS AND PROGRAM LISTINGS. 



LIST OF VARIABLES FOR VPAlC.FOR PROGRAM 


NAME DESRIPTION 


AO CONSTANT FOR THE ZEROTH ORDER NUMERATOR TERM OF THE PILOT 

TRANSFER FUNCTION. 

A1 CONSTANT FOR THE FIRST ORDER NUMERATOR TERM OF THE PILOT 

TRANSFER FUNCTION. 

AAO CONSTANT FOR THE ZEROTH ORDER NUMERATOR TERM OF THE AIR- 
CRAFT TRANSFER FUNCTION. 

AAOP MAGNITUDE OF THE PILOT OUTPUT SIGNAL ATTRIBUTED TO THE 
NUMERATOR OF THE PILOT TRANSFER FUNCTION. 

AA> CONSTANT FOR THE FIRST ORDER NUMERATOR TERM OF THE AIR- 

CRAFT TRANSFER FUNCTION. 

AA2 CONSTANT FOR THE SECOND ORDER NUMERATOR TERM OF THE AIR- 

CRAFT TRANSFER FUNCTION. 

AA3 CONSTANT FOR THE THIRD ORDER NUMERATOR TERM OF THE AIR- 

CRAFT TRANSFER FUNCTION. 

ABRTHT ALTITUDE ABOVE THE DECK MOTION MEAN AT WHICH THE PILOT 
IS TO ABORT 'CHASING' THE SHIP DECK AND RETURN TO THE' 
ASSIGNED HOVER ALTITUDE. 

B1 CONSTANT FOR THE FIRST ORDER DENOMINATOR TERM OF THE 

PILOT TRANSFER FUNCTION. 

B2 CONSTANT FOR THE SECOND ORDER DENOMINATOR TERM OF THE 

PILOT TRANSFER FUNCTION. SET = 0.0 IN CURRENT PROGRAM 
BECAUSE NOT IN CURRENT USE. 

BA1 CONSTANT FOR THE FIRST ORDER DENOMINATOR TERM OF THE 

AIRCRAFT TRANSFER FUNCTION. 

BA1P MAGNITUDE OF THE PILOT OUTPUT SIGNAL ATTRIBUTED TO THE 
DENOMINATOR TERM OF THE PILOT TRANSFER FUNCTION. 

BA2 CONSTANT FOR THE SECOND ORDER DENOMINATOR TERM OF THE 
AIRCRAFT TRANSFER FUNCTION. 

BA3 CONSTANT FOR THE THIRD ORDER DENOMINATOR TERM OF THE 
AIRCRAFT TRANSFER FUNCTION. 
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BEG INSTRUCTION TO PLOT3.FC>R SUBROUTINE INDICATING 

BEGINNING POINT OF THE GRAPH X-AXIS. 

BUF DIMENSIONAL DUMMY VARIABLE USED IN THE TIME SUBROUTINE. 

BUFF DIMENSIONAL DUMMY VARIAJ3LE USED IN THE DATE SUBROUTINE. 

Cl FIRST RECURSION CONSTANT USED IN A.I.L. CALCULATION OF 

PILOT MODEL OUTPUT. 

CIA FIRST RECURSION CONSTANT USED IN A.I.L. CALCULATION OF 
AIRCRAFT MODEL. OUTPUT. 

C2 REMNANT. FROM EARLIER PROGRAM. NOT USED IN THIS VERSION. 

C2A SECOND RECURSION CONSTANT USED IN A.I.L. CALCULATION. 

SAME VARIABLE NAME IS USED FOR BOTH PILOT AND A/C MODELS. 

C2B SECOND RECURSION CONSTANT USED IN A.I.L. CALCULATION OF 
PILOT MODEL OUTPUT USED AS INITIAL CONDITION INPUT TO 
PILOT MODEL. 

C3A THIRD RECURSION CONSTANT USED IN A.I.L. CALCULATION. 

SAME VARIABLE NAME IS USED FOR BOTH PILOT AND A/C MODELS. 

C3B THIRD RECURSION CONSTANT USED IN A.I.L. CALCULATION OF 

PILOT MODEL OUTPUT USED AS INITIAL CONDITION INPUT TO 
PILOT MODEL. 

C4A FOURTH RECURSION CONSTANT USED IN A.I.L. CALCULATION. 

SAME VARIABLE NAME IS USED FOR BOTH PILOT AND A/C MODELS. 

C4B FOURTH RECURSION CONSTANT USED IN A.I.L. CALCULATION OF 

PILOT MODEL OUTPUT USED AS INITIAL CONDITION INPUT TO 
PILOT MODEL. 

C5A, C5B, C6A, C6B, C7A,C7B,C8A,C8B ARE SIMIALAR TO THE ABOVE. 

CKPB REMNANT FROM EARLIER PROGRAM. NOT USED IN THIS VERSION. 

CTNH ’CHOP THROTTLE NOW HEIGHT'? HEIGHT ABOVE THE DECK. AT 

WHICH THE PILOT IS TO QUICKLY REDUCE THRUST FOR LANDING. 

D1 REMNANT FROM EARLIER PROGRAM. NOT USED IN THIS VERSION. 

D2 REMNANT FROM EARLIER PROGRAM. NOT USED IN THIS VERSION. 

D3 AN OFTEN. USED COMBINATION OF OTHER VARIABLES. 


i 
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DAMP MAGNITUDE OF VERTICAL VELOCITY DAMPING THROUGH THRUST 
TERM. 

OH PILOT PRECIEVED ERROR IN HEIGHT (COMMANDED HEIGHT - 

ACTUAL HEIGHT + NOISE). 

DHD PILOT PRECIEVED ERROR IN VERTICAL VELOCITY. 

DHDD PILOT PRECIEVED ERROR IN VERTICAL ACCELERATION. 

DHDTR TRUE ERROR IN VERTICAL VELOCITY. 

DHTR TRUE ERROR IN HEIGHT. 

ELC ENGINE LAG TIME OONSTANT. 

EN RANDOM PHASE ANGLE GENERATED AND USED IN THE GENSM.FOR 

SUBROUTINE. 

ERR ERROR TERM (RANDOM NOISE WITH MAGNITUDE FROM -0.5 TO 
0.5). 

ERROR ERROR TERM (RANDOM NOISE WITH MAGNITUDE FROM 0 TO 1 ) . 

FLAG 2 REMNANT FROM EARLIER PROGRAM. NOT USED IN THIS VERSION. 

FLT FLIGHT TIME (INITIAL TIME MINUS TIME AT TOUCHDOWN). 

FLT1 PREVIOUS FLIGHT TIME VALUE, HELD FOR USE IN STATISTIC 

CALCULATIONS. 

FPLS ' FLIGHT PATH LOGIC SLOT* ; INTEGER USED TO SHOW SECTION OF 
FLIGHT PATH COMMAND LOGIC IN USE WHEN MONITERING 
SIMULATION. 

FPLS1 VARIABLE USED IN FINAL PRINTOUT TO INDICATE IF THE FLIGHT 
PATH COMMAND LOGIC HAD CYCLED THROUGH THE 'CHOP THROTTLE' 
SEQUENCE AT LEAST ONCE. 

GAP VARIABLE WHICH DETERMINES THE AMPLITUDE FOR A COSINE 

SMOOTHING FUNCTION IN THE FLIGHT PATH COMMAND LOGIC. 

GAPP VARIABLE WHICH DETERMINES THE AMPLITUDE FOR A COSINE 
SMOOTHING FUNCTION IN THE FLIGHT PATH COMMAND LOGIC. 

GD DUMMY VARIABLE USED TO TRANSFER THE VALUE OF GHDOT 

BETWEEN VASCON.FOR AND VA1C.FOR. 

GHDOT GAIN USED IN CALUCULATING THE VERTICAL VELOCITY DAMP- 
ING THROUGH THRUST VALUE. 
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GPP VARIABLE WHICH DETERMINES THE AMPLITUDE FOR A COSINE 
SMOOTHING FUNCTION IN THE FLIGHT PATH COMMAND LOGIC. 

HLIM VARIABLE WHICH CAN BE USED TO SET A LIMIT ON THE LOWER 
BOUNDS OF COMMANDED THRUST. 

HTA THE NUMERATOR PORTION OF THE PILOT MODEL OUTPUT USED AS 
THE INPUT CONSTANT FOR THE ZEROTH TERM IN THE NUMERATOR 
OF THE AIRCRAFT TRANSFER FUNCTION AFTER ADDITION OF 
NOISE AND LIMITS. 

HTB THE DENOMINATOR PORTION OF THE PILOT MODEL OUTPUT USED 

AS THE INPUT CONSTANT FOR THE ZEROTH TERM IN THE DENOM- 
INATOR OF THE AIRCRAFT TRANSFER FUNCTION AFTER ADDITON OF 
NOISE AND LIMITS. 

HTPO MAGNITUDE OF THE PILOT MODEL OUTPUT BEFORE THE ADDITION 
OF NOISE OR LIMITS. 


HTPOA FIRST PORTION OF EQUATION USED TO CALCULATE HTPO. 

HTPOB SECOND PORTION OF EQUATION USED TO CALCULATE .HTPO. 

I DUMMY VARIABLE FOR DO STATEMENT. 

IDO DUMMY VARIABLE FOR DO STATEMENT. 

IFLAG10 LOGIC SWITCH USED TO CONTINUE 'CHOPPED THROTTLE' 

SEQUENCE IN THE FLIGHT PATH COMMAND LOGIC AFTER ITS 
INITIALIZATION. 

IFLAG11 LOGIC SWITCH TURNED ON IN THE 'RUN FAST' SEQUENCE OF THE 
FLIGHT PATH COMMAND LOGIC AND USED TO DETERMINE THE 
INITIAL AMPLITUDE FOR THE COSINE SMOOTHING FUNCTION IN 
THE 'RUN' SEQUENCE OF THE FLIGHT PATH COMMAND LOGIC. 

I FLAG 1 2 LOGIC SWITCH TURNED ON IN THE 'RUN' SEQUENCE OF THE 

FLIGHT PATH LOGIC AFTER INITIALIZATION OF TIME CONSTANT 
FOR THE COSINE SMOOTHING FUNCTION, TO PREVENT REINITIAL- 
IZATION ON CONSECUTIVE PASSES. 

I FLAG 1 3 LOGIC SWITCH TURNED ON IN THE 'ABORT TO 
HOVER ALTITUDE' SEQUENCE OF THE 

FLIGHT PATH LOGIC AFTER INITIALIZATION OF TIME CONSTANT 
FOR THE COSINE SMOOTHING FUNCTION, TO PREVENT REINITIAL- 
IZATION ON CONSECUTIVE PASSES. 
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I FLAG 3 LOGIC SWITCH TURNED ON IN THE 'ABORT TO HOVER ALTITUDE' 

SEQUENCE OF THE FLIGHT PATH COMMAND LOGIC AFTER- INITIALI- 
ZATION OF A TIME CONSTANT FOR THE COSINE SMOOTHING FUNC- 
TION TO PREVENT REINITIALIZATION ON CONSECUTIVE PASSES. 

IFLAG5 DATA OUTPUT SWITCH, IF IFLAG5=1 THEN THE STATISTICAL 
INFORMATION IS CALCULATED AND OUTPUT. 

IFLAG7 LOGIC SWITCH TURNED ON IN THE 'CHOP THROTTLE SEQUENCE' 

OF THE FLIGHT PATH COMMAND LOGIC AND USED IN THE 'CHASE' 
SEQUENCE OF THE FLIGHT PATH COMMAND LOGIC TO DETERMINE 
THE INITIAL. AMPLITUDE FOR THE COSINE SMOOTHING FUNCTION. 



IFLAG8 LOGIC SWITCH TURNED ON IN THE IN4.F0R SUBROUTINE AFTER 

INITIAL COMMANDED LETDOWN FLIGHTPATH REACHES THE CONSTANT 
HOVER ALTITUDE VALUE, PREVENTS REINITIALIZATION OF LETDOWN 
SEQUENCE DURING A RUN. 

IFLAG9 LOGIC SWITCH TURNED ON IN THE 'CHOP THROTTLE' SEQUENCE 
OP THE FLIGHT PATH COMMAND LOGIC, WHICH CAN BE USED TO 
CHANGE THE THRUST LOWER LIMITS FROM THE NOMINAL VALUE 
USED FOR THE REST OF FLIGHT PATH COMMAND SEQUENCES. 

INT DUMMY VARIABLE USED IN DO LOOP. 

ISEAS INTEGER VALUE OF SEA STATE BEING SIMULATED. IT IS 

PASSED TO VARIOUS SUBROUTINES AS A SIMULATION PARAMETER 
AND FOR PARAMETER PRINTOUT FOR THE SIMULATION ON DATA 
OUTPUT. 

IT DUMMY VARIABLE USED IN VARIABLE DIMENSION STATEMENTS. 

IT1 TIME INCREMENT USED' IN THE A.I.L. CALCULATIONS 

IT2 THROUGH IT8 ARE TIME INCREMENTS USED IN THE A.I.L. CALC- 
ULATIONS. IT2=IT1**2/2l , IT3=IT1 **3/31 , ETC. 

IT9 THROUGH IT1 2 ARE REMNANTS FROM EARLIER PROGRAMING. NOT USED 
IN THIS VERSION. 

ITOT TOTAL NUMBER OF ENTRIES USED IN DETERMING A HISTOGRAM OF 
STATISTICAL DATA. 

J INITIAL SEED FOR USE IN THE RANDOM HUMBER GENERATING 

SUBROUTINE. 

JAY DUMMY VARIABLE FOR PASSING THE NUMBER OF GROUPS BEING 

RUN IN A GIVEN SIMULATION SESSION TO VARIOUS SUBROUTINES. 
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JRAN INITIAL SEED FOR USE IN THE RANDOM NUMBER GENERATING 
SUBROUTINE. 

K NUMBER OF RUNS/GROUP OF A SIMULATION SESSION. 

KAY DUMMY VARIABLE, NUMBER OF RUNS/GROUP OF A SIMULATION 

SESSION. 

KDOO INCREMENTED VARIABLE USED AS A DIMENSION IN A DO LOOP. 

KMEAN NUMBER OF ENTRIES IN A STATISTICAL MEAN CALCULATION. 

KMN NUMBER OF ENTRIES IN A STATISTICAL MEAN CALCULATION. 

KP NOMINAL PILOT GAIN. 

KPA MODIFIED PILOT GAIN. 

KPB MODIFIED PILOT GAIN. 

KPC MODIFIED PILOT GAIN. 

L INTEGER VALUE OF THE PILOTS PURE LAG TIME MULTIPLIED BY 

100 . 

LIM CALCULATED VARIABLE WHICH MODIFIES THE dT/W VALUE OUTPUT 

FROM THE PILOT MODEL TO KEEP IT IN PHYSICALLY REALIZABLE 
LIMITS. 

MIVRD VALUE OF THE MAXIMUM INSTANTANEOUS VELOCITY RELATIVE TO 
DECK ENCOUNTERED DURING A SIMULATION RUN. 

MTOUT MEAN THRUST/WEIGHT RATIO OUTPUT DURING A SIMULATION 
RUN. 

POUT VALUE OF THE PILOT MODEL OUTPUT, AFTER ADDITION OF NOISE 
AND LIMITS. 

PRECYC NUMBER OF CYCLES THE SUBROUTINE TRANO. FOR IS TO RUN 

THROUGH TO OBTAIN STEADY STATE OUTPUT BEFORE VALUES FOR 
NOISE AND TURBULENCE INPUTS ARE RETURNED TO THE MAIN 
PROGRAM. 

RNP VALUE OF NOISE CALCULATED IN TRANO. FOR INSERTED INTO THE 

PILOT MODEL AS INTERNAL PILOT NOISE. 

RNQ VALUE OF NOISE CALCULATED IN TFANO. FOR ADDED TO THE 

POSITON ERROR (PILOT PRECEPTION NOISE). 

SI COMMANDED FLIGHT PATH ALTITUDE. 
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> -1 

- ... - • - 1*4 


S2 

COMMANDED FLIGHT PATH VELOCITY. 


S3 

COMMANDED FLIGHT PATH ACCELERATION. 


SA1 

COMMANDED HOVER ALTITUDE. 


SA2 

COMMANDED HOVER VELOCITY. 


SA3 

COMMANDED HOVER ACCELERATION. 


SGAPP 

INITIAL POINT FOR A COSINE SMOOTHING FUNCTION 
'RUN' FAST' SEQUENCE OF THE FLIGHT PATH LOGIC. 

IN THE 

SGP 

INITIAL POINT FOR A COSINE SMOOTHING FUNCTION 
'CHASE' SEQUENCE OF THE FLIGHT PATH LOGIC. 

IN THE 

SGPP 

INITIAL POINT FOR A COSINE SMOOTHING FUNCTION 
'RUN' SEQUENCE OF THE FLIGHT PATH LOGIC. 

IN THE 

SSI 

SHIP DECK POSITION.- 


SSI z 

RELATIVE DISTANCE BETWEEN SHIP DECK AND A/C DELAYED 
BY THE VALUE OF THE PURE PILOT LAG TIME. 

SS2 

SHIP DECK VELOCITY. 


SS2ZD 

RELATIVE VELOCITY OF THE SHIP DECK AND A/C DELAYED 
BY THE VALUE OF THE PURE PILOT LAG TIME. 

SS3 

SHIP DECK ACCELERATION. 


SSAMP 

AMPLITUDE OF THE COSINE SMOOTHING FUNCTION IN 

THE 


'ABORT TO HOVER ALTITUDE' SEQUENCE OF THE FLIGHT PATH 
COMMAND LOGIC. 


SSGAP INITIAL POSITION FOR THE START OF THE COSINE SMOOTHING 

FUNCTION IN THE 'ABORT TO HOVER ALTITUDE' SEQUENCE OF THE 
FLIGHT PATH COMMAND LOGIC. 

SX DUMMY VARIABLE IN CALL TO INPUT1 . FOR, NOT USED IN THIS 

VERSION. 

T TIME (SECONDS). 

T1 MODIFIED VALUE OF T2, NOT USED IN THIS VERSION. 

T2 REMNANT FROM EARLIER PROGRAM, NOT USED IN THIS VERSION. 

T4 MODIFIED VALUE OF T2, NOT USED IN THIS VERSION. 
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TC TIME INCREMENT FOR RUNNING THE SIMULATION. 

TE THROUGH TE4 TIME INITIALIZED AT THE BEGINING OF A SET OF 

CONSECUTIVE RUNS THROUGH ONE OF THE FLIGHT PATH COMMAND 
LOGIC SEQUENCES FOR TIMING OF THE COSINE SMOOTHING . 
FUNCTIONS. 

TF FINAL TIME TOR WHICH THE SIMULATION ABORTS A RUN IF 

TOUCHDOWN HAS NOT BEEN ACHEIVED. 

TFIN DUMMY VARIABLE USED IN GRAPHICS SUBROUTINE TO INDICATE 
X-AXIS MAXIMUM VALUE FOR GRAPH. 

TG PILOT LAG TIME CONSTANT. 

TGEE VALUE CALCULATED IN TRANO. FOR FOR INPUT AS TURBULENCE 
INTO THE A/C TRANSFER FUNCTION. 

Til THE INITIAL TIME A RUN WAS STARTED. AT AS OUTPUTTED FROM 

A RANDOM NUMBER GEN RATION SEQUENCE WITH BOUNDS FROM 0 
TO 3600 SECONDS. 

TL PILOT LEAD TIME CONSTANT. 

TOUT THRUST-TO- WEIGHT RATIO. 

TOUTMAX MAXIMUM THRUST- TO-WEIGHT RATIO COMMANDED BY PILOT (WITH 
A/C LIMITS) DURING A GIVEN RUN. 

TOUTMIN MINIMUM THRU ST-TO- W EIG HT RATIO COMMANDED BY PILOT (WITH 



A/C LIMITS) DURING A GIVEN RUN. 


TSDA 

PART OF THE SEQUENTIAL CALCULATION OF STANDARD 
FOR FLIGHT TIMES IN A GROUP OF RUNS. 

DEVIATION 

TW 

THE MAXIMUM THRUST- TO-WEIGHT RATIO ALLOWED FOR 
OF RUNS. 

A GROUP 

VSDA 

PART OF THE SEQUENTIAL CALCULATION OF STANDARD 
FOR TOUCH DOWN VELOCITIES IN A GROUP OF RUNS. 

DEVIATION 

VTD 

RELATIVE VELOCITY OF SHIP AND A/C AT TOUCHDOWN. 


VTD1 

TOUCHDOWN VELOCITY FOR PREVIOUS RUN USED IN STATISTICAL 
CALCULATIONS. 

VTZ 

DUMMY VARIABLE PASSED TO SUBROUTINE TRANO. FOR. 
THIS VERSION. 

NOT USED 
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YO AN INITIAL CONDITION OF POSITION FOR THE A. I. L. CALCULA- 

TION OF THE. DIFFERENTIAL EQUATION DESCRIBING THE PILOT 
TRANSFER FUNCTION. 

Y1 POSITION OUTPUT OF THE PILOT TRANSFER FUNCTION. 

Y1A FIRST PART OF POSITION CALCULATION FOR OUTPUT OF PILOT 

TRANSFER FUNCTION. 

Y1B SECOND PART OF POSITION CALCULATION FOR OUTPUT OF PILOT 
TRANSFER FUNCTION. 

Y2 VELOCITY OUTPUT OF THE PILOT TRANSFER FUNCTION. 

Y2A FIRST PART OF VELOCITY CALCULATION FOR OUTPUT OF PILOT 

TRANSFER FUNCTION. 

Y2B SECOND PART OF VELOCITY CALCULATION FOR OUTPUT OF PILOT 
TRANSFER FUNCTION. 

Y3 ACCELERATION OUTPUT OF THE PILOT TRANSFER FUNCTION. 

Y3A FIRST PART OF ACCELERATION CALCULATION FOR OUTPUT OF 
PILOT TRANSFER FUNCTION. 

Y3B SECOND PART OF ACCELERATION CALCULATION FOR OUTPUT OF 
PILOT TRANSFER FUNCTION. 

YA1 POSITION OUTPUT OF THE PILOT TRANSFER FUNCTION AS USED 

FOR INPUT TO A/C TRANSFER FUNCTION. 

YA1 A FIRST PART OF POSITION CALCULATION FOR OUTPUT OF THE 

PILOT TRANSFER FUNCTION AS USED FOR INPUT TO A/C TRANS- 
FER FUNCTION. 

YA1B SECOND PART OF POSITION CALCULATION FOR OUTPUT OF THE 

PILOT TRANSFER FUNCTION AS USED FOR INPUT TO A/C TRANS- 
FER FUNCTION. 

YA2 VELOCITY OUTPUT OF THE PILOT TRANSFER FUNCTION AS USED 

FOR INPUT TO A/C TRANSFER FUNCTION. 

YA : 3 ACCELERATION OUTPUT OF THE PILOT TRANSFER FUNCTION AS 
USED FOR INPUT TO THE A/C TRANSFER FUNCTION. 

YAA1 POSITION OUTPUT OF THE PILOT TRANSFER FUNCTION USED FOR 
MONITERING PILOT TRANSFER FUNCTION OUTPUT. 

YAA1A FIRST PART OF POSITION OUTPUT OF PILOT TRANSFER FUNCTION 
CALCULATION. 
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YAA1B SECOND PART OF POSITION OUTPUT OF PILOT TRANSFER 
FUNCTION CALCULATION. 

YAA2 VELOCITY OUTPUT OF THE PILOT TRANSFER FUNCTION USED FOR 
MONITERING PILOT TRANSFER FUNCTION OUTPUT. 

YAA3 ACCELERATION OUTPUT OF THE PILOT TRANSFER FUNCTION USED 
FOR MONITERING PILOT TRANSFER FUNCTION OUTPUT. 

YDO INITIAL CONSTANT FOR VELOCITY USED IN THE A.I.L. CALCULA- 

TION OF THE DIFFERIENTIAL EQUATION REPERSENTING THE 
PILOT TRANSFER FUNCTION. 

Z AIRCRAFT HEIGHT ABOVE MEAN DECK POSITION. 

ZA1 FIRST PART OF THE POSITON CALCULATION FOR A/C TRANSFER 

FUNCTION OUTPUT.- 

ZA2 SECOND PART OF THE POSITION CALCULATION FOR A/C TRANSFER 

FUNCT/ON OUTPUT. 

ZB1 FIRST PART OF THE VELOCITY CALCULATION FOR A/C TRANSFER 

FUNCTION OUTPUT. 

ZB2 SECOND PART OF THE VELOCITY CALCULATION FOR A/C TRANSFER 

FUNCTION OUTPUT. 

ZC1 FIRST PART OF THE ACCELERATION CALCULATION FOR A/C TRANS- 

FER FUNCTION OUTPUT. 

ZC2 SECOND PART OF THE ACCELERATION CALCULATION FOR A/C 

TRANSFER FUNCTION OUTPUT. 

ZD A/C VELOCITY OUTPUT. 

ZDD A/C ACCELERATION OUTPUT. 

ZW VERTICAL VELOCITY DAMPING DUE TO AIRFRAME TIME CONSTANT. 

ZWDTDT NUMERATOR CONSTANT FOR THE A/C TRANSFER FUNCTION. 
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LIST OF VARIABLES FOR IN4.FOR 


NAME DESCRIPTION 


A THE INITIAL TIME FOR THE 'A' SEQUENCE OF THE HOVER ALT- 

ITUDE COMMAND LOGIC. 

AMPA FINAL HOVER COMMAND HEIGHT. ALSO THE AMPLITUDE FOR THE 
'A' SEQUNCE OF THE HOVER COMMAND HEIGHT LOGIC (THE 
COSINE FUNCTION OF 'A' SEQUENCE IS NOT USED IN THIS 
VERSION). 

AMPB THE AMPLITUDE FOR THE COSINE FUNCTION OF THE 'B' SEQUENCE' 
OF THE HOVER HEIGHT LOGIC. AMPA+AMPB GIVES THE INITIAL 
STARTING HEIGHT ABOVE THE' SHIP DECK MEAN. IN THIS VERSION. 

B THE INITIAL TIME FOR THE *B' SEQUENCE OF THE HOVER ALT- 

ITUDE COMMAND LOGIC. 

I FLAGS LOGIC SWITCH TURNED ON AFTER THE INITIAL COMMANDED LETDOWN 
FLIGHTPATH ( ' B' SEQUENCE) REACHES THE CONSTANT HOVER ALT- 
ITUDE VALUE (AMPA) TO PREVENT REINITIALIZATION OF THE 'B' 
SEQUENCE DURING A GIVEN RUN. THIS VARIABLE CORRESPONDS 
TO IFLAG8 IN THE IN4 . FOR CALL STATEMENT FROM VPAIC.FOR. 

ISFLAG1 LOGIC SWITCH TURNED ON AFTER THE TIME INITIALIZATION IN 
THE 'A* SEQUENCE TO PREVENT REINITIALIZATION OF THE TIME 
DURING A GIVEN RUN. 

51 COMMANDED HOVER FLIGHT PATH ALTITUDE. 

52 COMMANDED HOVER FLIGHT PATH VELOCITY. 

53 COMMANDED HOVER FLIGHT PATH ACCELERATION. 

SX DUMMY VARIABLE USED IN CALL STATEMENT TO IN4.FOR, NOT 

USED IN THIS VERSION. 

T TIME (SECONDS). 

TAU TIME MINUS THE TIME OF INITIALIZATION OF EITHER THE 'A* OR 

'B' SEQUENCE. 

TI THE INITIAL TIME A RUN WAS STARTED AT, AS OUTPUTTED FROM 

A RANDOM NUMBER GENERATION SEQUENCE WITH BOUNDS FROM 0 
TO 3600 SECONDS. 
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VARIABLE USED IN THE INITIALIZING. THE TIME FOR THE 'A' 
SEQUENCE. 

WNS FREQUENCY OF THE COSINE FUNCTION IN THE ’A' SEQUENCE. 

WNS1 FREQUENCY OF THE COSINE FUNCTION IN THE 'B' SEQUENCE. 

LIST OF VARIABLES FOR TRANO.FOR 


NAME DESCRIPTION 

FLAG REMNANT FROM EARLIER PROGRAM, NOT USED IN THIS VERSION. 

I DUMMY VARIABLE USED IN DO LOOPS. 

IB I FIXED VALUE OF XB. 

IS IFIXED VALUE OF S. 

M THE NONINTEGER PART OF THE TIME (IN SECONDS) OF THE 

SYSTEM CLOCK. 

PRECYC NUMBER OF CYCLES THE TRANO.FOR SUBROUTINE IS TO RUN 

THROUGH TO OBTAIN STEADY STATE OUTPUT BEFORE VALUES FOR 
NOISE. AND TURBULENCE INPUTS ARE RETURNED TO. THE MAIN 
PROGRAM. 

RNO OUTPUT OF THE - FIRST ORDER FILTER USED AS INTERNAL PILOT 

; NOISE. 

RNOD FIRST DERIVATIVE OF RNO. 

RNP RNO MODIFIED BY A GAIN. 

RNQ OUTPUT OF THE FIRST ORDER FILTER USED AS PILOT PRECEP- 

TION ERROR. 

RNQD FIRST DERIVATIVE OF RNQ. 

RT SYSTEM TIME MINUS THE NONINTEGER PART OF THE SYSTEM 

TIME, IN SECONDS. 

S INDICATES IF XM IS EVEN OR ODD. 







; j 


SIGMA THE SIGMA VALUE FOR THE FILTER USED IH GENERATING RANDOM 
TURBULENCE. 

SIGMAN THE SIGMA VALUE FOR THE FILTER USED IN GENERATING THE 
'INTERNAL PILOT NOISE. 

SIGMAQ THE SIGMA VALUE FOR THE FILTER USED IN GENERATING THE 
PILOT P RECEPTION NOISE. 

TC THE PROGRAM TIME INCREMENT. ' 

I 

TGEE OUTPUT OF THE PkIST ORDER FILTER USED TO GENERATE j 

RANDOM TURP/JLENCE. j 

} 

TM SYSTEM TIME IN SECONDS. i 

j 

VTZ SAME AS TGEE. 1 

VTN FREQUENCY TERM FOR THE TURBULENCE SHAPING FILTER. I 

WNN FREQUENCY TERM FOR THE INTERNAL PILOT NOISE SHAPING 

FILTER. 

WNQ FREQUENCY TERM FOR THE PILOT PERCEPTION ERROR SHAPING. 

FILTER.- 

XB MODIFIED VALUE OF RT. 

XM FLOATED VALUE OF M. 

XS FLOATED VALUE OF IS. 

Y VALUE. OF 12 RANDOM NUMBERS ADDED TOGETHER IN THE PROCESS 

OF CREATING A GAUSSIAN DISTRIBUTION FROM A WHITE NOISE 

SOURQE. T v 


YA 

A 

RANDOM 

NUMBER 

GENERATED 

FROM 

A 

WHITE NOISE SOURCE. 

YB 

A 

RANDOM 

NUMBER 

GENERATED 

FROM 

A 

GAUSSIAN DISTRIBUTION 


SOURCE. 
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VARIABLES FOR GEN3M.FOR 


NAME 


DESCRIPTION 


EN 


RANDOM PHASE ANGLE. 


HI THROUGH H6 ARE THE CONTRIBUTIONS TO THE SHIP HEAVE IN FEET, 
FROM EACH OF THE SIX SINE CONPONENTS TO THE SHIP 
HEAVE MOTION APPROXIMATION. 

HD1 THROUGH HD6 ARE THE CONTRIBUTIONS TO THE SHIP HEAVE VELOCITY 
IN FEET/SEC, FROM EACH OF THE SIX SINE CONPONENTS TO 
THE SHIP HEAVE MOTION APPROXIMATION. 

HDD1 THROUGH HDDS ARE THE CONTRIBUTIONS TO THE SHIP HEAVE ACCEL- 
ERATION IN FEET/SEC**2, FROM EACH OF THE SIX SINE CON- 
PONENTS OF THE SIX SINE CONPONENTS TO THE SHIP HEAVE 
MOTION APPROXIMATION. « 


IB 

M 


IFIXED VALUE OF XB. 
I FIXED VALUE OF S. 


PI THROUGH P6 ARE THE CONTRIBUTIONS TO THE SHIP PITCH POSITION 
IN RADIANS, FROM EACH OF THE SIX SINE CONPONENTS TO 
SHIP PITCH MOTION APPROXIMATION. 

I 

PD1 THROUGH PD6 ARE THE CONTRIBUTIONS TO THE SHIP PITCH VELOCITY 
IN RADIANS/SEC, FROM EACH OF THE SIX SINE CONPONENTS. TO 
SHIP PITCH MOTION APPROXIMATION. 

PDD1 THROUGH PDD6 ARE THE CONTRIBUTIONS TO THE SHIP PITCH ACCEL- 
ERATION IN RADIANS/SEC**2, FROM EACH OF THE SIX SINE 
CONPONENTS TO SHIP PITCH MOTION APPROXIMATION. 




// 


RT 

S 

51 

52 
S2A 


SYSTEM TIME MINUS THE NONINTEGER PART OF THE SYSTEM 
TIME IN SECONDS. 

INDICATES IF XM IS EVEN OR ODD. 

SHIP DECK POSITION RELATIVE TO THE MEAN DECK POSITION. 
SHIP DECK VELOCITY. 

SUM OF THE SHIP HEAVE VELOCITY CONTRIBUTIONS TO THE 
SHIP HEAVE MOTION. 
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S 



S2B 

S3 

S3A 

S3B 

T 

T1 

TM 

XB 

XM 

XS 

Y 


\ / 



SUM OF THE SHIP PITCH VELOCITY CONTRIBUTIONS, MULTI- 
PLIED BY THE DISTANCE FROM THE SHIP C.G. TO THE SHIP 
LANDING PAD BULLSEYE TO OBTAIN THE CONTRIBUTION OF THE 
SHIP PITCHING MOTION TO HEAVE AT THE LANDING PAD BULLS- 
EYE. 

SHIP DECK ACCELERATION. 

SUM OF THE SHIP HEAVE ACCELERATION CONTRIBUTIONS TO THE 
SHIP HEAVE ACCELERATION. 

SUM OF THE SHIP PITCH ACCELERATION CONTRIBUTIONS, MULT- 
IPLIED BY THE DISTANCE FROM THE SHIP C.G. TO THE SHIP 
LANDING PAD BULLSEYE TO OBTAIN THE CONTRIBUTION OF THE 
SHIP PITCHING MOTION TO HEAVE AT THE LANDING PAD BULLS- 
EYE. 

TIME (SECONDS). 

THE INITIAL TIME A RUN WAS STARTED AT, AS OUTPUTTED 
FROM A RANDOM NUMBER GENERATION SEQUENCE WITH BOUNDS 
FROM 0 TO 3600 SECONDS. 

SYSTEM TIME IN SECONDS. 

MODIFIED VALUE OF RT. 

FLOATED VALUE OF M. 

FLOATED VALUE OF IS. 

RANDOM NUMBER WITH MAGNITUDE OF -0.5 TO 0.5. 
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r* * ********* ** ** * *** * * ** * * *** * 

* 

r * ** m* * ** * m** * * ** * ***** * ** 

* 

WmamW i ********** * 

a- 




\ 


2-FE8-1385 13:38:24.07 
2-FE8-19S5 13:38:24^07 
2-FE8-1S35 13:38:24.07 


FSOO:{ STEVENS. SHIPSTUFFJSHPREFLIN.CDMjl *** * ** ** * * 
FSDO: t STEVENS. SHI PSTUFF] SHPftEFUN.CGM;l ** * ** * * *** 
FSM:ISTEVENS.SHJPSTLfF]SHFF£FLIN.C0M5l ********** 



eS 


SHPREFL IN. COM ; 1 


i 

I 



i 

i 


i 
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\ 

2-FEB-1SS5 0:33:24.07 FSOO:ISTEVENS.SHIPSTlffF]SHPRFr IN.CDfljl *********** u-O' 


$ SET VERIFY 

$ ASSIGN/USER SYSSCOMMAND SYSSINPUT 

$ LINK SHPREF,DI SPLOT, SYS$LIBRARY : INTLIB/LI B , DI SSPLA/LIB, INTLIB/LI B 
* SET NOVERIFY 



( 




2-FEB-1985 13:23:47.97 FSOO : [ STEVENS . SH I PSTUr F ] SHFREF . LI S ; 1 
2-FEB-1985 13:23:47.97 FSDO: [ STEVD43 . SHI PSTlffT] SKPREF. LI S ; 1 
2-FE8-1S85 13:23:47.97 FSDO:[STEV£NS.SHIPSTirF]S«’REF.LIS;l 


SHPREF - l_I S ; 3. 


2-FEB-1985 13:23:4^.37 
2-FEB-1985 13:23:47.97 
2-FE3-198!' 13:23:47.97 


A40 SHIPSTLFF]SHPREF.L1S;1 
JHIPSTUFF JSHPREF.LISjl 
FSOO:[STEVEMS.SHIPSTIFFJSK?REF .LISjl 


/ 



i 


2-Frt-1385 13:21:49 UW-U FORTRAN V3.5-62 Pagt 1 

10-©ec-lS84 13:18:18 f SM :( STEVENS. SHI PSTUffJSHPJiJ.FM; 23 


0801 

C2345578301234367890...HC. 


0192 




0003 


DINBSI CH HIM) ,(£(301 ,1*1(30) ,34(30) 


0004 


REAL HUS 


0003 


OWACTW 12 run 


coot 

C 

VW2.17 142.17 !FT/$EC 


0097 

C 

HUS*2.0344 ! RAO 


0003 

C 

HS’12.0 !FT 


000) 


MINT*, 'DO YOU WISH TO OTTER DATA IN WJJ FILE? 1=YES' 


0010 


ACC£PT*,FLAGA 



0011 


PRINT*, 'DO YOU WISH A PLOT Cf THE DATA? 1*YES' 


i 

0012 


ACCEPT* ,F1AS 


j 

0013 


PRINT*, 'ENTER TO.VS.HUS.KS' 


! 

0014 


ACCEPT* , TO, VS, IRS ,HS 

i 

0013 

C 

T013.1 ! SETC 


S 

0018 




; 

0017 

C **t 

•ENTER FILEWS IN ACCEPT STATEHENT IN FORM Of 'films*' , 


i 

0013 

C 

i.*., IF YOU NAME IT STUFF. WI, INPUT IT AS 'STUFF .CAT' 



0013 


TYPE*, 'ENTER FILE WE' ! IffUT FILENAME AS 'filtftjs*' 


I 

0020 


ACCEPT*, FEW 


i 

0021 


IF (FLAGA .E0. 1) 60 TO 3 


0022 


OPEN (UIIT»l,FILE=fUH,TYPE*'CU') 


0023 


GO TO 10 


0024 

5 

CPEH 'lNIT»l,FILE=fU«,TYPE*'NEH') 


0023 

10 

J»M 


0025 


L*2C0 


0827 


K*500 


0028 


R S 1 


0023 

13 

DO 30 I*L,H,J 

I 


0030' 


H(R)*FUOT(I>/10C0. 

i 

• r 

0031 


IF (FIASA .E8. 1.0) GO TO 30 

| 


0932 

20 

READ (1,23) FAQ 

j 


0933 

25 

F0SmT(lX,FI0.7). 



0034 


GO TO 40 

i 


0035 

30 

TYPE*,'If>',H(R) 

5 


0035 


PRINT*, 'INPUT CCRRESPWDINO RAO UALUE' 

• 


0037 


ACCEPT*, RAO 

i 


0038 


WRITE (1,33) RAO 

1 


0033 

35 

FOSmT(lX,F10^) 

5 


0040 

40 

IE(R)=4J(R)-H(K)**2*VS*C0S(HUS)/32.2 


0041 


SiW(483.5*HS**2)/((T0**4)*(W(R)**3))) 


0042 


SWMXP(-1944.5/((H(K)*T0)**4)) 

» 

0043 


Sf(R)*S4Y*S43 

. j ■ 

0044 

C 

SH(K )=3tA*£3-£3 


0043 


ROKM*VS*COS(MUS)*HE(R)/32.2 

? 

0045 

C 

PRINT*, '4K444M4' ,CCHR 

» 

0047 


IF (RW .LT. 1.0) DO TO 43 


0043 


C»»E«1.0/((A8S(RC1K)-1)H0.5) 


0043 


GO TO 47 


0059 

45 

C«£*1.0/((1-RCHR)**0'.5) 


0031 

C47 

PRINT*, (SEWE 


0032 

47 

PHI(R)=S((K)*[S«NC*RAO 


0053 


IF (PKKR) .RT. PHimx) PHIKWHKR) 


0034 

C 

PRINT*, PHI (R),PHIMAX 


0033 


R«Rtl 


0935 

50 

CONTINUE 


0057 


IF (R .GT. 24) GO TO tt 
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sprefvwin 
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8053 


IF (R .0T. 15) 60 TO 55 

0053 


1*621 

0000 


K*750 

00(1 


J*25 

00(2 


GOTO 15 

00(3 

55 

L'fflO 

00(4 


jj-ipgo 

00(5 


3*50 

00(( 


GO TO 13 

00(7 

(0 

00 (5 1*1,24 

00(8 


PHI(1)*PKI(1)/PH]HW 

00(3 

65 

CONTINUE 

0070 


IF (FLAGS .15. 1.0) GO TO U 

0071 


CALL mRr.'lTfHE.PHl.It) 

0072 

£4 

PRINT*,' H HE PHI', 

0073 


J' SH' 

0074 


TYPE*,/ — ■ ' 

0073 


00 70 1*1,24 

0076 


TYPE*, N(I),H£(!),PH1<», 31(1) 

0077 

c 

TYPE*, 'WE*', WEI 1),' FHI«',ttfl(!),' SHU)*', 21(1) 

0078 

70 

CtWTINUE 

0073 


PRINT*, 'ramX*',FHIKW,' HUS*', HUS, ' IS*',H$,' TO*', TO 

COCO 


PRINT* ,FIW1 

0081 


acsEU) 

0082 


STOP 

0083 


DC 


PROGRAM SECTIONS 


Nas* 

Sylft 

Aur-ibsl** 



0 (CODE 

1208 

PIC CO> Fa Ltt 

SHR EXE 

R O IOST LC7IG 

1 4PCATA 

273 

PIC C£N Fa Ltt 

SKINQKE 

RD NCHST LOG 

2 HOWL 

732 

pic cm na Ltt lesm ncexe 

RD ; NRT LONG 


Tttal Spac* Allocate) 2233 


DITRY POINTS 

Mdwi Typt Njs* Wtrucn 

o-ooooaooa siprefitoin 

vwtiAais 


A&frm Type Nas» ftmitwtw Mtimcn 


2-00000229 R*4 805 

49* 

50* 

32 

2-OOOOOIFO R*4 FLAW 

10* 

a 

31 

2HJ0C001F4 R*4 FLAE3 

12* 

70 


2-OOOOOlEO CHW FUM 

3 

20* 

2 a 

2-00000200 R*4 HS' 

14* 

41 

73 


L 
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2-80800214 

1*4 

1 

29* 

30 

67* 

£8(2) 

73* 

7£(4) 



2-00000204 

1*4 

J 

23* 

23 

£1* 

£5* 





2-00000210' 

1*4 

8 

28* 

30 

35 

48(3) 

41 

42 

43 

43 




52(2) 

53(2) 

55(2)* 

57 

53 

71A 



2-00000203 

1*4 

L 

26* 

29 

59* 

£3* 





2-000002CC 

1*4 

M 

27* 

29 

CO* 

£4* 





2-OOOOOlFC 

8*4 

HUS 

4 

14* 

40 

45 

79 




2-00000220 

8*4 

PHIHtt 

53(2)* 

03 

79 






2-80000218 

R*4 

RAO 

32* 

37* 

33 

52 





2-00000224 

8*4 

ROW 

45* 

47 

43 

a 





2-0000021C 

8*4 

SKA 

41* 

43 







2-00000220 

8*4 

S3 

42* 

43 







2-OOOCOlFB 

R*4 

TO 

14* 

41 

42 

79 





2-000001FC 

8*4 

VS 

14* 

40 

45 







ARRAYS 










AdOrttt 

Typ* 

Nat 

Attributn 

Eytn DiMntion*- 

Rtfmncts 





2-OOOOOOFO 

8*4 

PHI 


120 (30) 

3 

52* 

53(2) 

68(2)* 

71A 






78 





2—00000163 

8*4 

34 


120 (30) 

3 

43* 

52 

76 


2-00800000 

R«4 

M 


120 (30) 

3 

30* 

35 

40(2) 

41 






42 

7 £ 




2-00000073 

8*4 



120 (30) 

3 

40* 

45 

71A 

76 


urns 


Addmi 

Libel 

(tfimcn 


0-00000106 

5 

21 

248 

8-00000*17 

10 

23 

258 

0-0000012A 

15 

298 

£2 

** 

20 

328 


1-00000108 

23” 

32 

338 

0-000081B8 

30 

31 

358 

1-00000)11 

33' 

33 

358 

0-00000236 

40 

34 

408 

8-00000263 

45 

47 

508 

C-80C002F1 

47 

49 

528 

H 

50 

29 

568 

0-00000336 

55 

53 

638 

0-00000369 

60 

57 

£71 

** 

63 

67 

698 

0-C00003SC 

66 

70 

728 

** 

78 

73 

788 
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r SffMTtfttlH 


Mrt- 1^3 Uiat«» uw-u f«ti!an v3.3-«2 p*s» 

ii-c^-ij 8< 13:19:1s rssoitsiatHS.SHiPsTirFiamr.fM;?! 


FlNCTItKS «* S96JOUIINES KJKENCFD 

Tyyy No* bfnnm 


n»»aosE 

81 


FOSMPOI 

22 

24 

KAROT.T 

71 


fm msecs 

44 

43 

»4 HIHKXP 

42 



t 

f j 

Kf TO FQTXDCE FIASS I ! 

■ • Vyltiy Hedifiyd I i 

♦ - Oyfioia? Ryfyrtnc* I * 

A - Actoil Ar$UE«t, pmibly codified I j 

D * Data loitialitttita I j 

(>) - Huo&yr tl tcttrrnctt »» liny I j 


(«) - Huo&yr ef tcarrfnctt »a liny 



CHWfi OUAUFIRS 



FtKTKK /C*0/US SWRFF.FO* 

/KES*(KW12rS,TlWaKa) 

/STA€OS»«(10Slia«(,MKaSE£ FOSSI) 

/sof»<Nors£p*oc£ssoa .noinouk^?) 

rm /ws.atATiia /u /optima minings /mbuiss /cubss tiraE«£ /rxfvamiEjxK /txNTiNi»Ti&5*i» 


CWttfTICN STATISTICS 

la Tis»: A, 37 syccitft 

tltju-l Tiay: 3. <7 itcondi 

Py$yFwl(i: 21$ 

Dynaic fta»ryi 130 p»$M 




» mmmu t nn +«ir inm i u « 

H 

WM i miiiiminmniHtH 

H 



** 


« 

w**44ww*m*wu * ww 


MfB-m 13i33i33.32 
2-fEMSS 13:39:33.32 
2-TQ-15« 13i35:J3.S 


r$8fli|STEV!}B.»ireturna«I»'IL.r«|3 »nmmmimimmu»i 
fSMiismtw.swpsnrntwsMKfiL.rcaii wmntww i)imnmn > 
fS60:lSTtVOfi.SH!fWf}«911ffll.rCR|3 MtHtHtt tt ******* ** **** 


QNSM I r^lF" I l_ . FOR 5 3 


2-rra-1535 13:39:33.32 rS«:iriIV0S.SHIPmfT!CN31IW’lL.r!«)3 mtHrtiHHmuHrtutntt 

■Km.tK* M.».» V) CUIKtirf MMUIUTM r «•.» 
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7 




CZ345S7S»Z345l?8J0...*tC. 

DWWCT0S*I2,FU« 

nujwc * 

MOTT*, FUJI 

OPO» (IWT *4,FIU"flWl t rifP£*'fS4' ) 

rm*,'iipvr wue fc* isms, tk sea state* 

ACOTT*,J 
WITE (4,5) J 
5 FGSKT (IX, 14) 

TYFEVMW WU£S FC* «,AI,{W2^TH, WO PHIO* i 

m*,'TTFE HKEXS .G£. 50 TO KIIT' i 

7 TYPE*,' IFfUTlS.AZ, PHIZ, AIK, FWO* ] 

«XOT*^£^Z,fflIZ,ATH,F«!0 1 

IF (AZ ,6£. 50) E3 TO 53 I 

WITE (4,10) WE^Z.PHIZ^TX.PHIO I 

1# FOBWT (1X,5(F8.4)) ; 

00 TO 7 J 

50 atE£(4) J 

stop ; 







* 

m i uiiii i iiimmummii 

* 

tWWWtMWWWmttH 

* 


2-FEB-iMS Hi57i54.i3 FSOOilSUVO«.»lftll»MVtoW*U.l.l>|» 

2-FQ-1S85 Ui57r54.13 rS53:[ST£VOB.SHIPST«T)WSI«-lL.LIS!l 

2-ftB-lSeS lli57:M.l3 r30i(STtVDiS.SH!PSHIT]y«ura.LlSil MUiniW t mmiM i MiHH 


UAS INFI L.UISjl 


tmmmHrtmHt+ttt+rttrt 


* 



2-FEB-1983 ll:57iM.13 
2-ftB'lJSS 11i57i54.13 


raO![STrVD6.SHIPSTlJFFIWIIflL.US!l I H l imHUHW l MMi i H tt 
rSOailSTtVtW.SHIKTirnVttSIMriUUSll WW W H tll t imuimm 
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2-F«iH93S 11:49:09 WX-11 F0PTR7W V3.5-62 p, ? , i 

10-0«c-19S4.16:23:10 FS00:[STEUmS.$mPSTlJFT!W$im.FC*;3 


o 


0001 

C234367890 234367890. . tc. 

0902 



0003 


PEAi rf 

0004 


CHNWaa*12,FlW1 

0005 


TYPC4/IWUT FILDm * 

0005 


accept*, fuh 

0007 


CPCN (miT*l.FILE*FU01,TYPE*'NQn 

0008 


TYPEt.'IffUT HAUJCS FC» J, Tit NU«£X OF CEMHTKNS, 

0009 


TYPE*, 'NO K, Tie NIMJ CF RLMS FCA EACH CXM51TICH, 

0010 


ACCEPT*, J.P 

0011 


MUTE (1,3) J,P 

0012 

5 

F0B14T (IX, 204)) 

0013 


TYPCVlim Wiles PC# LAG,ELC,TV.2H, tit) HP' 

0014 


dtevtype Ninons .ge. so to quit' 

0013 

7 

TYPEVIfFUT LA8.aC.TH.7H,W>' 

0015 


ACCEPT*,lA6,aC,TW,2>4.PP 

0017 


IF (LAG .GE. 50) GO TO 50 

0018 


MUTE 0.10) lAG.aC,TO,2H,M> 

0013 

10 

FORMAT OX,i4,4(F10.S)) 

0029 



0021 


GO TO 7 

0022 

S3 

CLCSE(l) 

0023 


STOP 

0024 


DC 


PPOGRArt- SECTIONS 


Nn* BytH Attribute 


* 

0 1CCOE 

446 

PIC COT pa ia 

SHU EXE 

KO POST ima 

t 

1 4P0ATA 

215 

pic cm pel La 

3R HXXE 

PO HOST L(HG 

* 

2 4L0CAL 

128 

pic cm sa La ncss noexe 

PO KPT Etna 

- 

TtUl Spjc t Allocitri 

790 





DTTPY POINTS 

Addrm Type Prftrtncts 

e-oooocooo wsinfilimain 


UNUAH.CS 

Addmt Typr ll» Attributu Peffrmcn 


2-ooooooie r* 4 ac 

16* 

18 


2-00000000 CHAP FUM 

4 

6* 

7A 

2-00000010 1*4 J 

10* 

11 


2-00005014 1*4 K 

10* 

11 


2-OOOOOOOC P*4 PP 

3 

16* 

18 

2-00000018 1*4 LAG 

16* 

17 

18 


C 


A48 


/ 

/ 


c 



WSIWILWW 2-feb-1933 11:49:09 UW-11 FCRTR/fi V3.5-62 • P*» 

10-D*c-»19S4 16:23:10 FSmSTMN$.StfPSTUFf)VASINFIl.FCR:9 


2-00000020 

R*4 TW 

16* 

18 

2-00000024 

R*4 ZM 

16* 

18 

LAsas 




Address 

Label 

References 


1-OOOOOOC1 

5* 

11 

120 

0-00000102 

7 

150 

21 

1-OOOOOOC8 

10' 

18 

13* 

0-000001AF 

50 

17 

22* 

FINCTHNS «e 

SU&KtWTIKS REFERENCED 



Type Nase 


References 



FORttLOSE 22 

FORWfEN 7 


REY TO REFERENCE FIA6S 
* - Value (tedif ied 

0 - Defining Reference 

A - Actual Arguaent, possibly codified 
0 * Data Initialization 

(n) - Number of eccurrtncts on lino 


CtmtO QUALIFIERS 

FORTRtfi aiS/CRO LKSINFJL.FW 

/OCCX=(N080lfes,O.tRfLOt,N(XKr£f.atH) 

/DEBOG»(HOSYmXS,TRACFEW3tl 
/ST«C«D=(NOSWTa(,NOSOUtCfJOSm 
/SKCU*t)mEPROC£SSOR .NOINCLUDE ,mf>) 

/F77 /NOGJLCATINO /M /OPTIMIZE /IWNINGS /NOD LINES /CROSS REFERENCE /WJttCHINE.CODE /CCKTINUATlCriS«19 


CtKPILATlON STATISTICS 


Run Time: 
Elapsed Ti«e: 
Pago Faults: 
Oynaaic Memory: 


1.68 seconds 
2.53 seconds 
133 

125 pages 


C 


A49 



I 



tttt tiHi i iUU t t i Utt tittilrt i i 

** 

« 

iU i ttiiti i kUiiiinkkiHtk ttt 

« 


2-FE8-l?$5 13:31:47.1? FS0a:(STEVD^.SH!m5FlGN^iNra.Ll$;i 

2-FFB-1935 13:31:47.17 FSW:(Sm’0«.SHlPSTlFF]GN91!HFIl..US;l 

2-FEB-1985 13:31:47.17 FSOflitSTWD<S.SHlPSnrFlGNS«lMFlL.USjl 


GN SM IMF" IL.LIS;! 


2 -HB- 15 S 5 13 : 31 : 47.17 FSCOtlSTEVFHS.SHlPSTtTFlWIIff Il.LIS;! 

«44U>«ttl>>Hl>ltllUHlW** *>.CrO-»fto* < 0 . 91,47 17 retsft.f CTpCVJC CUl C^T»CC>«MCWUC1 1 »«e.1 

A50 



umtHtwHiutuuium it 




i -r 


/ 


2-Feb-1983 13:21:30 UW-I1 FOaTRAN V3.5-62 Pl?e 

10-Dec-1934 13:51:03 FSOOUSTEVDiS.SHIPSTWTlMliNFfL.FOSiS 


0001 

C234567890234567899...ete. 

0002 



0003 


CH*ACTE**12,FU« 

0004 


TYPEVIWT FILENAME ' 

0005 


ACCEPT*, Flirt 

OOOC 


OPEN (tNIT»4 l FIL£*Flirt,TYPE»'NF l, '» 

0007 


TYPE*,'ltfVT VALUE FOR ISEAS, Ti£ SEA STATE' 

0009 


ACCEPT*, 3 

0009 


WRITE (4,3). 4 

0010 

3 

Fomtr (ix , id 

0011 


TYP£*,'II#W (BLUES FOR it, A2, PHIZ, ATM, WO 

0012 


TYPE*, 'TYPE MISERS .GE. 50 TO CUIT' 

0013 

7 

TYPEt.'IWW HE,AZ,PHIZ.ATH,PHIO' 

0014 


ACCEPT*, HE, AZ,PHIZ,ATH,PHIO 

0015 


IF (AZ .GE. 50) GO TO SO 

001S 


MUTE (4,10) HE,AZ,PHIZ,ATH,PH10 

0017 

10 

FORMAT (1X,3(F8.4)) 

0018 



0019 


GO TO 7 

0020 

50 

aOSE(4) 

0021 


STOP 

0022 


DO 


PROGRAM SECTIONS 


Htst 

Bytes 

Attributes 


0 tCOOE 

403 

PIC CON REL LCl SHR EXE 

RO i»KT LONG 

1 4PTATA 

1G1 

pic ccn Ra La sm icexe 

RD MUST LOG 

2 usm. 

lit 

pic con Ra La noshr noexe 

RO HRT LONS 

Total Spice Allocitrd 

680 




ENTRY POINTS 

Addrett Type Nsae References 

0-00000800 ffiSNINFILtmiN 

VARIABLES 


Address Type Nj*t Attributes References 


2-ODOCOOlC 

R*4 ATH 

14* 

If 


2-00000014 

R*4 AZ 

14* 

15 

IS 

2-00000000 

CTAR FLfrt 

3 

5* 

6A 

2-OOOOOOOC 

1*4 J 

8* 

9 


2-00000020 

R*4 FHIO 

11* 

IS 


2-00000018 

R*4 PHIZ 

14* 

16 


2-00000010 

R*4 HE 

14* 

IS 



A51 



CNSNHflUMAlN 


lakes 

Addrost Label 
1-00000091 S' 

0- ocoooooi ? 

1- 00000096 10' 
0-00009104 SO 


2-Feb-19S3 13:21:30 UW-11 FORTR/W V3.S-02 F«* 

IO-Ok- 1584 15:34:03 FSM:[STiVEN5.SHIiWF]3tSmtfIL.F()*t3 


Rrftrtncn 

9 10# 

13# 19 

IS 17# 

13 20# 


rwcriats AM) SUBROUTINES REFftENCEO- 


Type Net 


Wbmch 


FCttOXSE 20 

FOMOPOi 6 


i -4 

I LEY TO REFERENCE FLAGS I 

I • - Valor Modified I 

I # - Dtf inin? Reference I 

I A - Actual Ar^went, po«ibly wdified I 

I 0 ■ Data Initialization I 

I (n) - NucSir of accarrncn'nn lint I 

^ + 


OBffS OflllFIERS 

FCSTMN /CRO/LIS ffaiJNFlL.FCa 

/CHECWKBISiflS , IRJERFLC4 .KHfflERFLOfl 
/KRG> (N8SVK3aS,TRACEEfta) 

/Sr«3MD*(N0SfflT« t NaS0«;CE FORM) 

/sHcn«(M)?MPi(oassa,faiKaibE,ttv>) 

im />03_FLCAT1HS /14 /OPTIMIZE /VVWIU5S /WO LINES /CROSS REFERENCE /T4CmCH!tC_CCD£> /CO(TllLATIOiS»19 


C0MPILAT1CN STATISTICS 


Run Tier: 
Elwiod Tier: 
P*j« Fault! : 
Dynaic Ntasry: 


1.5$ seconds- 
2.33 seconds 
123 

123 pa^es 



wmtWHW tH HW t HtHl 2-FE8-1S33 13:12:33.4$ FSDO : [ STEVENS . STORAGE I G£NSM . DrtT ; 4 

t uummmmmwmm i 2-FE8-1583 13:12:50.45 FSM:fSTEVOfi.STQIW0£)GDfi(1.MT;4 


QENSM - D#=»T s ^ 


UrtHtHtitHtit H l W H tWK W 


2-FEB-13S3 1 3:12:33.43 FSOO : ( STaOfi . STCFAGEl GTNSH. MT ; 4 

2-FEB-1S35 13:12:50.43 rSOO : ( STEVENS . STCRACE 1 GEX'Sfl . DAT ; 4 

2-FES-1S23 1 3:12:50.45 FSOO:(STEVQfi.$TCF/VElGD131.MTi4 



l l ll l limilHHi ii tltmiHM 


mUrtti HtHUUH +titttrtlt 


A53 


r r 

/> 



.1 



DATA FOA 00333 CLASS. SCA STATE 4 
W£ AZ FK12 ATH PMO 


4 

0.3714 0.1234 0.0035 0.0004 -1.1447 

0.5435 0.5354 0.0035 0.0331 -1.1131 

0.7200 1.5111 0.0132 0.0033 -0.3325 

0.8314 1.5322 0.1242 0.0073 -0.7343 

1.0323 0.8724 0.5887- 0.0043 -0.2552 

1.2400 0.1511 1.4183 0.0018 0.27lf 




A54 






n 


ittiitti H rt tit i tiii i nnti t H k 


2*fTM?35 13:14:^.21 rS00:(STP^S.ST«WGC]G0O1.WTi5 
2-fIB-l?85 13:14:50.21 FSOO : [ STEENS . STOSAG£} G£KSM . OAT : 5 

2-4XB-1535 13:14:50.21 FSOOsISTSVEMS.STOfiMnGChSrt.Mr;S 



GENSM - O^T ; 5 


H 2-FEB-1S85 13:14:50.21 FS00:(STEVWS.STWUG£1GINSM.WT;5 

2-FE8-1393 13:14:50.21 FSOO:!ST£VO«.STOWGE1GEN5M.»T;5 

2-F£B-1535 13:14:50.21 FS00:(STW3iS.STC?AGE]GDO1.WT;5 






**** * H HI t m 


A55 




MTA res DD963 CLASS, SEA STATE 5 
ME Id PHIZ ATH PH10 

0^3543 0.5917 -0.0035 0.0319 -1.1439 

0.5134 2.7355 -0.0035 0.0086 -1.1467 
0.7029 2.4986 0.0053 0.01C8 -1.0103 

0.8743 2.1737 0.1031 0.1000 -0.7674 

1.0457 1.1326 0.5160 0.0052,-0.3152 
1.2229 0.2533 1.3551 0.0019 0.2378 





A56 



n 


ttWHW t WWIMWWWW 


2-FEB-1935 13:15:35.98 
2-FE8-1935 13:15:35.99 
2-FEB-1935 13:15:35.93 


fSDO:lSTEVENS.STORfi5EJCOCM.WT;J 
FSDOit ST EVENS. STQRAGE)G£3<SM.TAT;6 
FS00ilSTMKS.ST0BAKlGD<9t.#ftT;( 


>n«mtut»» »***Httm*****« 



GEN SM - DAT ; <s 


/ 


mniiiimmmw wtwM 


2-FEB-1933 13:15:35.99 
2 -FES-1933 13:13:35.93 
2-FEB-19S5 13:13:35.93 


FSOO:lSTEVENS.STOMGClGElS1.MTjS 

FEDO :( ST EVENS . ST WAGE 1 GSN31 . DTiT ; 6 
r 500 : ! ST EVENS. STCfACE 1 GDttl . OAT ; 6 


tftti*ii*t***MiH*Mtnk**** 

itirH*t********t*t*tH t*t*tt** 


r 


A57 



/ 

/ 


/ 


- '7r> 


DATA rCS M963 CLASS, SEA STATE 6 
HE AZ PHIZ ATH PHIO 


6 

0.3486 1.0005 -0.0035 0.0013 -1.1443 

0.4686 4.4116 -0.0035 0.0113 -1.1492 

0.6343 3.6532 0.0002 0.0143 -1.0743 

0.0143 2.0503 0.0546 0.0135 -0.8660 

0.3823 2.0422 0.2384 0.0034 -0.3174 

1.1257 0.5776 0.8234 0.0033 -0.0384 


i 

i 


■r' 


V 


i j 


; i 

j 

t • 






r 



2-FCB-15S3 13:40:23.50 rSOO:ISTEV£><S.SHIPSTUrF)VASLNX5.CCM;3 

2-FE8-1935 13:43:23.50 FSS0:tSTEVEW.SHlPSTlFn«fta>«5.CB1:3 

2-FEB-13S5 13:40:23.93 fSOO:lSTEVOS.5HIFSTUfr]'A«Utt5.Cai;3 





*» 




VASLNKS . COM ; 3 


ttHWHi ii i mi i imm ttw 
ttn+it ww mmwtwttHH 


2-FE3-19S3 13:40:23.50 FS»:[STaTO.SHIPSTtFF]UASU*3.C01;3 i wmmH wtHwwWtt: 

2 -FEB-1 933 13:40 : 23.93 FSM:[STaTNS.$HIPSUlFF)VftSLW5.CCM:3 mtWHitmttmwwwi 

2-FE8-19S3 13:40:23.50 FSCa:ISTWEWS.SKiFSIUFFlWafM.CtMi3 




7l59 





* SET VERIFY 

* LINK VA3CON.VPA1C,IN4.GENSM.TRANO.VASSTAT,MPLT2,PLOT3,SYS*UBRARY: INTLIB/LIB*D 

* SET MOVER I FY 





/ 

/ 

/ 

/ 



2-FEB-1S35 11:05:18.63 rSOOsISTEVEJfi.SHlPSTlFnWftlC.LlSil 

2-FE8-1333 11:05:13.68 FSD0:[STEV3IS.$MPSnifF)VFAlC.L!S:l 

2-FEB-19S5 11:05:18.68 FSDO:t$TEWNS.SmPSTUFFjtf’AlC.LIS;l 



i 




i ! 

j , 
i ; 



C . L. I S ; X. 



2-ftB-I335 11:05:18.63 FSOO :{ STEVENS.SHtF'STUFFJVPftlC.LIS;! HWWWW ihHHit+ntttiit 

W HttWHwm HtiwttitiiH 2-FE3-1385 11:05:18.68 FSOO : l STEV€HS.SHI PSTUFFIVPA1C.H S; 1 

lUDiii i mHHt i imuium 2-FEB-1385 11:05:18.68 FSOO: ( STEVENS. SHlFSTUFFl'vTAlC.LI Sjl 






/*' i • . ' •/• 

' ' i 


21 -Jm-1935 12:03:35 UVM1 FORTRAN V3.5-62 PJ?» 1 

IJ-Jm-1335 21:06:43 FSDO:[STEVD 4 S.S 14 .'PSnrr]VPAlC.FOR;l? 


0001 

C12345678S012345S7890...HC 

0002 

C 


0003 


amount uaswi iseas.k, jran.l, tic, w, so, rf\he.az, phiz, 

0004 


1 ATH.FHIO.JW 

0005 



0006 


0IHENSIO4 aS(30),aV(50).SlA(59).Zl(50).2r)U50).ZPOU50) 

000? 


DIK04SICN HE(10),AZ(10).PHIZ(10),ATH(10),PHIC(10), 91(500) 

0008 


OIHENSKN AC(500),T11S(500) 

0009 


DIMENSION acrr(50).TWTR(33),IS(30), 04143(50) 

0010 


DIMENSICN Wt4<500),V$0(500>,TW4t500>,TSt>(500),Tim(WO) 

0011 



0012 


REiV. W.!T1.!T2.1T3.IT4.IT5.IT6,IT7 1 ITB,1T9,IT10,IT11 

0013 


PEAL IT12,IT13,LIM,HTGUT,MIVRO 

0014 


REAL KPA,KP3,R?C 

0015 


CMRACTER*8,8UF 

0016 


CHARACTER*3,BUFF 

0017 



0018 


PR1NT*,'RIH4INS VPA18 AS MAIN PRPGRAF1' 

0019 


PRINT*, 'SEASTATE*',IS£AS 

0020 


FRINT*.'RP*',RP 

0021 


KPA=RP 

0022 


KPB=RP*1.3 

0023 


RPC=RP*2.0 

0024 



0025 


GO =61 DOT 

0026 


ZH=8.01 

002? 


IFIA6M 

0028 


RAY -6 

0029 


DO 1 1NT*1,30 

0030 


IS(1NT>=3 

0031 


DLK3(INT)*0 

0032 

1 

CU4TINDE 

0033 


OO 175 1D0*1,K 

0034 



0035 



0036 

C* 

TYPE*. 'DO YOU HOT A PLOT? NON)' 

0037 

C* 

ACCD’T«.FIA04 

0033 

C* 

TTPE*. 'INPUT INITIAL TIME' 

0039 

C* 

ACOEPTM 

0040 

c* 

TYPE*,.' INPUT T/H, ELC, ZW.TLAGh' 

0041 

c» 

ACCEPT*,n4,aC,ZH.L 

0042 



0043 


t«3600«PJY4(JIWf4) 

0044 



0045 


tll't 

0046 


TOO .01 

0047 


dHOO.Q 

0048 


dHOD=0.0 

0049 


OHTR=0.0 

0050 


(ffflTIWI.O 

0051 

04 

Ht«T 

0052 

04 

Wi=0.25 

0053 


TFMIIF360 

0054 

C* 

CRPHT’t 

0055 

C* 

0RPIf)»0.2 

0056 


A8RTHO6.0 '4.5 

003? 


CTNF4»-1.0 


A62 



V r 



W 


21-JW-1985 12:05:35 
19-JU-1585 21:06:« 


'.Wt-ll fORTWW V3.5-62 






Pj^t 


2 


0058 

0059 
0000 
00(1 
00(2 
00(3 

00(4 CW 
00(5 
00(6 
00(7 
00(3 
00(9 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 
0073 

0079 

0080 
0081 
0082 

0083 

0084 

0085 
COES 

0087 

0088 
0089 
0030 
0051 
0092 
0053 

0094 

0095 
C09S 
0097 
0093 

0099 

0100 
0101 
0102 

0103 

0104 

0105 
OIOS 

0107 

0108 
0109 
OlU 
0111 


0112 

0113 

0114 


Y1*0.0 

Y2=0.0 

Y3=0.0 

Y0*0.0 

YC0=0 .0 

T2*0.0l 

KP*0.01 

Tl*0.75 
TG*O.01 !0.01 

f lW M.O ! l*TURBaEMX IS ON 

IFIAG3=0 

IFIAG7*1 

lflfiG8-8 

1FIA69*0 

IFLAG1Q=0 

|flA611*0 

1FLA612=« 

IFIAG13=1 

TE=0.0 

FPIS1=0.0 

2W0T0T=32.2 

$1=40 

S2=0.0 

S3=0.0 


T1=T2*2 

T4=T2«2 

IT1=TC 

lT2=TC«C/2 

1T3=IC«3/S 

lT4-"TC**4/24 

1T5=ICHV120 

IT(=TC«(/720 

lT7=TC**7/5040 

1T8=TC** 3/40320. 8 

dHTR-0.0 

<w-o.o 

ZMO.O 

ZD=0.0 

700=0.0 

IT-l 

n€wn 

Toumw=o.o 
TOUTNUMO.O 
MMD=0.0 
J=79856423 
DO 5 1=1.50 
aS(l)=«0.0 

avd)=o.o 

SW(l)=0.O 
ZL( 1 1*40.0 
ZDUI)=0.0 
7COt(I)=0.0 
aoTd)=o.o 
TWtTR< 11=0.0 


A63 



wsw* 


0113 

one 

0117 

0118 
0113 
0179 
0171 

0177 

0173 

0174 
0173 

0178 
0177 
0173 

0179 

0130 

0131 
0137 

0133 

0134 
0133 
0130 
0137 
0133 
0133 
0149 
0141 
0147 

0143 

0144 
0143 
0140 
0147 
0149 

0149 
0130 
0151 
0137 

0153 

0154 

0155 

0150 
0137 
0159 
0133 
0100 
0101 
0107 

0103 

0104 

0105 
0100 
0107 
0109 
0103 

0170 

0171 


/ 



' ■ I 


■l '• I 


7 


31-Jin-198S 12:09:35 VAX-11 FMTIWJ U2.5-07 Pjo* 

19-JM-1395 71:00:49 FSDO;(ST£VEJ4S.SHIPSTUEF)VPA1C.FO(<:17 


3 CtNTIMJE 

DO 0 1=1,300 
AC(l)=0 
T1ME(I)*0 
0 COfflNUE 

04 401110(5,10) 

0410 Few (SX.'t'.ex.'si'iOx/Ht'.ex.'dH'.ax/TG'.ex.'TL', 

04 1 8X,'FJ>') 

04 HUTE(0,15) 

0415 FOfWtf (1X,10(' '),//) 


CALL INPUT1(&U,SA7,SA3,MII,IFLA53,SX) 

CALL INPW2<SSl,SS2,SS3.t,tn,EN,)«,AZ.PMZ,ATH,PHI0) 
CALL TUS9ULD«(TGE£,IC,VTZ,PR£aC,K4( , ,FNO) 

60 TO 149 

79 CALL INPOTl(SAl,SA7,SA3,t,tlI,)FlAC3,SX) 

CALL lHrUT?(SSl,SS7,SS3,t, til, EN.WE.AZ, FIIIZ, ATH.PHI0) 
CALL TUR£aOC£(TG£E,TC,VTZ,PSCCrC,RfP,»;0) 

C tmtttti* FllEHTPATH C01WO LOS1C ********** 

30 SS1Z=SS1-Z 

SS2ZD=SS2-ZD 

IF (IFLAG10 .CO. 1 .A440. TC4 .LT. 1.5) CO tO 59 

IF (SS2 .GT. 5.5) GO TO 35 

IF (SS1Z .GT. -9.0 .AK>. SS2 .GT. 4.5) GO TO 49 

IF (SS1Z .LT. -3.0) 60 TO 45 

IF (Z .LT. ASSTH) GO TO 43 - 

IF (SS1Z .GT. CTJ34) CO TO 50 " 1 ' 

IF (1FLAG3 .CD. 0) THD4 
TC1=0.0 
!FIAG3=1 
1FIAG7=9 
1FIAG9--0 
1CLAG10=0 
1FLAG11=0 
1FLAG13=0 
IFIAG13=9 
35P=S1 
OT4H=-1.0 
EM) !F 

S1=SS14((5CP-SS1 )4(S6P-SS1 )*C0S(6.2332*TE1) )/2.0 
S2=S37 
S3=SS3 
KP=XPB 

IF (CXP8 .EQ. 1) KP=XFA 
TE1=TE14TC 

IF (TCI .GT. 0.5) TE1--0.5 

FPLS=3 

GO TO 65 

35 IE (IflAGll .EQ. 0) THEN 


7164 



/ 


/ 






/ 


/ 

/ . 


. 1 .' 

/•■/ 


* 17 

!/. 




-7^ 

f 



i 

i 




/ 

7 ' 




7 

/ 

/■ 


.-/ 

/: 


WSM> 


21-Jao-l?S3 12:09:35 UAX-11 F08TVH 03.5-62 Pa 9 » 4 

19— Jan-1533 21:Cb:43 FSD0:[STD.DIS.SHIPSTLFF)OPAlC.F0I!;17 


0172 

0173 

0174 

0175 
0170 

0177 

0178 
0173 
0100 
0181 
0132 

0183 

0184 
0183 

oiee 

0187 

0183 

0183 

0190 

0191 
0132 

0193 

0194 

0195 

0196 

0197 

0198 

0199 
0290 
0201 
0202 

0203 

0204 

0205 

0206 
0207 
0203 

0209 

0210 
0211 
0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 
0221 
0222 

0223 

0224 
0223 
0226 
0227 
0223 


TE3*0.0 
iriASM 
IFW67*0 
1FLA5S*0 
IFIAC10=0 
IFIAG11*1 
IFLAG12=0 
IFLAG13*0 
SGAPP=S1 
GWP«6.0 
CTMF«-1.0 
DO IF 

Sl*S&W((^-«PP*CCS(6.2832tTE3))/2.0 

S2»SS2 

S3*SS3 

KP*F?C 

TF3‘TI3+TC 

IF (TE3 .CT. 0.5) TE3*9.5 

FPli=l 

60 TO 65 

40 IF (IFIAG12 .Efl. 0) TIE31 
TE=0.0 
m*G3*o 
IFLfiB7=0 
IFU£3*0 
IFLA610=O 

IF (IFEAG11 .EQ. 1) TKN 
SGPP*S1 
GPP-3.0 

as£ IF (IFLAG11 .EQ. 0) TIEN 
SGPP'Sl 
6PP=3.0 

ENO IF 
IFLAGll^O 
IFLAG1M 
IFIAGIM 

nw»-i.o 

DO IF 

Sl*S6PM(CPf-GPP*CCS(6.2832*TE))/2.0 

S2=SS2 

S>SS3 

FP=KP8 

TE-TE+TC 

IF (TE .CT. 0.5) TE=0.5 

fpes*2 

CO TO 65 

45 IF (IFEAG13 -E0. 0) TIEN 
TE2=0.0 
SSWP.S1 
SS<WP«SA1-S$«P 
1FLAG3=0* 

IFl/W9*0 

IFlAGlOiO 

IFIAG11=0 

IFLAGl^O 
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i 


1 


r 

v_ 


(i 


0229 

0238 

0231 

0232 

0233 

0234 

0235 
0234 
023? 

0233 

0239 

0240 

0241 

0242 

0243 

0244 

0245 
0244 
0247 
0243 

0249 

0250 

0251 

0252 

0253 

0254 

0255 
0254 
0257 
0253 
0253 

0240 

0241 

0242 

0243 

0244 

0245 
0244 
0247 
0243 
0249 

0270 

0271 

0272 

0273 

0274 

0275 
0274 

0277 

0278 

0279 

0230 
0281 
0232 
0283 

0234 
0285 


IFIAG13=1 
CINK=-1,0 
BO IF 

IF (IFLA67 .EQ. 1) GO TO 46 
Sl*SSaPF(SS(flP-sam£OS(4.2832*TE2) 1/2.0 
GOTO 47 

46 S1-SA1 

47 S2=SA2 
S3=M 

IF (Z40.5 .IT. SI) TICK 
KP=KP3 

ELSE IF (240. 5 .G£. SI) TUN 
KP*KPA 
0® IF 
TE2=TE24TC 

IF (TE2 .GT. 0.5) THEN 
TE2=0.5 
IFLAG7=1 
ENO IF 
FPL$*0 
GO TO. a 

50 IF (IFV.S10 .EQ. 0) THEN 
k?=wa 
FPIS= 5 
]FUC2=0 
IFLAS7*1 
lFlAGM 
IFLAulO-1 
1FUG11=0 
IFIAS12=0 
IfLAG13=0 
C1H4-1.7 
TE4=0.0 
GAP=S1 
END IF 
HLIN=-1.0 

IF (TE4 .LT. 0.25) THIN 

S1=GAP449.WCO$(12.5664*TE4)-9.0)/2.0 

FPLS=4 

as£ IF (TE4 .GE. 0.25 .»3. TE4 .LT. 1.25) THEN 
Sl=GAP-3.0 
FF1S-? 

ELSE IF (TE4 .GE. 1.25) THEN 

Sl=£AP4(9.C*COS(12.5464*(TE4-l))-9. 01/2.0 
FPLS=3 

END IF 
TE4=TI44TC 

45 IF ((SSI-2) .GE. 0.0) GO TO 150 

DO 70 1=1, 2,-1 

SLS(!)=SIS(I-1) 

aV(I)=SLV(I-l) 

SIA(I)=SU(I-1) 

70 caninijE 
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0286 

SLS(1K1 


0287 

SLV(1)*S2 

1 

0288 

SLA(1)=S3 


0289 



0250 C 

HtUMMi PILOT T.F.F. 


0291 



0232 

A6 J dH 


0293 

A1*Y0 


0299 

B1«1/TG 


0295 

P2=0.0 


0236 



0297 

C1--A1 


0238 

C23*-atCl 


0239 

C3S»-E1*C28-B2*C2 


0330 

C9S»-B1*C38-B2*C39 


0331 

CSB»-B1*C9B-B2*C98 


0302 

C68«-B1»C58-J2*C5B 


0303 

C78=-B1*C68-82*C68 


0309 

C3S=-B1*C7B-B2*C7S 

j 

0305 


1 

0306 

C2A«1.0 


0397 

C3A*-B1*C2A' 


0308 

C9A*-B1*C3A 

j 

0303 

C5A»-B1*C9A-52*C4A 

f 

0310 

C6A*-Bl*OA-B2*C5A 


0311 

C?A»-B1*C6A-B2*C6A 


0312 

C8A*-B1*C7A-B2*C7A 

/ 

0313 


( 

0319 

YlA=Cl+C2S*ITl+C53*IT2+C9B*IT3tC5B*lT94C€B*IT5 


0315 

Y1B=C7B* 1 T 6+CSS* I T7 


or, 6 

Y2A*C23FC389|THC9B*IT29C5B*IT39C6B*1T9 

• 

0317 

Y2S=C7B»IT5K8MIT6 

1 

0318 

Y39*C33+C93*mtCSS*IT24C6BW3 

£ 

0313 

Y33=C7B* IT4+CSB* 1 T5 


0320 C 

Y9A=C9BtC58»imC63*!T29C7B*lT39C8B*IT9 


03 a 


i 

0322 

YAlA<2A*lT14C3A*lT2F04A*IT3tCSA*IT9+C6A*IT5FC7A*lT6 


0323 

YA1B<8A*!T7 


0329 

YAWA1A4YA1B 


0325 

YAaC2A+C3A*ITHC9A*lT2tCM*lT3FC&9*IT9(C7A*n54t8A*IT6 


0326 

YA3=C3A9C9A9ITlFC5A>IT2FCfA9IT3iC7A*IT9+03A*IT5 


0327 C 

YA9*C9AFC53*imC6A*!T24C7A*lT3K8A*lT9 


0329 



0323 

C2A=A0 

- 

0330 

C3A*-B1*C2A 

! 

03 a 

C9A-B19C39 


0332 

C5A*-fil*C4A-E2*C9A 


0333 

C6A*-61*C5A-B2*C5A 


0339 

C7A>-B1*C6A-B2*C£A 


0335 

C8A*-B1*C7A-B2»C7A 


0336 



0337 

YAMA<2AHTl»C3A*IT2+C9A*IT3FC»<lT99C6A*lT5tC7A*IT6 


0338 

YM19=CaWIT7 


0339 

YAA1*YFA1AFYAA18 


0390 

YM2<2Att3A*m«9AW2K5AWaC«MT99C7A*ITMCOft*IT6 


0391 

YAA30A4C9Aini9C5A*IT2FCCA*IT3tC?A*iT9IC8A*IT5 


0392 C 

YM9<9AFC5WimC6A9ITaC7A*IT3FC9A*IT9 
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034] 

0344 

0343 

0346 

0347 C 

0348 

0349 

0330 

0331 

0332 

0333 C 

0334 

033$ C 

0356 C 

0357 C 

0338 

0339 C 

0360 

0361 C 

0362 

0363 

0364 

0365 

0366 

0367 
0363 

0369 

0370 

0371 C 

0372 C 

0373 

0374 74 
0373 

0376 

0377 73 
0373 
0379 

0303 

0381 C 

0332 

0383 80 

0384 

0335 

0336 
0387 
0383 

0389 

0390 
0331 
0392 

0333 
0354 
0353 
0396 

0337 

0333 C 
0333 C 


Y1*Y1A4Y18 

YM2MY23 

YJ>Y3A+Y3B 

Y4«Y4AfYA4 

01*(1/Tt-1/T2) 

DM/(Tl*I2) 

03?(KP*Tl/TG) 

Ht8*YA2+(l/TU*YAl 

HtB»(YA3F01*YA2-02*YAl) 

HtA 3 (Y2faS(L)*YA2)t(l/TL)*(YlfaS(L)4YAl) 

HtA=(Y34a$(U*YA3H01*(Y2+aS(l)JtYA2)-D2*(Ylta$(l)*YAl) 

HtP03*-D2*(Y14YftU) 

HlP0A*(Y3+YAA3)401*(Y24YAA2) 

mPO^(im)*(YH(aS(L)-ZUU)*YAl) 

mpoe*-o24 (Y»( asi u -zuu j*yai ) 

HtP«=r2t($lS(l)-ZUl)>*YA2 

HtPOA=(Y34(aS(l)-ZUU)*rA3)tO!*(Y24(aS(U-2Ul))*Y«) 
HlP0=03*6 Ht POA+HtFOS) 

Y0*YHYAA14fHP 

YD0*Y24YAA2 


LW*1.0 

IF (HtPO .60. 0.0) HtPO’O.OOOCl 
IF (1FLAG3 .63. 1) GO TC 74 
HUH=-((TO-1)*1.5) 

IF (aiM ,GT. -0.05) aiM=-0.05 

aim-i.o 

IF (HtPO ,GT. HW) GO TO 73 

LW=aifVH(?a 

GO TO 80 

IF (HtPO .LT. (W-l)) GO TO 80 
lW=(W-l)/HtPO 


tlumiutl A/C T.F.F. Hit ti kHU 

AA0=!(2W)TDTf((TGF£+[VW)/LW))/aC)*UW03*HtA 

AAl=Z*ZH/aCtZW(l/aC't2W)CTD 

AA2=Z*(l/ELC43J)tZ0 

M3=Z 

tm=((2M)T0T)((TaC4im , )/lW))AlC)*LW»03»HtB 

BA2«ZH/ELC 

&0=ZW1/8LC 

AA0P=U»03»HtA 

0AlP=LIH*D3*HtB 

poor^ptayptf-zLai) 

T0UT=P0lfTfI 

mo'Ji * ( ( » 1 1 wcumtrji ) /r.>ow 

IF (TOOT .GT. T0U7M0X) TOUTWX^TOUT 
IF (TK/T .LT. TCOmiN) T0UTXIN*TtWT 
IF (TOUT .LE. 1.21 MIS. TOOT .£T. 1.101) T|£N 
CO TO 83 
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0400 

0401 

0402 

0403 

0404 
0403 

0406 

0407 

0408 C83 

0409 C 

0410 C 

0411 C 

0412 C 

0413 C 

0414 C 

0415 C 

0416 C 

0417 C 

0418 C 

0419 E 

0420 C 
04a C 

0422 C 

0423 C 

0424 C 
0423 C 

0426 C 

0427 C 

0428 C 

0429 C 

0430 90 

0431 

0432 

0433 

0434 

0435 

0436 

0437 
0433 

0439 

0440 

0441 

0442 

0443 

0444 
0443 

0446 

0447 

0448 

0449 

0450 

0451 

0452 93 

0453 

0454 
0453 
0436 


IF (TOUT .U. 1.101 .ft®. TOUT .GT. 1.001) THEN 
GO TO 30 

asc IF (TOUT .LE. 1.001 .(H). TOUT .GT. 0.901) TKN 
GO TO 93 

ELSE IF (TOUT .LE. 0.901 .(H). TOUT .GT. 0.800) THEN 
GO TO 100 
DO IF 
GO TO 105 

IF (TO UT .LE. 1.201 .ft®. TOUT .GT. 1.191) TKN 
IS(l)*IS(l)fl 

ELSE IF (TOUT .IE. 1.191 .(H). TOUT .GT. 1.181) TUN 
!S(2)»!S(2M 

ELSE IF (TOUT .LE. 1.181 .ft®. TOUT .GT. 1.171) TKN 
IS(3WS(3)H 

ELSE IF (TOUT ,l£. 1.171 .(H). TOUT .GT. 1.161) THEN 
IS(4)=IS(4)41 

RSE IF (TOUT .LE. 1.161 .(H). TOUT .GT. 1.151) THEN 
1S(5)»!S(5)H 

ELSE IF (TOUT .LE. 1.131 ,(H>. TOUT .GT. 1.141) TKN 
IS(6)-IS(6)41 

ELSE IF (TOUT .LE. 1.141 .ft'®. TOUT .GT. 1.13) TKN 
IS(7)*!S(7M. 

ELSE IF (TOUT .LE. 1.131 .(H). TOUT .GT. 1.121) TKN 
IS(B)=IS(8)« 

ELSE IF (TOUT .LE. 1.121 .(H). TOUT .GT. 1.111) TKN 
ISC9)=IS(9)+I 

ELSE IF (TOUT .LE. 1.111 .ft®. TOUT .GT. 1.101) THEN 
I S( 1 0 ) =1 S( 1 0)4-1 
D® IF 
GO TO 105 

IF (TOUT .IE. 1.101 .(tO. TOUT .GT. 1.091) THIN 
IS(ll)=IS(ll)tl 

ELSE IF (TCUT .LE. 1.091 MO. TOUT .GT. 1.081). THEN 
ISU2)=IS(12)U 

ELSE IF (TOUT .LE. 1.081 .ft®. TCUT .GT. 1.071! THEN 
IS(13)=1S<13)+1 

as£ IF (TOUT .IE. 1.071 .ft®. T'JUT .GT. 1.061) THDI 
IS(14)=IS(14)H 

ELSE IF (TOUT .LE. 1.061 .ft®. TOUT .GT. 1.051) TO! 
IS(15)-IS(15)41 . 

ttSE IF (TOUT- .LE. 1.051 .ft®. TOUT .GT. 1.041) THEN 
IS(16)«IS(16M 

aSE IF (TOUT .LE. 1.041 .ft®. TOUT .GT. 1.031) TKN 
IS(17)«IS(17)U 

aSE IF (TOUT ,LE. 1.031 .ft®. TOUT .GT. 1.021) THEN 
IS(18)=IS(13)fl 

ELSE IF (TCUT -LE. 1.021 .ft®. TOUT .GT. 1.011) THEN 
IS(19)=IS(19)F1 

ELSE IF (TOUT .LE. 1.011 .M. TOUT .GT. 1.001) THEN 
1S(20)-IS(20)H 
DO IF 
GO TO 105 

IF (TOUT .LE. 1.001 .ft®. TOUT .GT. 0.991) THEN 

isiajMstam 

ELSE IF (TCUT .LE. 0.931 .ft®, TOUT .GT. 0.931) THDI 
1S(22)’1S(22)F1 

ELSE IF (TCUT .LE. 0.931 .ft®. TOUT .GT. 0.371) TKN 
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0457 


IS(23)*IS(23)F1 

0458 


QSE If (TOUT .LE. 0.971 .HO. TOUT .GT. 0.961) THEN 

0459 


IS(24)-IS(24)tl 

0460 


ELSE IF (TOUT .16. 0.961 .AND. TOUT .GT. 0.951) THEN. 

0461 


IS(25)*IS(Z5)U 

0462 


ELSE IF (TOUT .IE. 0.951 .AND. TOUT .GT. 0.941) THEN 

0463 


!S(2S)-!S(2G)fl 

0464 


asE if (Tan .le. 0.941 .ho. tout .gt. 0.931 ) then 

0465 


IS(27)=IS(27)tl 

0466 


ELSE IF (TOUT .LE. 0.931 .AND. TOUT .GT. 0.921) TUN 

0467 


IS(23NS(28)H 

0463 


aSE IE (TOUT .LE. 0.921 .AND. TOUT .GT. 0.911) THEN 

0469 


IS(23)=IS(29M 

0470 


ELSE IF (TOUT .LE. 0.911 .HO. TOUT .GT. 0.901) THEN 

0471 


IS(38)=IS(30)I1 

0472 


END IF 

0473 


GO TO 105 

0474 

100 

IF (TOUT .LE. 0.901 .HO. TOUT .GT. 0.891) THIN 

0475 


IS(31)=IS(31)+1 

0476 


ELSE IF (TOUT .LE. 0.891 .AND. TOUT .GT. 0.881) THEN 

0477 


IS(32)=IS(32)tl 

0473 


ELSE IF (TOUT .LE. 0.881 .AND. TOUT .GT. 0.871) THEN 

0479 


1S(33)=IS(33)F1 

0430 


ELSE IF (TOUT .LE. 0.871 .AND. TOUT .GT. 0.861) THEN 

0481 


!S(34)=IS(34)fl 

0482 


aS£. IF (TOUT .IE. 0.861 .HO. TOUT .GT. 0.851) THEN 

0433 


iS(35)’iS(35)tl 

0434 


ELSE IF (TOUT .LE. 0.851 .AND.. TOUT .GT. 0.841) THEN 

0435 


IS(3G)=IS(36)tl 

0486 


ELSE IF (TOUT .LE. 0.841 .AND. TOUT .GT. 0.831) THEN 

0487 


IS(37)=IS(37)fl 

0438 


ELSE IF (TOUT .LE. 0.831 .AND. TOUT .GT. 0.821) THEN 

0489 


IS(3S)=lS(33)tl 

0490 


ELSE. IF (TOUT .LE. 0.e2I .HO. TOUT .GT. 0.811) THEN 

0491 


IS(39)=IS(39)+1 

0492 


ELSE IF (TOUT .LE. 0.811 .AND. TOUT .GT. 0.801) MN 

0433 


IS(40)=!S(40m 

0494 


END IF 

0495 



0496 

105 

C1A=AA3 

0497 


C2A=AA2-£&3*C1A 

0498 


C3A«AA1 -BA2*C1A-BA3*C2A 

0499 


C4A=AA3-BAl*ClA-fi‘.2*C2A-B43*C39 

0500 


C59»-BA1*C»-BA2*C3A-EB3*C4A 

0501 


C6A ? *&41*C3A-fiA2*C4A-8A3*C5A 

0502 


C7A*-E«<C4A-EW2<C5A-GA3*C6A 

0503 


C3A*-&A1*C5T-GA2*C6A-GA3*C7A 

0504 



0505 


2A1 =C1AH C2A* F T1 +CiA* I T2+C Vrt I T34 C5A* 1 T44C£A* ITS 

0506 


ZA2*C7A4|T6+CaA*IT7 

0507 


ZBl=82MC3ArtimC4A*|T2+C5A*IT3fC«*IT4 

0509 


ZB2=C7A*IT5+C£M!T6 

0509 


2C1 =C3A+C4M ITHCM* IT2FC6A4 1 T3 

0510 


2C2-C7A4IT4+C8ft*IT5 

0511 



0512 


Z»ZAI4ZA2 

0513 


ZD’ZBIFZBZ 
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0514 


ZtD*ZClfZC2 

0515 

0516 


DO 110 K.2,-1 

0517 


zujkui-i> 

0518 


zoun*zot(M) 

0519 


ZOOUlWOOUI-l) 

0520 

no 

C34TIISJE- 

osa 


ZL(1)*Z 

0522 


ZDUll'ZD 

0523 


ZOOUD-ZOD 

0524 


DaMP=(Zr*0.45*GHD0T) 

0525 

0526 


<HrMa.s(u-zuu 

0527 


ERRCR=Rtf4(J> 

0523 


ERJt*0.5-E?J40R 

0529 


(W’dHTBlRNP 

0530 

c 

d5+=dHTfi+0 . 1*ERH 

0531 

0532 


dK>TR*(aV(L)-ZDUL)) 

0523 


(K>=tMOTR 

0534 

c 

cK5=dHOTR40.12*UX 

0535 

0536 

c 

(H»=<S3-HtDD) 

0537 


10SN 0W4GIHG TO THE C04TitlOUS PR1MTOVT MODE 

0533 

C04 

0539 

COT 

SUBSTITUTE 6U4KS FOR THE Ol'S CWNGE 69 

0540 

C04 

04 L11C 59 TO' 71, THE 65 04 LIKE 256 TO 71 

0541 

04120 

IF <t .IT, HI) 00 TO 71 

0542 

04123 

TYPE 130. t.Sl.Z, SSI, TGEE, POUT, H<PO 

0543 

04130 

FORMAT (1X,?{F10.5,X)) 

0544 

04 

HT44Tttft 

0545 



0546 

C*133 

1 IF (t .LT. GRPHT) GO TO 71 

0547 

C* 

S4(1T1=SS1 

0543 

c* 

ACCT)=Z 

0549 

C* 

TIH£EUT) s t 

0550 

C* 

IMT41 

0531 

c* 

GRPKT=GRPH7EGRPW) 

0532 

140 

t*t4TC 

0553 


KHEAN*KKBt4tl 

ay 


IF {(SS2-2D) .GT. M1VRD1 MIVRD=( SS2-ZD) 

0535 


IF «t-TC> ,CE. TF) GO TO 150 

0556 

145 

GO TO 20 

0557 

150 

PRINT*, 'T/14 >',W,' ac =',ac,' 34 ='.34 

2353 


PRINT*, 'T.D. VU -',(SS2-Z0),'TIKE =',<t-tll) 

0539 


PRINT*, 'HTOUT=' .HTOUT,' FPIS*' ,FPLS,' RUN NO.:', 100 

0360 


PRINT*, 'ISia)*',IS(21),' FPIS1*' ,FflSl. 

0561 

C* 

IF (FL6G4. .EQ. 6.0) GO TO 1E0 

0362 

C* 

TFI1MII438 

0563 

CC* 

TFIN=IFIX(t+l> 

0364 


WRITE (2.153) TII ,T-TI I ,(SS2-ZD) ,t11VRD,MT0UT ,TOUTMAX, 

0565 


1 TCUIHIM 

0566 

155 

FORMAT (lX,2(3,45.:i.2(2X,F5.1),3(2X,F6.3)) 

0367 


IF (1FIA65 .EO. 0) GO TO 175 

0563 


aW-TI! 

0569 


VTD=(S32-ZD) 

0570 


A71 


/ 


n 


• < 


C 


-I 



mw 






13-JaiH 


0571 


R4(=IC0 


0572 


m« Kmi=t (oth itriTi/wo 


0573 


vw(ntD=((iitH)*i*t«(m-i)«T0)/itf»4 


0574 




0575 


IF (OV .60’. 11 GO 70 165 


0578 


IF (IW .GT. 2) GO '0 160 


0577 


TSC(Bt(W(FlTl-Tm' W))**24(aT-Ht4(PW»H21**0.5 


0578 


VS0(I34()*< (Vr01-tft«tMl) )**2i(VIW4W(»N) )«2)**0.3 


0573 


GO TO 165 


0580 

0581 

160 

TS0ft*((TS0tl0fH)H2)*(WH*2)4(FlT*TlW( 67401*42) 


0532 


TSS(IW)=<TSOV(»4H))**0.5 


0583 


V$CA«((VS0< 044-1 )**21* ( ttt*-2)+tVTt>-M-M( R44) 1 **2) 


0584 


VS0(Ht4) s (VS0A/(BN-l>)**8.5 


0585 




0586 

165 

Tiro(f344)=»44 


0587 


aii*ar 


0583 



7 

0583 


W1I76 ( 3,170) K744,D44( B44) ,TSP( B44) ,U44(KI44) ,VSD( 6144) 


0530 

170 

FORMAT (1X,:?,4(2X,F8.4)) 

* 

0531 

175 

CCNTINUE 


0532 


IF ([FLAGS .60. 0) GO TO 180 

- 

0533 


6£G=0 


0534 


TFD4=a»T(t) 

I 

0535 


IT«R 


0538 


CALL a0T3(TI.‘tt,TH4,T$0,1444 l VSD l E66.TFIN,IT t JAY,KA'f) 


0537 

188 

rxoo=i 


0538 


DO 185 1=1135,735,-10 

.V 

0533 


MTOH»0)=amT(l)/1000 


C5S0 


KCOO=HX)OH 


0801 

185 

CCt(TI!4UE 


0802 


ITOT=0.0 


0803 


DO ISO 1=1,40 


0804 


1T0T*1T0UIS(I) 


0805 

130 

COITINL? 


0808 


DO 135 1=1,40 


0807 


SICT(I) =a»T( IS( 1 1 1/aOATUTOT) 


0609 

135 

C0KTI7VJE 


0693 


TF!N=1.2 


0610 


TI1'0.8 


0611 


IT=40 


0612 


CALL aOT2(MTR, SUIT, SL0T.TI 1 ,TF1N, IT , JAY.KAY.TVl. 31,1 


0813 


7 FPfl.lSMS) 


0614 

C* 

CALL fn.0Tl(TI7CE,S1,AC,TIl .TF1N.IT) 


0615 


GALL TIHCtBUF) 


0618 


GALL DAT6(BUFF) 


0817 




0618 


MUTE (2,200) KJF.EUFF 


0613 

200 

rCRTWT (1X,2(A121) 


0620 


«n*n 


0621 


00 
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PROGRAM SECTItHS 


0 tcooc 

1 4P0ATA 
ZtlOCW. 


BvtM Auntotn 


c-*q pjr «w ^CL $Wt EXE RD N32.T LOJG 
Z "cSSia S*«*X£ rs«wtlo« 

15116 pic cm na ia noshr wsxt » L(W 


Tgul Sp*c* Allocitui 


ENTRT MINTS 

Addrns Typ* Na» 
C-0002CQ03 W31? 


Rtfrrtncn 


VARIABLES 

Addm» Typ» Hjm 

2-000047A8 R‘4 . AO 
2-000047AC R*4 PI 
2-00004853 R*4 MO 
2-005048/8 R‘4 MOP 
2-00004860 R*4 Ml 

2-00034864 R‘4 M2 
2-30004868 R‘4 M3 _ 
2-0000476C R‘4 ABRTH 
2-00004693 R‘4 ABRTHT 
2-00004780 R‘4 81 


2-00004784 R*4 82 

2-OC00486C R‘4 f»l 
2-O000487C R*4 BUP 
2-00004370 R‘4 BA2 
2-00004874 R*4 BA3 

2-000048C8 R*4 KG 
2-00004600 CHOR Elf 
2-00004608 CHAR DUFF 
2-00004783 R‘< Cl 
2-00004833 R‘4 CIA 

2-000047C4 R*4 C2 
2-0000 4 70C R“ C2A 

2-C00047EC R“ C23 
2-000047EO R*4 C3A 


2-00004700 R*< C33 


Attribute R»f*r«icM 


232* 

325 

253* 

237 

383* 

433 

351* 

333 

334* 

433 

335* 

437 

386* 

436 

146 


56* 


254* 

238 

307 

303 

332 

333 

255* 

233 

310 

311 

387* 

453 

352* 

333 

383* 

438 

383* 

43? 


239 

309 

334 

300 

310 

335 

301 

311 

302 

312 

303 

330 

304 

331 

300 

312 

500 

301 

332 

501 

302 

333 

502 

303 

334 

503 

304 

335 

309 

499 

433 

500 

433 

501 

500 

502 

501 

503 

502 

503 

618 






618 






314 

458 

493 

505 




322 

499 

325 

509 

323* 

505 

330 

507 

337 

340 

314 

322 

498* 

316 

325 

433 

326 

500 

330* 

31 

331 

505 

337 

507 

314 

316 

318 





■v=y 



UA9? 




21-J»n-lS35 12:09:33 
19-J»-19S5 21:06:48 
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2-000047E4 

R*4 

C4A 

308* 

309(2) 

322 

325 

326 

331* 

332(2) 

337 




340 

341 

499* 

500 

501 

502 

505 

50? 




509 








2*00004708 

RM 

C48 

300* 

301(2) 

314 

316 

318 




2-00004718 

R*4 

CSA 

309* 

310(2) 

322 

325 

326 

332= 

333(2) 

337 




340 

341 

500* 

501 

502 

503 

505 

50? 




509 








2-000947CC 

RM 

CSB 

301* 

302(2). 

314 

316 

318 




2-00004 7EC 

RM 

C6A 

310* 

311(2) 

322 

325 

326 

333* 

334(2) 

337 




340 

341 

501* 

502 

503 

505 

507 

509 

2-00004700 

R*4 

C6B 

332* 

303(2) 

314 

316 

318 




2-C00047FO 

RM 

C7fl 

311* 

312(2) 

322 

325 

326 

334* 

335(2) 

33? 




340 

341 

502* 

503 

506 

508 

510 


2-000047M 

R*4 

C78 

303* 

304(2) 

315 

317 

313 




2-00004764 

RM 

CM 

312* 

323 

325 

326 

335* 

338 

340 

341 




503* 

506 

508 

510 





2-00004703 

RM 

C88 

304* 

315 

317 

319 





2-00004778 

RM 

CKPB 

166 








2-00004S5C 

RM 

CTNH 

57* 

147 

159* 

162* 

209* 

230* 

262= 


2-00004833 

RM 

01 

349* 








2-00004S3C 

RM 

02 

350* 








2-00004840 

RM 

03 

351* 

362 

383 

387 

391 

392 



2-0000485C 

R*4 

DAMP 

383 

387 

524* 






2-00004704 

RM 

OP 

56* 

292 

529* 






2-00004684 

RM 

DHD 

47* 

533* 







2-00004688 

R*4 

DHOO 

48* 








2-00004690 

RM 

CHOIR 

50* 

532= 

533 






2-0000468C 

RM 

DHTR 

49* 

95* 

526* 

529 





AP-00000014# RM 

ELC 

3 

383 

384(2) 

3S5 

337 

383 

389 

557 




612A 








2-00004748 

R*4 

m 

123A 

135A 







2-0000<8A8 

R«4 

ERR 

528* 








2-00004844 

R*4 

ERROR 

527= 

523 







2-00004680 

RM 

ELAG2 

67* 








2-000049AC 

RM 

ar 

568= 

572 

577 

581 

587 




2-00004623 

R*4 

ail 

577 

587* 







2-0000477C 

RM 

pas 

169= 

190* 

217* 

249* 

254* 

269= 

272* 

275’ 




559 








2-00004668 

RM 

FPLSl 

78= 

560 







2-C00047AO 

RM 

6V> 

264= 

268 

271 

274 





2-00004783 

RM 

GAPP 

101= 

184(2) 







AP-0000001C* RM 

GO 

3 

25= 







2-00004660 

RM 

GKiOT 

25 

524 







2-00004750 

RM 

GPP 

201= 

204= 

211(2) 






2-C00047M 

R*4 

aw 

266* 

373= 

374 

375 





2-00004848 

RM 

HTA 

354= 

383 

391 






2-00004844 

RM 

KI8 

352= 

387 

392 






2-00004854 

Rt4 

KIPO 

362= 

363(2)= 

374 

375 

377 

378 



2-00004850 

R*4 

KIPOA 

360= 

•362 








A74 



3553 


I 
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2-00004340 

R*4 

HTP08 

358* 

362 







2-00004728 

1*4 

I 

106* 

107 

103 

109 

110 

111 

112 

113 




114 

116* 

117 

118 

281* 

282(2) 

283(2) 

284(2) 




516* 

517(21 

518(2) 

519(2) 

598* 

599 

603* 

604 




606* 

607(21 







2-00004674 

1*4 

100 

33* 

559 

571 






2-00004604 

1*4 

1F16610 

72* 

142 

154* 

176* 

1 sa- 

226* 

252 

258* 

2-00004608 

1*4 

1FL6C11 

73* 

155* 

171 

177* 

ns 

202 

206* 

227* 




259* 








2-00004600 

1*4 

IFL6G12 

74* 

156* 

178* 

193 

207* 

223* 

269* 


2-00004660 

1*4 

IFLA513 

75* 

157* 

179* 

208* 

220 

229* 

261* 


2-00004604 

1*4 

IFL6G3 

68* 

149 

151* 

173* 

195* 

224* 

255* 


2-00004663 

1*4 

IF16G5 

27* 

567 

532 






2-00004608 

1*4 

1FL657 

69* 

152= 

174* 

196* 

233 

247* 

256* 


2-OOC646CC 

1*4 

I FLAGS 

70* 

1284 

134A 






2-00004620 

1*4 

IFLA63 

71* 

153* 

175- 

197* 

225* 

257* 

370 


2-00004670 

1*4 

IMF 

29* 

30 

31 






ftP-000300040 1*4 

ISEAS 

3 

19 

612A 






2-00004714 

1*4 

IT 

100* 

595* 

596A 

611* 

612A 




2-00004614 

8*4 

IT1 

12 

67* 

314 

316 

318 

322 

325 

326 




337 

340 

341 

505 

507 

509 



2-00004633 

R*4 

IT10 

12 








2-0P00463C 

R*4 

mi 

12 








2-00304640 

R*4 

m 2 

13 








2-00004644 

R*4 

1T13 

13 








2-00034518 

R*4 

IT2 

12 

88* 

314 

316. 

3)8 

322 

325 

326 




337 

349 

341 

505 

597 

509 



2-OOOOi61C 

R*4 

IT3 

12 

89* 

314 

316 

318 

322 

325 

326 




337 

340 

341 

505 

507 

509 



2-00004620 

R*4 

1T4 

12 

90* 

314 

318 

319 

322 

325 

326 




337 

340 

341 

505 

507 

510 



2-00004624 

R*4 

ITS 

12 

91* 

314 

317 

319 

322 

325 

326 




337 

340 

341 

K5 

508 

510 



2-00004623 

R*4 

IT6 

12 

9?* 

315 

317 

322 

325 

337 

348 




506 

508 







2-0000462C 

R*4 

IT7 

12 

93* 

313 

323 

338 

506 



2-00034630 

R*4 

118 

12 

94* 







2-00004634 

R*4 

ITS 

12 







2-00004BJ4 

1*4 

nor 

602= 

604(21= 

607 






2-09004724 

1*4 

j 

105* 

527A 







AP-C0000033I 1*4 

JAY 

3 

596A 

61 2A 






AP-OOOOOOOCO 1*4 

JIW4 

3 

434 






AP-000000084 1*4 

K 

3 

23 

33 

594 

595 




2-0000466C 

1*4 

FAY 

28* 

596A 

61 2A 






2-00004800 

1*4 

RDOO 

597* 

599 

600(21= 






2-00004718 

1*4 

F71M4 

101* 

395(2) 

553(21* 






2-00004684 

1*4 

Kti 

571* 

572(4) 

573(4) 

575 

576 

577(3) 

578(3) 

581(3) 




582(2) 

583(3) 

584(2) 

586(2) 

589(5) 




62-000000204 R*4 

RP 

3 

12 

20 

21 

27 

23 

165* 

166* 


A75 
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■n 


VftSMP 





CO 

'vj 

M 

214* 

240* 

242* 

253* 

331 



2-00004654 

6*4 

CM 

14 

21* 

166 

242 

253 

612ft 



2-00004658 

8*4 

m 

14 

22* 

163 

214 

240 




2-G000465C 

R*4 

RPC 

14 

23* 

187 






AP-OOOCOOICI 1*4 

L 

3 

281 

354(2) 

358(2) 

360(2) 

333 

516 

526(2) 




532(2) 








2-00004648 

R*4 

L1M 

13 

368* 

373* 

378* 

383(2) 

387(2) 

331 

332 

2-00004650 

R*4 

(11VR3 

13 

104* 

534(2)* 

564 





2-0000464C 

R*4 

KTOUT 

13 

335(21* 

553 

564 





2-OCOC4SSO 

R*4 

POUT 

333* 

334 







2-00004754 

R*4 

PRECYC 

130A 

136A 







2-00004753 

R*4 

RN? 

13CA 

136ft 

364 

523 





2-C000475C 

R*4 

m 

130A 

136A 







2-O00046F0 

R*4 

SI 

CO* 

153 

162* 

180 

184* 

200 

203 

211* 




222 

234* 

236* 

233 

241 

264 

268= 

271* 




274* 

286 







2-00004SF4 

R*4 

S2 

81* 

163* 

185* 

212* 

237= 

287 



2-000046FB 

R*4 

S3 

82* 

164* 

186* 

a3* 

233= 

283 



2-0000472C 

R*4 

SA1 

128ft 

134A 

223 

236 





2-00004730 

R*4 

342 

128ft 

134ft 

237 






2-00004734 

R*4 

$A2 

128A 

134A 

233 






2-00004784 

R*4 

SGftPP 

180* 

184 







2-00004774 

R*4 

SCP 

158* 

162(2) 







2-000047SC 

R*4 

SGPP 

200* 

203* 

211 






2-00004 ?3C 

R*4 

SSI 

123A 

135A 

140 

162(3) 

273 




2-00004760 

R*4 

SS1Z 

140* 

144 

145 

147 





2-00004740 

R*4 

SS2 

129ft 

135A 

141 . 

143 

144 

163 

185 

212 




554(2) 

558 

364 

563 





2-00004764 

R»4 

SS27D 

141* 








2-00004744 

R*4 

SS? 

123ft 

135ft 

164 

186 

213 




2-000047SC 

R*4 

SSrttf 

223= 

234(2) 







2-00004733 

R*4 

SS64P 

2P2* 

223 

234 






2-00004738 

R*4 

sx 

128ft 

134ft 







2-00004678 

R*4 

T 

43* 

45 

128ft 

123ft 

134ft 

135ft" 

552(2)= 

555 




558 

564 

568 






2-000046FC 

R»4 

T1 

85* 








2-00004684 

R*4 

T2 

63* 

85 

66 

343 

353 




2-00004700 

R*4 

T4 

86* 








2-00004660 

R*4 

TC 

46* 

87 

88(2) 

83 

SO 

31 

32 

93 




34 

13CA 

136A 

167 

183 

213 

244 

277 




552 

555 







2-00004664 

R*4 

76 

77* 

134* 

211 

215(2)* 

216(2)= 




2-C0004770 

M4 

761 

150* 

162 

167(2)* 

168(21* 





2-00004734 

R*4 

T62 

221* 

234 

244(2)= 

243 

216= 




2-00004780 

R*4 

763 

172* 

184 

183(2)* 

133(2)= 





2-00004768 

R*4 

T64 

142 

263* 

267 

268' 

270(2) 

273 

274 

277(2) 

2-0000(634 

R*4 

TF 

53* 

555 







2-000040CC 

R*4 

TFIN 

534* 

536ft 

CCS* 

612ft 





2-030046SC 

R*4 

7G 

66* 

294 

351 







1 r 

. 11 
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A76 
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2-0000474C 

R*4 

TG££ 

130A 

136A 

383 

387 





2-00004S7C 

(1*4 

ill 

45* 

53 

129A 

12M 

134A 

135A 

558 

564(2) 




568- 

• 610* 

61 2A 






2-00004683 

8*4 

TL 

65* 

343 

350 

351 

352 

354 

358 


2-00004884 

8*4 

TOUT 

334* 

335 

336(2) 

337(2) 

400(2) 

402(2) 

404(2) 

430(2) 




432(2) 

434(2) 

436(2) 

433(2) 

440(2) 

442(2) 

444(2) 

446(2) 




448(2) 

452(2) 

454(2) 

456(2) 

458(2) 

460(2) 

462(2) 

464(2) 




466(2) 

468(2) 

470(2) 

474(2) 

476(2) 

478(2) 

480(2) 

482(2) 




484(2). 

486(2) 

488(2) 

430(2) 

492(2) 




2-0000471C 

R*4 

Touirn 

102* 

336(2)* 

564 






2-00004720 

RM 

TOimilN 

103* 

337(2)* 

564 






2-00C048CO 

R*4 

ISM 

581* 

582 







AP-000000180 R*4 

TW 

3 

377 

378 

557 

612A 




2-000048C4 

R*4 

VSOA 

583* 

584 







2-00004830 

R*4 

vro 

563* 

573 

578 

583 





2-0000488C 

R*4 

VTD1 

578 








2-00004750 

R*4 

VIZ 

130A 

136A 







2-000046AC 

R*4 

YO 

61* 

233 

364* 






2-00004SM 

R*4 

Y1 

58* 

344* 

354 

358 

364 




2-C00047F8 

R*4 

yia 

314* 

344 







2-000047TC 

R*4 

Y1B 

315* 

344 







2-000046M 

R*4 

Y2 

53* 

345* 

354 

360 

365 




2-00004800 

R*4 

Y2A 

316* 

345 







2-OOG04S84 

R*4 

Y2B 

317* 

345 







2-000C46A9 

R*4 

Y3 

60* 

346* 







2-00004303 

R*4 

Y3A 

318* 

346 







2-0000430C 

R*4 

Y38 

319* 

346 







2-00004818 

R*4 

YA1 

324* 

352 

354 

353 





2-00004910 

R*4 

YA1A 

322* 

324 







2-00004814 

R*4 

YA18 

323* 

324 







2-C000431C 

R*4 

YA2 

325* 

352 

354 

360 





2-00004820 

R*4 

YA3 

326* 








2-0000482C 

R*4 

YAA1 

333* 

364 







2-00004824 

R*4 

YAA1A 

337* 

339 







2-00004828 

R*4 

YW19 

338*' 

339 







2-00004830 

R*4 

YAA2 

340* 

365 







2-00004834 

P*4 

YAA3 

341* 








2-00004680 

R*4 

YDO 

62* 

365* 







2-00004708 

R*4 

Z 

97* 

140 

146 

233 

241 

279 

394 

385 




386 

512* 

521 






2-0000483C 

R*4 

ZA1 

505* 

512 







2-00004830 

R*4 

ZA2 

506* 

512 







2-00004894 

R*4 

281 

507* 

513 







2-00004833 

R*4 

ZB2 

508* 

513 ' 







2-00004S3C 

R*4 

2C1 

503* 

514 







2-00004830 

R*4 

ZC2 

510* 

514 







2-0000470C 

R*4 

2D 

93* 

141 

384 

385 

513* 

522 

524 

554(2) 




558 

564 

569 







A77 



m2 



l-i/H / / 

; > /• / / 


21-JW-19S5 12:09:35 WX-U FCf.TWN V3.5-62 Ps ? f 17 

19-J»-13S5 21:05:43 FSM:(STEVD(S.SHlPSTUfF]VPA3C.rOR;17 


2-00004710 R‘4 ZOO 
2-00004564 R*4 TA 
2-000046EC R*4 ZWOTDT 


99* 

384 

514* 

26* 

384(2) 

335 

79* 

333 

387 


Mir til Type 

2-00000030 R*4 AC 
AP-000000304 R*4 ATM 
AP-000000284 R*4 AZ 
2-00001E28 R*4 WhB 
2-00001D50 1*4 IS 


AP-C0000034# R*4 PHIO 
AP-0000002C? R*4 PHIZ 
2-00000140 R*4 SLA 
2-00001800 R*4 SLOT 
2-00000000 R*4 SLS 


2-00000073 R*4 SLV 
2-00000460 R*4 91 
2-00003E30 R*4 TIMA 
2-00001400 R*4 TINEE 
2-0C002E90 R*4 Tftl 


2-00003660 R*4 TSO 

2-00001C93 R*4 TWRTR 
2-OOOOlEFO R*4 OW 

2-C00026CO R*4 VSD 

AP-00000024P R*4 HE 

2-00000338 R*4 ZDOL 
2-00000200 R*4 ZDL 
2-00000208 R*4 ZL 


Bytes DiMfisions 

2000 (500) 

40 (10) 

40 (10) 

200 (50) 

200 (50) 


References 


40 (10) 
40 (10) 
206 (50) 
200 (50) 
120 (30) 


200 (50) 
2000 (500) 
2000 (500) 
2000 (500) 
2000 (500) 


2000 (500) 

200 (50) 
2000 (500) 

2000 (500) 


200 (50) 
200 (50) 
200 (50) 


8 

117* 




3 

7 

129A 

135A 


3 

7 

123A 

135A 


9 

31* 




9 

30* 

431(2>* 

433(2)* 

435(2)* 

437(21= 

439(21* 

441(2)* 

443(2)* 

445(2)* 

447(2)* 

449(2)* 

453(2)* 

455(2)* 

457(2)* 

459(2)* 

461(2)* 

463(2)= 

465(2)* 

457(2)* 

469(2)* 

471(2)* 

475(2)* 

477(2)* 

479(2)* 

461(2)* 

483(2)* 

485(2)* 

487(2)- 

489(2)* 

431(2)= 

433(21* 

560 

£04 

607 

3 

- 7 

129A 

135A 


3 

7 

129A 

135A 


£ 

109* 

284(2)= 

288* 


9 

113* 

607* 

612!2)A 


£ 

107* 

282(2)* 

286* 

354(2) 

358 

360 

526 



£ 

108* 

283(2!= 

287* 

532 

7 





10 

586* 

536A 



8 

118* 




10 

572(2)* 

577(2) 

581 

589 

596A 





10 

577* 

581 

582* 

589 

596A 





9 

114* 

539* 

612A 


10 

573(2)* 

578(2) 

533 

583 

526A 





10 

578* 

583 

584* 

589 

596A 





3 

7 

129A 

135A 


6 

112* 

513(2)= 

523* 


£ 

111* 

518(2)* 

522* 

532 

£ 

110* 

358 

360 

333 

517(2)* 

S21* 

526 




Address l»t>el 


A78 


mJ 




n 





r 


WSMP 


0-000003C3 20 

tt 30 

0-00000^34 33 
0-00000480 40 
0-00000365 45 
0-00000501 46 
0-00000506 47 

0-00000622 50 
0-000006F4 65 
It 70 
0-C00C0BS5 74 
0-00000&C3 75 

0-00000662 60 
0-00000017 SO 
0-00000611 S3 
O-OOOOOFCB 100 
0-00001002 105 - 

tt 110 

0-00001225 140 
tt 145 

0- 00001234 150 

1- 00000066 155' 

0-00001451 160 

0- 00001519 165 

1- OOOOOOSF 170' 
0-00001590 175 

0- 00001560 160 

tt 165 
tt 150 
tt 155 

1- 0000005E 200' 


21-Jsii-1535 12:05:35 
15 -Jan-lSflS 21:06:43 


0ax- 11 FOSTIWi V3.5-62 


1 rvwwi' — 

r,..r>rTncriit9Air FfR! 


P»5» 13 
17 


1344 556 

1404 


143 

1714 


144 

1934 


145 

146 

2204 

233 

2364 


235 

2374 


142 

147 

2524 

170 

191 

218 

261 

2854 


370 

3744 


374 

3771 


376 

377 

3834 

401 

4301 


403 

4524 


405 

4744 


407 

451 

473 

516 

5204 


132 

5564 

5524 

5574 

279 

555 

564 

5664 


576 

5014 


575 

579 

5364 

539 

5904 


33 

567 

5514 

552 

5574 



333 6014 
(03 6054 
(06 tO 0 * 
618 6134 


230 


4364 


2794 


( 


t 



{ * 


(UNCTIONS tw SUBROUTINES FXF6POIC6D 


Typt Nat 

FORtMTEJJOS 

FOR$TI1*-6_T_OS 

iNPini 

1NPUT2 

Rt4 KTH4C0S 

R*4 KTHtPTNOOl 
PICT 2 
PLCT3 

TURBU.DK6 


Rtftrtncts 

616 

615 

128 l 04 

125 l 05 

162 1M 

43 527 

612 
596 

130 136 


211 234 


274 


■’S 



A79 





/ 


/; i . 


• f. 

- >ri v 




•/ / 

/ 


■ / * 




it 




0001 


21 -Jan-1335 12:03:33 WX-11 FORTRAN V3.S-62 Pas* l 3 

13-Jan-1335 21:06:43 FSM:tSTEVCW.S8!P$TlfF]VPAlC.FCR:l7 


KEY TO REFEREJ4CE FLAGS 
» - Value Modified 
t - Defining Reference 
A • Actual Argument, possibly »odified 
D - Data Initialization. 

(n) - Muster of occurrences on line 


cams oualifiers 

FORTRAN /US/CRO VPA1C.FDR 

/CHECH* INOOOCNDS , CWERFIOH .NOIMJERFLCH) 

/DE6«G*<t«Sni20LS,TRACE3Aai 

/STffmO*(NDSr:fTAX,faSOURC£_FO(!H) 

S'S®™.'” Mi -.-.IMS «»'»»«•'> 


COMPILATION STATISTICS 

Run Tint; 

Elapsed Time: 

Page Faults: 

Dynaaic Neoory: 


26.45 seconds 
23.33 seconos 
432 

333 pages 


> 

I 

t 

f. 




r 


// 

/, 

•/ 


r 


/ : 




miiuiumHiuiiiimHni) 


2-FEB-13S3 11:37:52.33 
2-FEB-1985 11:37:52.33 
2-FEB-1335 11:37:52.33 


FSMiismwK.SHiPsnrriiM.usii 

FSDO:[STO®S.SHIFSrWFIlW.LlSil 

FSM:IST£VOfi.SHIPSTUFFllM.US|l 


luxwiuuiHHimiimui 






IN4 . L_ I 


5 X 


•M 




2-FCB-1335 11:37:52.33 
2-FEB-1985 11:37:52.33 
2-FE8-1985 11:37:52.93 


F500 : [ ST EVEJ« . 5H 1 PSTIFF ) I IH . U S : 1 
F$D0:( STEVENS. SHIP$TUFFJIN4.1.1S;1 
FSOO: 1ST EVENS. SHI PSTUFF]IN4.US;1 



A81 





2-F«b-19S5 11:27:41 UW-11 F08TIW4 V3.5-62 Page 1 

2-frt>-1585 10:52:12 FSOO : I STEVENS . SHI PSTUFF] IN4 . FOR ; 44 


0001 


SUBROUTINE INFUTKS1 ,S2,S3, t , 1 1 , IFLAG5.SX) 

0002 



0003 


«1PA‘9.0 

0004 


AflPB*31.0 

0003 



0006 


1KS«1.0 

0007 


MISMU2 10.14 

0008 


Sl*0i0 

0009 


S2=0.0 

0010 



0011 


IF (IFIAG5 .EQ. 0) GO TO 20 

0012 



0013 


IF {ISFLAS1 .18. 0) T1HAM 

0014 


MIMA 

001S 


IF (t .GE. A+3.14159/MIS) 60 TO 10 

0016 


TAIM-A 

0017 


Sl«rtlPA*(9.5f0.5*C0S(W«*TAUt3.14159)) 

0018 


S2 3 -FTWi*0 .5ttfIS*SIN WS*TAUt3. 14153) 

0013 


S3*-A<1?A*fl.5*WSit*2*tO$WS4TAUT3. 14139) 

0020 


ISFUAGl=l 

0021 


60 TO 30 

0022 

10 

■ S1=7T1FA 

0023 


S2-0.0 

0024 


S3=0.0 

0025 


60 TO 30 

0026 



0027 

20 

BMl 

0028 


IF (t .GE. 8T3.14159AMS1) CO TO 25 

0023 


TAU=+(T-B) 

0030 


Sl=#mrf«1PB*(0.5M.5*C0$(l«Sl*TAU>) 

0031 


S2=-m°8*0.5*WlSl*SIH(mSl*T«I) 

0032 


S3=-(We*0.5*ttJSltt2tC0S(lNSl*TAU> 

0033 


I SFLAG1=0 

0034 


60 TO 30 

0035 

25 

S1=AWA 

0036 


S2=0 

0037 


S3*0 

0038 


!SFLA61=0 

0033 

30 

RETURN 

0040 


DID 





/ 




P}\ 


INfl/Tl 

mcm sections 

Kw 

0 1CCK 
2 iLOCAL 

Total Space Allocated 


2-Fet>-19B5 11:27:41 
2-Fe6-1965 10:52:12 


WC-11 FGSTR/W V3.5-62 
FSD0:($T£VO<$.$HlPSTUfFJIN4.FG8;44 


Paw 2 


Bytes Attributes 


292 

36 

328 


pic con pa' ia sk* eke 

PIC CCN *a la NOSiS NOEXE 


ro war long 
RD WiT KN2 


entry points 

Address Tvpe Name 
0-00000000 INPUT1 


Mermen 

1 


UWIABLES 


Address Type 

N soe 

Attributes References 





2-00000018 P.*4 

A 

14= 

15 

16 



2-0000.0000 R*4 

em 

3= 

17 

18 

13 

22 

2-00000004 R*4 

(m 

4= 

30 

31 

32 


2-00000020 R*4 

B 

27* 

28 

29 


AP-0000001B* 1*4 

1 FLAGS 

1 

11 




2-00000010 1*4 

ISFLAG1 

13 

20* 

33= 

33* 


AP-00000304* R*4 

SI 

1 

8* 

17= 

22= 

30= 

AP-00000008* R*4 

S2 

1 

9= 

18= 

23= 

31= 

AP-OOOOOOOCf R*4 

S3 

1 

19= 

24= 

32= 

37= 

AP-0000001CI R*4 

SX 

1 





AP-000000101 R*4 

T 

1 

13 

15 

16 

23 

2-0030001C R*4 

TAU 

16* 

17 

18 - 

13 

23= 

AP-00000014J R*4 

TI 

1 

27 




2-00000014 R*4 

T1W 

13= 

14 




2-00000008 R*4 

UiS 

6= 

15 

17 

18(2) 

13(2) 

2-OOOOOOOC R*4 

msi 

7* 

28 

39 

31(2) 

32(2) 


35 


35= 

36* 


29 

30 


32 


uoas 

Address Label 

0-000003A3 10 
O-O0C0CCAF 20 
0-00000116 25 
0 -00000123 30 


References 

15 

11 

28 

21 


221 

27* 

35* 

25 


■M 


39* 


A 83 


* F 

i . | 



INPUT! 


2-Feb-1985 11:27:41 WX-11 FORTRAN V3.5-62 

2*Feb-19B5 10:52:12 FSOO:[STEV0tS.SHIPSTUFF]lN4.FOR;44 


FUCI10NS AN 0 SUBROUTINES REFERENCED 


Type 

Hat- 

References 




RM 

HTHSCOS 

17 

19 

30. 

32 

Rt4 

KTHJS1N 

18 

31 




KEY 10 REFERENCE FLAGS 

* - Value Modified 

♦ - Defining Reference 

A - Actual Argument, possibly codified 
0 - Coil lnttialization 

(n) - Nuuber of occurrences *"i line 


cam® clkufiers 

FORTRAN AIS/CRQ- IN4.F0R 

/O£CR=(N030UNDS,O7£RFL0H,imi0£Ra(H) 

/DEBUG 2 (N0SYM8OLS,TRAC£E3ACK) 

/STANOARD=(fiOSrNTA3<,NOSOURCE_FOR.*i) 

/sucH=(tm£PRoassw,noitauDE,m?) 

/F77 /HOG FLOAT INS /I4 /OPTIMIZE /WARNINGS /WOO LINES /CROSSJEFESENOE /MTttCKINE.COOE /CtNTINlWTICHS 2 19 


CCKPILAT1CH STATISTICS 

Run Time: 2.16 seconds 

Elapsed Time: 2.B4 seconds 

Page Faults: 127 

Dynamic Memory: 117 pages 



s. ’" v 



H ' linnmiiimmumiiii 

miiiiniiiiiiH i inmiiHrt 

wiimuinmiiiiuiiumii 


2-FF8-1385 11:32:46.8? rSOO:ISTEVDfS. SHIPSTIFFJGENSM.LIS;! HlUlnm > »imimiiunii 

2-FFB-1S85 11:32:46.87 FSM:[STEV£NS.SH!PS71JFFJ(»BH.US:1 t tiiiiiiimmi.uiuMDixi 

2-fE8-1335 11:32:46.87 FSSD:'STEm.SMPSUFF]G£N3UI$;l iwuimaunininmiim 


GENSM -LIS;! 


i r tttHtik H -H r t t i ****** 



2-FE8-1885 11:32:46.87 FSM;|STEVB«.SHIPSTlfFJGEJ«UISi] 

2-FE8-1JS5 11:32:46.87 FS00:(3Ta£M.SMFSW]GEN91.US:l 

2-FEB-1885 11:32:46.87 FSM:rCTEVC«.SK;p$Tl*T]GEN9UIS;l 


*«****«4«UihHtm***+*m* 

WrtWtKW HHUIIIIlUimn 


A85 


' KZJ' 

4 

I 

I 

'I 






2-f»fc-l?33 11:23:19 WA-11 fCCTVN V3.3-£2 Pi$ f 1 

lO-fcc-1334 17:13:40 f3a):tSTittN$.$HlP$TUff)G£mfOS;3 


r- 


0001 


SU8P.CUTM INPUT2I $1 ,S.-.,S3, t,t! ,CN,R£,AZ,FtiIZ,ATH,PHJ9) 

0002 


DlMCNS!O4lI(10).A2(10),fHI2(lC),ATK(10),FW0(10) 

CC03 



0004 


If (1 .GT.lI) CO TO 10 

0003 



ooot 



0307 



0008 

C 

444WCCH fmSf OKU GEMJATED44* 

000) 



0019 


TW » SECNOS(O.C) 

0011 


H * NINTOT1) 

0012 


»i » n«T(«) 

0013 


ST « TT1 - X31 

0014 


S « N001S1.2.) 

0515 


IF ( S .CO. 0.0 ) THIN 

00U 


M1« »!+ 1. 

0012 


ETC IF 

0013 


M * IFIXOffl) 

001! 


XB * -2147433F13.0 * RT 

0029 


S > X3 / 2. 

0921 


IS * IfIXIS) 

0922 


XS * FlCftT(IS) 

0923 


S • S -XS 

0024 


IF ( S J£. 0.0 ) T1CN 

0025 


X8 * X8 4 1. 

cars 


DO If 

002? 


IB * IFIX(XB) 

0029 


N « IB -N 

0023 


Y * RRNltir - 0.5 

0033 


D4*«.2S32*Y 

0031 


04=4.0 

0032 



0033 

c 

444 SNIP 1D4.T (FT,)444 

0034 



C03S 

10 

Hl=AZ(l)4COS(«fl)4t4P«lZ(lHEN) 

COSE 


H2=A2(2)4COS(«(2!4ttPHI2(2HEN) 

003? 


H3=A2( 3)*COSO,t( 3)4Hf«IZ( 3HD4) 

0033 


H4=AZ(4)*IOS(kI-(4)*t4PHIZ(4HEN) 

0039 


KS=AZ(3>*CO$(irt5)4ttf!N!Z(3HEH» 

0044 


H£=A2(£)4CCS(kt(E)4itFHIZ(F)tEN) 

0041 



0042 

c 

444 SNIP PITCH (RAOJ444 

0043 



0044 


Pl=ATH(l)*CK(WC(l)4t4FHIO(lHOO 

0045 


P2=ATH(2)*CO$(l£(2)4t*PH!0(2HEN) 

0940 


P3=ATH(3)4CCS!lf(3)*t+PHI0(3)4EN) 

004? 


P4=ATH(4)4CCS(k£(4)»t(PHlO(4)tCil 

0043 


P5=ATH(5)4CC$(Rtt3)4t»FH10!5»EN) 

0049 


PS=ATH(£)4C0$(WC(6)*l4FHIC(£)KN) 

0033 



0031 



0032 

c 

4HHW.I AT WCING FfiWt* 

0053 



0054 


Sl=HlfH24H3*H44N34H£4(]f0.34SIN(FltP24P2IP44PS4F6)) 

0033 

c 

SI— SI 

C03i 



0057 

c 

444 snip Hc«.t mam in/sim 


A86 



' '\ 


2-frt-lH5 11:29:19 IttX-ll TtMPM V3.5-62 
10-0«c*lJ94 17:13:40 fS»:[SmtW.SHIPmfF]G0fiM.ra;3 

0858 

0359 K)l*-«<l)tWEU)*SIN<t£(l)4tmfim)«3<> 

COM roZ*-«(2)*UEfl)*S!WW£(2)*ttf>HmiMOO 

C3£l K>3*-*Z(3)*«:<3]»SINOf(3)Ut?HIZ(3]tO<) 

0002 HM=41Z<4)*H£<4J*SIN<HE(4!«t+PWZ(4HDO 

080 K>5-«(5]»*<3)*$If«kC< 5) • ttPHIItJUDO 

COM HMWa6)«£(6)*SINU£(£)4ttFH!Z({)tOi) 

0303 

oo« c «*smf> mot velocity «wo/si«* 

0507 

0308 PS:«4»TH(l)*rt<l)*SIN(«(l)*ttPHWHEN) 

0509 P02='-AIK(2)*«(2)*StN(UE(2)»(tPHI0(2n£N) 

0028 PD3»iVTH(3)*WC(3)*$IN(HE(3)*ttFXIC(3HPO 

0871 PM*-ArH(4)#£(4>*SIN(«EM)*4tPH!0t4HEN> 

0872 P05«-ArH(3)»W(3)«S!N(W(3)'>r*PHI!!<5>-£N) 

0873 fK»T3TH(0)«tf(6)*SiN(H£(6!*It?H!0(il‘0« 

0874 

0075 C trtSHjp utQlW P « WWE VELOCITY ((7/3!** 

0370 

08 77 S»=h01««24H>34HMW»W>S 

0079 $2M8B .8*SIN( P0HP02PP03WI!4tP05tFC*) 

0079 

0088 S2=$2MS29 
0 0 83 C S2* S2 
00?2 

00S3 C tttSHIP KALE ACCaEMTItN <PT/S*42)t« 

COM 

0085 

0080 M»l*-WEcn«2*A2(l)*CQS(8ff(l)*t+P«3rtlHEN) 

0037 K3D2«4£(2)H24AZ(2)*C£!Sri£(2)»t+PH!2i2)KN) 

0089 H»3=4£(3)«2*A2(3)*CC$(«(3)4ttmi:(3)»EN) 

0089 K»4v4(E(4)M2*AZ(4)*OS(l«4)*t4FHIZ(4HOO 

03M rc03*-+S(5)M2»A7(51»CSS(«:(3)*t4PHi:(5)4£N) 

0051 K06=4iE(6)H2*A2(£)*C)S(t£<6)»ttPHi2(8)l'50 

0052 

0033 C SKIP PITCH ACtEUMTICN PAa/SH2 

0094 

0095 FTt'i=-H£(l)M2»9TH(l)tCOS»fE(l)«t4FWC!«m[N) 

0058 PM2=-HtJ.-!4»WTH(2)*£OSO«£(2)*t4PH;W2)*£X) 

0097 PK3>H»t(3)*i2»ArH(3)*C0S(«(3)*HPH!0(3!4Di) 

0098 Pft,3»-4fE(4)«2*9TH(4)*C0S(l4E(4)*t*PWS(4I«J4( 

0099 Ft«5=-HE(5)M2*AIH(3)»CaS(H£(5)*I4PHI8(5»EN) 

0108 PC©S*4C(6)M?tATH(S>4C0S<WE(£)ttlPHiEt8)ID4> 

0101 

0102 C ***SHIP L8M1JNS m ISAE ACCELERATION (TT/8H2)*** 

0103 

0104 S»»H0QHKl02tK)034HDC44IO054H»S 

0103 S^=160.e*Sm(P0D14POO2+PCO3fPOC4tPCe5+fCOi) 

0108 

0107 S3=S3A+S3B 

0108 C £3— S3 
0103 

0118 HCTISN 

0111 BO 


Pi9» 2 


A87 


IKTUT 2 


/ 


/ 


program nitus 

Km 

0 tccoe 

1 4PMTA 

Tsl»l So set Allscsttd 


2-frtrl?33 11:28:15 VAX-1 l FORTRAN V3.5-62 

lO-Ctc-1534 17:13:48 FSOO:(STM<S.SH!P$n!fFlG£X9UOR;9 


Evtn Attnbutn 


1667 

pic cm rel La 

SHR at 

R 0 MMT LOS 

8 

PIC CSN REL La 

SWNCexr 

R 0 MIST LOS 

388 

PIC C£N REL La 

NOSH! NQ 6 XE 

RO 1 ST LCNG 

2043 





DffKr POINTS 


Mdrrit T*pt Kw Rpftruicn 

(-00083003 IffUTZ 1 

W 0 UAa.ES 

AOdrtsi Type taw Attribute Rpftrmcti 


AP-OOOOOOIEO R* 4 - 

04 

1 

30 * 




43 

44 




63 

61 




71 

72 




SI 

33 

2 H 5 C 0 SC 323 

R *4 

ra 

33 * 

34 

2 - 0003002 C 

R *4 

H 2 

36 * 

54 

2 - 00 CCO 030 

R *4 

H 3 

37 * 

54 

2-00000074 

R *4 

B 4 

38 * 

54 

2-00050039 

R *4 

HS 

3 S* 

54 

2 -OOCC 033 C 

R *1 

W 

40 * 

34 

2-00000059 

R *4 

101 

35 * 

77 

2 -C 000005 C 

R *4 

W 2 

60 * 

77 

2 - 0003 C 0 EO 

R *4 

103 

61 * 

77 

2-00000004 

R *4 

104 

62 * 

77 

2-00000068 

R *4 

H 05 

63 * 

77 

2 - 0000306 C 

R ‘4 

106 

64 * 

77 

2-00000055 

R *4 

HCOl 

86 * 

104 

2-00003034 

R *4 

1002 

87 * 

104 

2-00030053 

R »4 

1 CC 3 

68 * 

104 

2 - 0000005 C 

Ri 4 

Et >4 

BJ* 

104 

2 -OOCOOOA 3 

R *4 

1003 

? 0 * 

104 

2 -OOOCOOA 4 

RH 

1006 

31 * 

104 

2 -OOCS 3 C<S' 

1*4 

IB 

27 * 

28 

2-00033018 

1*4 

IS 

21 * 

22 

2 - 030 C 0004 

1*4 

It 

11 * 

12 

2-00030040 

R *4 

PI 

44 * 

34 

2 -OOOC 5344 

R *4 

P 2 

45 * 

54 

2-00000143 

R *4 

P 3 

46 * 

54 


R *4 

M 

47 * 

34 


31 * 

33 

36 

37 

33 

39 

43 

46 

47 

48 

43 

39 

62 

63 

64 

68 

69 

70 

73 

86 

87 

63 

89 

30 

36 

57 

38 

33 

100 



16* 26(21* 2SA 


- A- 


P*?t 3 


A88 



■X:':x 

— j . i 

, / ’ ■' 

A 

\ « 
' \ 

r ' 

■ \ 


V 



> 

. n ; 

t;' 

>' 











, i 

* 

( 

v * 

INPUT2 






2-F>b-19G5 I! 

1:23:19 

WW-ll FORTRAN V3.5-62 










ID-Dk-1584 17:13:40 

f SOO :( STBKNS.SHl PSTurF)G0i$M.fOR 

\ 


2-00000050 

RM 

P5 


49« 

54 







2-00000054 

R*4 

PE 


49* 

54 





\ 


2-00000070 

RM 

POl 


E8« 

78 







2-00000074 

RM 

P02 


69* 

78 





• ” 


2-00000073 

RM 

FD3 


70* 

78 







2-0000007C 

KM 

P04 


71* 

78 







2-00000080 

RM 

P05 


72* 

78 







2-00000064 

RM 

PD6 


73* 

78 







2-OOOOOOP.8 

RM 

POM 


95* 

105 





• ...-V . 


2-OOOOOOAC 

RM 

P002 


96* 

105 







2-ooooooeo 

RM 

PM3 


97* 

105 







2-00000084 

RM 

P004 


98= 

105 







2-00000083 

RM- 

POOS 


59= 

105 







2-0000008C 

RM 

POOS 


100* 

105 







2— OCOOOOOC 

RM 

RT 


13* 

19 







2-00000010 

RM 

S 


14* 

15 

20* 

21 

23(21* 24 




AP-OOOOOOQ4# R*4 

SI 


1 

54* 







AP-00000008# R*4 

S2 


1 

80* 





\ 


2-00000088 

RM 

S2A 


77* 

80 







2-0000008C 

RM 

S2B 


78* 

80 





. , ■ 


AP-OOOOOOOCf RM 

S3 


1 

107* 







2-00000000 


S» 


104* 

107 





1 

\ 

2-OOOOOOC4 


S38 


105* 

107 





r^" 

\ 

i 

tfM)0000010# RM 

T 


1 

4 

35 

36 

37 38 39 

40 


1 





44 

45 

46 

47 

48 49 59 

60 







61 

62 

63 

64 

63 69 70 

71 


1 ~ 





72 

73 

86 

87 

88 89 33 

91 

-■» •. 

1 # 





95 

96 

97 

98 

99 100 




AP-000000144 R*4 

TI 


1 

4 






T “ 

2-00000000 

RM 

m 


10* 

11 

13 





j 

2-00000014 

RM 

X8 


19* 

20 

25(21* 

27 




1. 

2-00000008 

RM 

Ml 


12* 

13 

14 

16(21* 

18 




2-OOOOCOlC 

RM 

XS 


22* 

23 






i 

* 

2-00000024 

R*4 

Y 


29* 

30 






I 


X- 



ARRAYS 

Addrtss Tvpt Nane 
AP -00000028# RM ATH 


AP-00000020# RM AZ 


Attributes 


8ytes Omeosiom References 


40 (10) 

1 

2 

<4 

45 

46 


47 

43 

49 

68 

69 


70 

71 

72 

73 

95 


96 

97 

98 

99 

100 

40 (101 

1 

2 

35 

36 

37 


33 

33 

40 

59 

60 


61 

62 

63 

64 

B6 


87 

83 

69 

90 

91 

40 (101 

1 

2 

44 

45 

46 


47 

49 

49 

6v 

69 


70 

71 

72 

73 

95 



A89 



\- 





AP-0000002C# RM PHIO 





IWW2 




2-f>b-1985 11:28:19 

UW-il FORTRAN V3.5-62 






l#-Oec-1984 17:13:40 

FSD0:|STEVEN3.SH1PST(JFF1GEN31.FCR;9 





96 

97 

98 

99 

100. 

AP-OOCOOC24# RM 

PHIZ 

40 

(10) 

1 

2 

35 

36 

37 





33 

39 

40 

59 

60 





61 

62 

63 

64 

86 





87 

88 

89 

90 

91 

AM000001C4 RM 

HE 

40 

(10) 

1 

2 

35 

36 

37 





38 

39 

40 

44 

45 





46 

47 

48 

49 

59(2) 





60(2) 

61(2) 

62(2) 

63(2) 

64(2) 





68(2) 

69(2) 

70(2) 

71(2) 

72(2) 





73(2) 

86(2) 

87(2) 

88(2) 

89(2) 





90(2) 

91(2) 

95(2) 

96(2) 

97(2) 





98(2) 

99(2) 

100(2) 




/• 


LABELS 

Mdress Label 


Kef trencts 


o-oooooocf 10 


350 


FlUCTICNS WO SUBROUTINES RUERDICED 



Type No* Rtf unices 


RM rORSSECK* 10 

RM HTHMKJO 14 


RM KIHSCOS 

35 

36 

37 

38 

39 

40 

44 

45 


46 

47 

43 

49 

86 

87 

88 

89 


90 

91 

95 

36 

37 

38 

99 

100 

RM HTH4RWCM 

29 








RM KTHSSIN 

54 

59 

60 

61 

62 

63 

64 

68 


69 

70 

71 

72 

73 

78 

105 






mJ. 








CWPILATICN STATISTICS 

Run Tilt; 8.42 seconds 

Elipsed Tiw: 9-33 seconds 

Pj?t Faults: 244 

Dynanic Hano'y: 178 pages 










A- 



iimHhm t iumi t moaurt 


2-FEB-19S5 Jl:«r26.97- FSD0:(STEVENS.SHIPSIlff]TIW0.LI3il 

2-FEB-1985 11:40:26.37 FSOO:[SIEVDfi.SHlPSTUFF]IfWNO.LIS:l 

2-FE3-1385 11:40:26.97 FS00:[STE\'ES’S.SHIPSTUFF1TF-««).L1S!1 






?*’ 

V. 

V 


r- 


TRAN O . l_ I S ; 3 - 



ammaamiuatw aaa 


2-FEB-19S5 11:40:26.97 FSOO:lSTEUOB.SH!PSTUFF]rfANO.LISil 
2-FE8-19S5 11:90:26.57 FSOO : [ STE'.’ENS.SHl P5TUF F1TIW0.LI S;1 

2-FE8-1S3S 11:40:26.97 FSOO : t STEVO. 1 ? .SHI F>STUFF JTRAf.t). L IS;1 


ltw«4W rt ut t m tytwm: 

tttitttHrtiu mm i nHtHn ** 

Hw nuiii u wnrt utwitt 


e 


A 92 








I 



2-feb-1935 11:28:01 UW-11 FORTRAN V3.5-62 Pag* 1 

2-Ftb-I985 10:48:16 rSDO : ( STEVENS . SHI PSTLf/ ) ; 56 


0001 

C2345678901234567890...*te. 

0002 


SUBROUTINE RR8ULB4CEOKE,TC,Vr2,PR£CYC,RNP,RNQ) 

0003 


104=7.79 

0004 


Slfftt=0.032 

0005 



0006 


1004=2.8 >1.8 >2.8 

0007 


SIOW4=0.0003 >0.005 .IJ003 

0008 


1040=0.0 fl.fi ! 2.6 

0009 


SiffttQ=0.0 !24.0 <0.6 

0010 



0011 

5 

PR£CrC=PRECTCH 

0012 


tm * sEcwsto.o) 

0013 


H » HINTON) 

0014 


)N * FUMT(M) 

0015 


RT * TH - >N 

0016 


S * MOWN, 2.) 

0017 


IE ( S .EQ. 0.0 ) D04 

0018 


XH * XH 4- 1. 

0019 


DC IF 

0020 


H * IFIXOH) 

0021 


XB * -2147483648.0 * RT 

0022 


S = X8 / 2. 

0023 


IS = lEIX(S) 

0024 


XS * aOAT(IS) 

0025 


S = S - XS 

0026 


IF ( S .NE. 0.0 ) THEN 

0027 


X8 = XB 1 1. . 

0028 


D4D IF 

0029 


16 = IFIX(XB) 

0030 


M = IS - M 

0031 


FLAG=0 .0 

0032 


Y8=0.0 

0033 


r=o.o 

0034 


DO 10 1=1,12 

0035 


YA=RA»4(m 

0036 


y=wa 

0037 

10 

CONTINUE ' 

0038 


Y8=(Y-6.0) 

0039 


VJ70=-AOWT2+SI ffVW(2*104)rt24 (1/00*0 .5) )*Y8 

0040 


VTZ=(VT71TC*VJZD) 

0041 


TGEE=VT7 

0042 



0043 


Y8=0.0 

0044 


Y=fl.O 

0045 



0046 


DO 15 1=1,12 

0047 


YA=RA«(H) 

0048 


Y=Y!YA 

0049 

15 

COfTINUE 

0050 


YB=(Y-6.0) 

0051 


17400= HOOT* PMHSK0YY4*(2N044) **2* (1/(TC**0. 5) )*Y8 

0052 


RJ4O=(W4O+TC*SN00) 

0053 


RJP=Pf40*0.066 

0054 



0055 


YB=0.0 

0056 


Y=O,0 

0057 





A93 



/ 


/ / 
;/ 




i * 

c 

TURBIUME 2-Feb-1385 11:28:01 VWC-11 FORTRM V3.5-62 Pine 

2-Feb-1985 10:48:1 b FSMs(STEVENS.$HIP$TUFF)TR»IO.FOR:36 

0058 DO 20 1=1,12 

0055 YA=R/W(H) 

' 0060 Y*YfYA 

> 0061 20 CONTINUE 

0062 TB*(Y-6.0) 

0063 SNS0«-Wf«RNI»SIGHA0*(2<Hil)«3t(l/(TCH0.5))*Y8 

0064 RNQ=CRNO+TC*RNQO> 

0065 RNO--IW0.066 

0066 

0067 

0068 IF (PRECYC .11. 100) GOTO 5 

0065 50 ' RETURN 

0070 END 


PROGRAM SECTIONS 


Naoe 

Bytes 

Attributes 



0 WOOE 

437 

PIC C£N REL LCL 

SHR EXE 

RO NCWtT LONG 

1 4P0AIA 

8 

PIC cm REL La 

SKR NOEXE 

RD NOsST l(N5 

2 4L0CAL 

124 

PIC CON REL La NOSHR NOEXE 

RO « T LOJG 

Total Space Allocated 

569 





ENTRY POINTS 

Address Type Hast References 

0-00000000 TURBULENCE 2 

VARIABLES 


Mdttss Type Nait- Attributes References 


2-0000003C Re 4 

FLAG 

31= 

2-00000048 1*4 

1 

34= 

2-00000033 1*4 

IB 

29= 

2-00000030 1*4 

IS 

23* 

2-0000001C 1*4 

H 

13* 

AP-00000010# R*4 

PRECYC 

2 

2-C0000059 R*4 

RNO 

51 

2-00000054 R*4 

RNOO 

51= 

AP -000000144 R*4 

RNP 

2 

AP-00000018# R*4 

RNO 

2 

2-0000005C R*4 

R.NQO 

63= 

2-00000024 R*4 

RT 

15= 

2-00000023 R*4 

S 

16= 

2-00000004 R*4 

SItW 

4= 

2-OOOOOOOC R*4 

SIGMON 

7* 


46= 58= 

30 

24 

14 20= 30(2)= 35A 47A 594 

11 ( 2 )= 68 

52(2)= 53 

52 
53= 

63 64(2)* 65(2)* 

64 
21 

17 22* 23 25(2)’ 26 

33 

51 



A94 


TUR6U.ENCE 



/ 




I 



2-40000014 11*4 SI3W1 
AP-00000008* R*4 TC 
AP-00000004J R*4 TGEE 
2-0000001 S R*4 m 
AP-OOOCOOOC# R*4 VIZ 

2-00000050 R*4 TOD 
2-00000000 R*4 W 
2-00000003 R*4 P4f| 
2-00000010 R*4 WO 
2-0000002C R*4 XS 

2-00000020 R*4 
2-00000034 R*4 XS 
2-00000044’ R*4 Y 

2-0000004C R*4 YA 
2-00000040 R*4 Y8 


2-Feb-1985 11:28:01 4.WC-11 FCRTRW V3.5-62 Pa?* 3 

2-Feb-1905 10:48:16 fSOOllSTEVENS.S.HIPSTUFf jrRW0.F0R;5S 


3 * 

£3 







2 

2 

33 

41* 

40 

51 

52 

63 

64 


12* 

13 

15 






2 

39 

40(2)« 






39* 

40 


1 





3* 

39(2) 







6* 

8* 

51(2) 

63(2) 






21* 

22 

27(2)* 

23 





14* 

24* 

15 

25 

16 

18(2)= 

20 




33* 

62 

36(2)* 

33 

44* 

48(2)* 

50 

56* 

60(2)» 

35* 

36 

47* 

48 

59* 

60 



32* 

£3 

33* 

33 

43* 

50* 

51 

55* 

62* 


LAEELS 


Address Label References 


0-d000002£ 

5 

111 

£3 

H 

10 

34 

37#- 

** 

IS 

46 

491 

« 

20 

58 

£11 

** 

50 

£31 



FUNCTIONS tttD SUBRdfTlNES REFERQJCED 


Type Name 

References 

R*4 FORSSEOiDS 

12 

R*4 HTWAftOO 

16 

R*4 NT KS RANDOM 

35 




K£Y TO REFERENCE FLAGS I 

* - Value Modified | 

♦ - Definin? Reference | 

A - Actual Aryjnsent, possibly modified I 

D - Data Initialiiation | 

(it) - Number of occurrences on line I 


A95 






2-Feb-1935 10:43:15 


rS00:l5TEWEMS.SHlPSTUFFJTRPN0.F0R : « 




catwa oualifius 
F 0RIIW4 /llS/CttO TMNO.fOR 


/C xa-.(MW*S,WTLW,mM*SlW 
/0€6UG= (NOSttlSCt. S . TRACEBf ! 4C10 

FsSsS 



COIPILATHH STATISTICS 


Hun Ti»t: 
Elapwd Tisi: 
Pa« Faults: 
Dyn»:c He»«ry: 


3.09 ttMcds 
3.84 iKondi 
133 

13? pas* 






t ktttkttnt U ti i lttttH U kl K t 

M tttttt H tU iiHit t UnHI t t 


2-FEB-1S85 11:54:35.32 FSM:[ST£Vafi.SHJPSTUFF)MPI.T2.LIS;l MHMrt t Mntmtinttittt 

2-FEB-1S85 11:54:35.32 FS00:(STEVENS.SHIPSTUFF1MPLT2.US;1 tHM t HttitMttmttitn ** 

2-FEB-1935 11:54:35.32 FSD0:(STCVQts:SHlPSnrFmT2.LIS;l iimiiimmiiiimiBH i H t 




n mmmm tWH iww m 


2-FEB-1985 11:54:35.32 FSDO : J STEVEWS.SH1 PSTUFF JMPLT2.U S;I 

2-FEB-19S5 11:84:35.32 FSOMSTEVO.'S.$NlPSTlTFmT2.US;l 

2-FEB-1SS5 11:54:35.32 FSM:[$TEVENS.SNIPSUrFlKaT2.US:l 




A97 



V 


! 


t 

V ' 




IN 


/ 


C 


r 


0001 

0002 

0003 

0000 

0003 

0000 

0007 

0003 

0009 

0010 
0011 
0012 
0013 
0010 
0013 
0010 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
0020 


0023 

0020 

0027 

0028 

0029 

0030 

0031 

0032 

0033 
0030 
0033 
0030 
0037 
0033 

0039 

0040 
?041 
0002 
0003 
0000 
0005 
0040 
0007 
0048 
0009 

0030 

0031 

0032 
C053 
0050 

0033 
0050 
0057 


02345078901230307090. . .etc. 

StESOUTINf ftOT2(X,Y,l(,Xt .Jff, 101 , JAr,MY,TM,2M,aC,lt? ISC) 
OI«NSI(MX<IOI),YdOI),*dDI) .".U-l.W.ISE) 

REAL U> 

CHARACTERtO, BUT 
cmiwaDi*8,6irf 
charactej!*i tAoa(20) 

CtWWCTEWl 1012(08) | LEL3(59) 
omaa*43 i/«a2 
CHARACTEW59 LASEL3 
EOURWUNCE (tfEEL2,LBL2(l)) 

EOUICWLENCE (LA££L(1),SUF) 

E0U1U4LENCE OA6£L(10).$ft4C£) 

Eoumocc oAsntni.pjfF) 

EOUIWliTO (U8a3.LBUil)) 


2-feb-lJ85 11:48:23 
29-3in-lS65 18:00:53 


WW-I1 FORTRAN '73.5-02 Pjje 1 

fSM:(STn'0fi.SHlPSIUrnfPLT2.f(«:43 


NO. kUHS: SEA 

2M: TAUCKS: 


MTA LAB£L2/'GR0CP NO.: 

1 STATE: V 

data LABav-r/w mx : 

1 KP: V 

“TOOE (3,5,LBL2(12)),IAY 
O'CSOE (3,5,La2(29))XAY 
ETCOOE (3,5 1 ia2(46))ISE 

5 Fosmr d3) 

DCKJE (4,10 1 L£S.3(10)ITM 
tWCOOC (4,10,LSL3(23))ZM 
EHCCCE (4,10,LEL3(40))£LC 
MXK (6,i5,LEL3(54))XP 

10 FORMAT (F4.2) 

15 FORmr (F 0 . 4 ) 

CALL TIHE<eOFF) 

CALL OATE(SUF) ' 

space*' 


10 * 101-1 

CAa ctwfis 

<ACL PAGE (II. ,8.5) 

CAU XNAME ('T/H RATIC*',100) 

CALI mt ('TIME (SECll',100) 

CAU ww ('PERCENT CF TIME!', 100, 

CALL )W€ ('HEIGHT (FT)*', 100) 

CAU ASEA2D (9., 5. 75) 

CALL HEAOIN ('T/M RATIO STATISTICS*' 29 1 4) 

CAU BASIN (3EREFCWEELJ.20.0.7 4) ’ 

CALL HEAOIN WIF(LAE£L2),48,0.7,i) 

CALL HEAOIN <*R£F(IAE£L3),59,0 7 4) 

m^^r ,CTwffa:pwA/c ™'^i..i) 

CALL GRAF (XI,2..XF,-10.,3. l 20d 
FOLLCMIW GRAPH CALL IS FCR WUC.FCR 
CALL GRAF (XI,.02.XF,0.,.1,0.4) 

CAU. GRJO (1,1) 

CALL MARX£R(1£) 


t 


i ■ 

u 

: "*> 

• ) 


A98 



•V • .. . 


■ \ ^ 

c- A" ^ 

\< • 


2-Feb-1985 11:48:39 
23-Jai»-1985 18:03:55 


i.WC-U FORTPj=« V3.5-62 Pace 

FSCff:CSTEW£NS.SHlPSIUfT]«PLT2.F0R;43 


CALL. CURVE (X.Y.JD.l) 
CALL HARKERIC) 

CALL CURVE (X,R,I0,1C> 
CALL ENDPL(O) 

CALL MNEPL 
PRINTVID1*',I0J 
RETURN 
DC 


PROGRArt SECT11HS 


Attributes 


0 SCttE 

548 

PIC CON REL La 

SHR EXE 

RO NCJtRT LONG 

1 SPOATA 

141 

PIC CCN REL la 

SHRNOEXE 

RD N34RT LCM3 

21LCCAL 

528 

pic cm Ra La noshr ncexe 

RO 1ST LING 

Tttal Space Allocated 

1217 





ENTRY POINTS 


Address Type Nam 


References 


0-00006000 PL0T2 


Mims Type Naae 


Attributes References 


2-00000038 CHAR EOF 

E0UIV 

5 

12 

35A 


. 2-00000047 CHAR 8Uff 

EQU1V 

6 

14 

34A 


AP-0000C02C# R*4 ELC 


2 

29 



1 2-00000080 1*4 ID 


39* 

58A 



AP-COCOOC184 1*4 IDl 


2 

3(3) 

39 

63 

j AP-000000341 1*4 ISE 


2 

24 



! AP-C000001C4 1*4 JAY 


2 

22 



■ AP-CC0000204 1*4 MY 


2 

23 



AP H30C000304 R*4 RP 


2 

4 

30 


: 2-0000004P CHAR LABEL2 

4 

i ■ . ■ 

EtSJIV 

9 

11 

183 

49A 

1 

! 2-00000000 CHAR LABaj 

E0U1V 

10 

15 

209 

5CA 

2-00000044 R*4 SPACE 

EOUIV 

13 

36* 



AP -000000244 R*4 TH 


2 

27 



J AP-C800C0144 R*4 Xf 


2 

55A 



j AP-000030104 R»4 XI 


2 

55A 



AP-000000284 R*4 34 


2 

28 




A99 



P10T2 


2-f«(>-1985 11(48:33 um- 11 FC*rw< V3.5-»2 P*j» 

29-J*«-l?8J 18:03:15 fSOOUSTCVOtS.SHIttTirr ;tflT2.rM;43 


tom 


AMrm Tr»* Ke» 


AltriCdio Dttti Dteeautat £»f«rt*c*» 


2*8(804933 cxa U30. 

count 

20 (23) 

2-4303:049 CHS* 1&.2 

court 

40 (48) 

2-45500033 00* LEI 3 

ccun 

59 (51) 

AM9MC04C4H4 * 


M (!) 

tf-OOKOCOW £14 X 


*♦ (*) 

W*8S»0«8O *14 r 


H (1) 


usaj 



Mmi 

lt£«l 

Wuncn. 

1*8988(982. 

5* 

22 23 

1*80400025 

It' 

27 20 

1-852SSJ9 

15' 

38 339 


fUCTlWS *© Fjmi!t£S tOTKJCO 


T>f* Hct 

liftmen 

A3XA20 

4( 

ttms 

<9 

CUSVC 

53 

CC£Pl 

a 

OCA 

a 

rOMCrtTt T OS 

35 

ratine t os 

34 

ttff 

55 

BUD 

58 

toow 

47 <8 

renn 

57 

PAEC 

41 

TVWW 

52 

Ktti 

42 

YWrt 

44 


7 

8 
8 

38 

2 

2 

2 


2 < 2 * 

» :.i 


43 19 


12 

13 

14 

48A 

11 

22 

23 

24 

15 

27 

28 

21 


1 

] 508 

3 3» 


«Y 10 tfTCKDtt CLAES | 
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B1 



A. TEST VARIABLES 


The following table gives a listing of the variables used in 
the test matrix. 

Table B1 s Fixed-Base Simulation Computer Variables 
and Related Variables and Values 


LIST Of VARIABLES USED I!f flXEO-BASE SIMULATION 

TEST VARIABLES 



Sea State 

0 , 4 , 6 


BUD 

i. j 


T/U 

1.1, 1.07, 1.05, 1.03, 1.01 1 

Damping 

0.0, -0.2, -0.4 aec“ 


Engine Lag 

0.3. 0.7 .k 


HEl 

MULATTO ITEMS 

C0m*T7* VAST ABLE 

Engine Lag 



T tng (a«c) 
0.3 


TAL’WC (i.c) 
0.3 

0.7 


0.7 

Vertical 

Velocity 

Duping 

5* 

A/C CAIK 

T (sec) (ft/»*c 2 )/°PLA 

CHDOT 

°PLA/(ft/*ec) 

0.0 

- 0.45 

0.0 

-0.2 

5.0 - 0.45 

0.43 

-0.4 

2.5 - 0.45 

0.87 

mUST/VTTCHT 



T/W 


WAIT* (lb) 
S.S.O S.S.4-6 

1. 00 

• ' 

10574 18739 

l.Ot 


18786 18554 

1,03 


18422 18193 

1.05 


non 178*7 

1.0? 


17733 17513 

1.10 


177*9 17035 

* Different weights are used, due to. the changing eeaa vlnd as 
a function of a as state, in order to tlLalnaea the effects of 
the aerodynaslc forces. 
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The mean aerodynamic force for a given sea state was determined 


by: 


FAZN = 


WAIT 

HKTW 


FGS 


(B1) 


Where: FAZN - Mean aerodynamic force 

WAIT - Aircraft weight 
HKTW - Scale factor 
FGS - Gross thrust 


The following table shows the values used for the above variables 


as a function of sea 3tate. 


Table B2: Variable Values Used in Correcting 

for the Mean Aerodynamic Forces 


TRIM FOR T/W - 1.0 

Sea State 

HKTW 

WAIT 

FGS 

FAZH 

0 

0.99045 

17035 

17530.3 

-330.25 

4 

0.99045 

17035 

17353.3 

-163.23 

5 

0.98734 

17249 

17509.7 

- 39.48 


B. HUD DYNAMICS 


The following equations describe the hover point dynamics used 
in the HUD system: 


e = (s.f.) K (X ) - K ,t(T X + 

x pos position ,f el h 


+ V ) + k ,<5, _ J 

x 6 long stick 

(Expressions for e are similar to those for e . ) 

y x 


( B2 ) 


Where: 

s.f. 

K „ 
pos 


= Scale factor 
= 1 °/1 7 ft 


k 5 = 1.0 

X = True X position 
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K vel = 0.29°/( ft/sec) 
Tjj “■ 1.1 sec 


•• •• 
“• T x + X 

X => -T— -! 

(FsTT) 


True longitudinal 
velocity 

( B3 ) 


-g T 2 s9 

r^FT-TT 


Where: 


(B4) 



X = True longitudinal acceleration 

X = Estimated true acceleration 

X, = Estimate of high frequency acceleration 


Figure B1 shows a graphical definition of the 


error e x and e y . 



Figure B1 : HUD Hover Point Dynamics 




C. DATA OUTPUT OF THE PILOTED SIMULATION 


B2 is an example of to. output data fro. the flx . d . bas0 
simulation. Information used include, the following: 


Page 1 
Page 2 
Page 3 


All of the title information 

- Start and stop times, and T/W information 

- Maximum gear z velocity, and the position 

error when the side task (attitudfcoSand 
system) was being flown. nd 


NOTE: 

Since the simulation was primarily concerned with the vertical 
axis, the highest main gear, designated N (nose) or T (tail), z 
velocity was used, and the outrigger readings were ignored. 

D. PILOTS 

Pilot A - Vernon K. Merrick 
Pilot B - Glenn G. Ferris 


nc 


r.uii 



r.:n r.rct r.nrj rrm r.^ 


rrrm na - c ct a - rzzn 


CVSONG 15:47 JAM 23. '85 RUN 4P7 PILOT: V. MERRICK 


IlfiR AVPA DATA PRINTOUT 


FINAL PRINTOUT- HOVER APPROACH - VFR TEST 

SHIP MOTION ON . IISTATE - I. LOADING • CLEAN , LCIGHT • I75I3.0LB3 

LHLE TUR8ULENCE OH , BACKGROUND TUPB • .CO FT.'SEC 

ICSVS - G ISUESVS • 2 IDCASE ■ 5 IhVPB - 0 


Figure B2: Example Data Output from the Piloted Simulation 
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PHASE in- VERTICAL DESCENT 




r . fyifi - 


TIHE- START AT 

.ee 

END AT 

23.41 

DURATION OF 

28.48 SECS 











ALTITUDE 

AIRSPEED 

ALTITUDE 

H IF SPEED 

DESCRIPTION 

UNITS 

NEON 

STD DEV 

Mill 

MAX 

(FT) MIN 

(F'Sl MIN 

(FT' MAX 

CF tW 

H0R2 VEL 

F.'S 

-.7796E-0 3 

.39 I8E-01 

1804E Oe 

.leioE ee 

.5714E 02 

. 4495E 02 

.4 174E 02 

• 6019b A*. 

VERT VEL 

F.'S 

-.I504E 01 

. 1062E 01 

-.3493E PI 

.usee ei 

.5998E 02 

.3164E e2 

.S063E 02 

.50102 02 

X PILOT ACCEL. 

G 

.1129E 00 

.81 18E-02 

•6923E-01 

. 1434E 00 

•5759E 02 

.5253E 02 

. 57 1 JE 02 

.44996 02 

Y PILOT ACCEL. 

G 

. 1554E-P2 

.2276E-01 

-.5983E-01 

.E765.E-P1 

. 4S20E 02 

. S54SE 02 

•4302E 02 

. 4665c 02 

Z PILOT ACCEL. 

G 

-.99 JOE 00 

.4344E-01 

-. 1 12 IE 01 

-.7716E Ce 

.5024E 02 

•5347E 02 

4009E 02 

.4264-: e*;- 

PC DEFLECTION 

% 

•2271E 02 

. 3267E 01 

.1165E 02 

•3462E 02 

■5976E 02 

■3372E 02 

. 4044E 02 

.2362F 0.3 

RC DEFLECTION 

% 

. 1S5PE 00 

. I2S4E 02 

-.4S79E 02 

.5 19 IE 02 

. 42S8E 02 

.50:'7E 02 

. 449 1 E 02 

•STOOL 02 

YC DEFLECTION 

V, 

. 1863E Pi 

•. 126SE 02 

-.3888E e2 

.3736E 02 

•5549E 02 

. 3333E 02 

.5M1E 02 

.3*3 1 6E 02 

AERO 

R.'S2 

•5195E-0! 

. 20 I0JZ-0 1 

. 1502E-P1 

. 193 ee ee 

.4920E 02 

. 3S36E 02 

.57S9C 02 

. 5"'5;£ t 2 

RC AERO 

R'.'S? 

-.203 IE-04 

. 17S8E-01 

-.S3S3E-P1 

•8653E-P1 

•5774E 02 

.71 ME 02 

.4431E 02 

.5700? 02 

YC AERO 

R.'S? 

-.3597E-P3 

. 1943E-02 

-.6337E-P2 

. 1C91E-P1 

.5125E 02 

. 4266E 02 

.5T5SE 02 

.5253c 73 

PC ENGINE 

R.'S2 

I02BE 00 

.3509E-01 

-.217 IE 00 

.519IE-01 

. 4043E 02 

.4220E 02 

•5354E 02 

.507.1= 07. 

RC ENGINE 

R'S? 

.4S2GE-P2 

.257 IE 00 

-.7964E 00 

.8C39E 

■4274E 02 

•330CE 02 

. 4476E 02 

. 479Sc 02 

YC ENGINE 

R.'S? 

-.7042E-P2 

.4G3y£-& 1 

-. 1839E OP 

•1750E PC 

.5 MCE 02 

.4220E 02 

. 5544C C2 

.3437C 0.3 

TOT ENGINE BLEED 

LB'S 

.3948E 01 

. 1S94E 01 

.7099E 00 

■1214E 02 

.5042E 02 

.2957E 02 

.5 M1F 02 

•3S16E 0.3 

Rpri 

PCT 

. 1002E 03 

. 1 1S0E 01 

.9157E 02 

. 102SE 03 

.3994E 02 

.4476E 02 

. 4033E 02 

. 4506E 02 

T.'U 


•9334E 00 

.3284E-01 

.739GE eP 

. 1C48E 01 

.393-1E 02 

.4476E 02 

•4033E 02 

.5MTE 0.3 

NOZZLE ANGLE 

DEG 

.7837E 02 

.3S78E 00 

. 773SE e2 

.eC02E 02 

.57 me ea 

.4435E 02 

. 4174£ 02 

-6019E 03 

THETA 

DEG 

.SSOOE 0! 

.3220E-01 

•5428E ei 

.eeioE ci 

. 4053E 02 

.412GE 02 

. 4044E 02 

. 23C2F 0.3 

PHI 

DEG 

-.94SSE-ei 

•68SSE 00 

-.1634E 01 

. ieesE oi 

.449 IE 02 

-S700E 02 

.4797E 02 

.42622 0.3 

FLAP DEFL 

DEG 

.500BE e2 

.S840E-0S 

.5000E 02 

.5C00E 02 

. 8250E 02 

•4330E 02 

. 02 SOE 02 

. 4330E 0.3 

T.'U 


.00 

.18 

.34 

.43 

.60 

.70 

.79 

.84 

TIME ABOVE T.'U 

SEC 

•2835E e2 

•283SE 02 

. 283SE 02 

.2035E 02 

. 2835E 02 

.2835E 0? 

•2822E 02 

. 230“E 0.3 

T.'U 


.08 

.90 

.92 

.9-1 

.96 

.97 

.93 

.99 

Tir^E ABOVE T.'U 

SEC 

■279CE P2 

■2784E 02 

•2776E 02 

. 274CS 02 

.2592E 02 

. 2035E 02 

. 1440E 02 

. 9S56E •*! 

T.'U 


1. 00 

1.01 

1.02 

1.03 

1 .04 

1.06 

1.03 

1.10 

TIME ABOVE T.'U 

SEC 

.6400E 01 

.4SS4E 01 

. 3264E 01 

.21 12E 01 

.7630E OP 

.OOOOE 00 

■.OOOOE 00 

.OOOOE 

RELATIVE DESCENT 

RATE 









.TO DECK 

F-'S 

. I532E 01 

•6356E 01 

-.9101E 02 

. 1 193E 02 

. 8250E 02 

.4618E 02 

.62T5E 02 

. 3S49E 0.3 

TO SEA LEVEL 

F.'S 

. 1S04E 01 

. 1062E 01 

-.11S6E 01 

.3493E 01 

. 50b3E 02 

.5010c 02 

.539SE 02 

.3I64E 0.3 

UATER USED * 

.0 LBS. 

FUEL USED - 

113.3 LBS. UE1GHT • 

17399. 1LBS 






Figure B2, continued: Example Data Output from the Piloted Simulation 
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TOUCHDO'JI 

at TIT£ ■ 2e.4l 







VEL - DECK AXES 

position error 

RX ATTITUDE 

SHIP ATTITUDE 

SHIP VElOCir,' 

X- .2P4SE PO F-S 

X- .*?FSE CO FT 

PHI- , 

, 1I44E PI DEG 

PHI- 

.27C6E PI DEG 

y* 

.2I07E 03 F X 

Y- .2-455t-02 

V- -.3?32£ 01 

THT- . 

•EEEOE PI 

THT- 

.6D5ZE PP 

T ■ 

. I I2.*E 02 

Z- .I3ZIE PI 

BIST’ .3E63E 01 

PSI- 

, 3CM3E PZ 

PSI- 

. 4S00E PP 


-.34 LIE 02 

WX OLEO FGRCES 

ATTITUDE-PEL. DECK 

!«/ GEAR X VEL. 

HA/ CEAP Y VEL. 

isc: 

GC»P 3 'TL 

II- .7701E P4 LBS 

PHI- J3ISE PP DEC 

H* 

.3u5'E PP F'S 

(1- 

.HOSE PI F S 

Ih 

. 3V3'S PI F S 

R- . I6S3E P4 

THT- . ^E35E PI 

R- 

. 2202E TP . 

R- 

.3IS4E PP 


.f,;::e pi 

t- . 1364E C3 

PSI- -.3CZ4E P2 

L- 

. 1 1 ODE PP 

L- 

•3IS4E PC 

L - 

.2: :;F 01 

T- .HIDE P5 


T* 

.228CE PP 

T * 

.90rCE PP 

T- 

.4IS2E PI 

the niNitun mss distance las low egual to 

U.4 FT. 
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Figure B2, continued 


Example Data Output Erora the Piloted Simulation 
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E. LISTING OF PILOTED SIMULATION RUNS AND RESULTS 
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PILOTED SIMULATION TEST MATRIX 


T/U -1.10 

eng 0.3 0.7 

v X 

T/V - 1.07 

eng 0.3 0.7 

0.0 

0 5-B 
A 5-B 

□ 5-A, 10-B 

OS-B 

OlO-B 

O 5-B 
A 5-B 

□ 5-A, 10-B 
O 5-B 
O 5-B 

0.0 

A 5-B 

□ 5-A. 10-8 
O 5-A, 5-B 

A 5-B 

□ 5-A, 5-B 
O 5 -A, 5-B 

-0.2 

05-B 
A 5-B 

O 5-A, 10-B 
O 5-B 
OlO-B 

O 5-B 
A 5-B 

□ 5-A, 10-B 
O 6-3 
O 5-B 

-0.2 

A 5-B 

□ 5-A, 10-B 
O 5-A, 5-B 

A 5-B 

□ 5-A, 5-3 
O 5-A, 5-B 

-0.4 

O 5-B 
A 5-B 

O 5-A, 10-B 

05-B 

OlO-B 

O 5-B 
A 5-B 

O 5-A, 10-B 
O 5-B 
O 5-B 

-0.4 

A 5-B 
O 5-A, 5-B 
O 5-A, 5-B 

A 5-B 
O 5-A, 5-B 
O 5-A, 5-B 

\ 

T/U - I. OS 

eng 0.3 0.7 


T/U - l.C 
eng 0.3 

n 

0.7 

0.0 

06-3 
A 5-B 

Q 5-A, 5-B- 
O 5-B 
OS-B 

0-5-3 
A 5-B 
□ 5-A, 5-B 
O 5-B 
OS-B 

0.0 

A 5-B 

□ 15-A, 10-B 
O 6-B 

A 5-B 

□ U-A, 10-B 
O 5-B 

-0.2 

O 5-B 
A 5-B 

□ 5-A, 5-B 
OS-B 

O 5-B 

O 5-3 
A 5-3 

□ 5-A, 5-B 
O 5-B 
O 4-B 

-0.2 

A 5-B 

□ 15-A, 15-B 
O 5-A, 10-B 

A 5-B 

□ 15-A, 15-B 
O 5-A, 10-B 

-0.4 

O 5-B ’ 

A 5-B 

D 6- A, 5-B 
O 5-B 
OS-B 

O 5-B 
A 5-B 
□ 5-A, 5-B 
O 5-B 
O 5-B 

-0.4 

A 5-B 

□ 15-A, 10-B 
O 5-B 

A 5-B 

□ 10-A, 10-B 
O 5-B 


T/U - 1.01 

eng 0.3 0.7 




0.0 

O 5-B 
A 5-B 

□ 5-A, 5-B 
O 5-A, 5-B 
O 5-B 

O 5-B 
A 5-B 
□ 5-A, 5-B 
O 5-A, 5-B 
O 5-B 


O n-l S.S.O 
A n-l S.S.4 

BUD3 V.C. 
HUD3 V.C. 

-0.2 

0 5-3 
A 5-B 

OS-A, 5-B 
O 5-A, 5-B 
OS-B 

O 5-B 
A 5-B 
□ 6-A, 6-B 
O 5-A, 5tB 
O 5-B 


□ n-l S.S.6 11UD3 V.C. 

O n-l S.S.6 HUD1 V.C. 
O n-t S.S.6 HUD3 A.C. 

n- Number of runs 
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-0.4 

0 5-3 
A 5-B 
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05-5 
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A 5-B 

□ 5-A, 5-B 
O 5-A, 5-B 
O 5-B 





Figure B3: Piloted Simulation Test Matrix 
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PILOTED SIMULATION RUN LISTING 


SEA STATE 0 
HUD3 

VELOCITY COMMAND SYSTEM 


T 

eng 

Z 

w 


- (sec) 

- (sec -1 ) 


T.D. v j- (ft/sec) 







Flight Time - 

(sec) 

T/W 

max 

t 

eng 

z 

w 

Pilot 

Run No . 3 

T.D. , 
vel 

FLICHT 

TIME 

P.R. 

1.01 

0.3 

0.0 

B 

916-920 

1.87 

26.0 

3 

II 

II 

0.2 

II 

911-915 

1.70 

22.9 

2 

II 

II 

0.4 

It 

906-910 

2.17 

25.0 

1 1/2 

II 

0.7 

0.0 

11 

891-895 

3.15 

20.2 

3 1/2 

II 

II 

0.2 

11 

896-905 

1.42 

27.8 

2 1/2 

II 

• 1 

0.4 

II 

901-905 

1.58 

31.5 

1 1/2 

1.05 

0.3 

0.0 

II 

921-925 

2.04 

23.3 

3 

II 

II 

0.2 

II 

926-930 

1.61 

21.6 

2 

II 

II 

0.4 

19 

931-935 

1.36 

28.0 

1 

II 

0.7 

0.0 

II 

946-950 

2.04 

17.5 

3 1/2 

It 

II 

0.2 

It 

941-945 

1.14 

28.4 

2 1/2 

11 

II 

0.4 

It 

936-940 

1.57 

28.0 

1 1/2 

1.10 

0.3 

0.0 

II 

976-980 

1.95 

15.1 

3 

II 

II 

0.2 

It 

971-975 

1.26 

23.8 

1 1/2 

II 

11 

0.4 

II 

966-970 

0.68 

42.1 

1 

II 

0.7 

0.0 

1* 

951-955 

1.49 

24.2 

3 1/2 

II 

II 

0.2 

II 

956-960 

1.29 

33.2 

3 

II 

II 

0.4 

• 1 

961-965 

1.21 

28.4 

2 

SEA STATE 4 







11UD3 








VELOCITY 

COMMAND 

SYSTEM 






1.01 

0.3 

0.0 

It 

861-865 

2.9 

26.1 

4 1/2 

II 

ft 

0.2 

11 

866-870 

4.0 

23.2 

3 1/2 

It 

II 

0.4 

11 

871-875 

2.8 

30.0 

3 

II 

0.7 

0.0 

It 

886-890 

3.5 

27.2 

6 

II 

II 

0.2 

II 

881-S85 

3.7 

24.1 

5 

II 

II 

0.4 

<1 

876-880 

3.4 

33.1 

4 

1.03 

0.3 

0.0 

II 

856-860 

3.0 

27.7 

3 1/2 

11 

It 

0.2 

II 

851-855 

2.3 

24.0 

3 1/2 

II 

II 

0.4 

H 

846-850 

1.9 

20.1 

2 1/2 

II 

0.7 

0.0 

• 1 

831-835 

3.9 

17.2 

4 1/2 

II 

II 

0.2 

II 

836-840 

2.8 

19.7 

3 1/2 

II 

II 

0.4 

II 

841-845 

3.5 

23.0 

3 1/2 

1.05 

0.3 

0.0 

II 

801-805 

4.6 

14.7 

3 1/2 

• 1 

II 

0.2 

11 

806-810 

3.4 

17.0 

3 1/2 

II 

II 

0.4 

II 

811-815 

5.3 

23.3 

4 

II 

0.7 

0.0 

II 

826-830 

5.4 

17.3 

5 

II 

II 

0.2 

II 

821-S25 

4.0 

27.7 

4 1/2 

II 

II 

0.4 

11 

816-820 

3.7 

21.9 

4 
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PILOTED SIMULATION- RUN LISTING Cont. 


SEA STATE 4 Cont. 

HUD3 

VELOCITY COMMAND SYSTEM 


T/W 

max 

T 

eng 

Z Pilot 

w 

Run No. 3 

T.D. , 
vel 

FLIGHT 

TIME 

P.R. 

1.07 

0.3 

0.0 

B 

796-800 

2.3 

15.9 

3 

tl 

II 

0.2 

II 

786-790 

2.8 

17.3 

3 

It 

•I 

0.4 

If 

791-795 

3.0 

20.9 

2 1/2 

ft 

0.7 

0.0 

It 

771-775 

4.2 

21.3 

5 

•1 

II 

0.2 

II 

776-780 

3.7 

20.5 

4 

II 

II 

0.4 

II . 

781-785 

4.4 

21.2 

4 

1.10 

0.3 

0.0 

tl 

741-745 

3.5 

14.7 

4 

II 

II 

0.2 

II 

746-759 

3.5 

18.3 

3 

II 

II 

0.4 

If 

751-755 

3.5 

21.3 

1 1/2 

•1 

0.7 

0.0 

II 

766-770 

3.6 

18.6 

4 

It 

»• 

0.2 

ft 

761-765 

4.0 

22.3 

3 

II 

II 

0.4 

It 

756-760 

3.7 

22.6 

2 

SEA STATE 6 







HUD1 








VELOCITY 

COMMAND 

SYSTEM 






1.01 

0.3 

0.0 

A 

506-510 

8.0 

26.2 

7 

II 

II 

II 

B 

616-620 

11.6 

26.0 

10 

ft 

It 

0.2 

A 

511-515 

6.6 

29.3 

6 1/2 

ir 

II 

II 

B 

611-615 

9.4 

22.7 

7 1/2 

»» 

II 

0.4 

A 

431-435 

6.2 

39.0 

5 

•it 

If 

It 

B 

606-610 

6.9 

31.0 

7 

tl 

0.7 

0.0 

A 

526-530 

6.5 

24.9 

7 

if 

II 

II 

B 

591-595 

7.2 

18.2 

8 1/2 

if 

II 

0.2 

A 

521-525 

8.6 

25.1 

6 3/4 

ft 

II 

II 

B 

596-600 

9.1 

15.4 

8 1/2 

ii 

If 

0.4 

A 

516-520 

5.3 

29.0 

5 1/2 

ft 

II 

If 

B 

601-605 

7.9 

24.5 

8 1/2 

1.03 

0.3 

0.0 

B 

561-565B 

9.3 

15.4 

8 

II 

II 

0.2 

A 

411-415 

6.2 

23.0 

5 

II 

II 

' tl 

B 

416-420 

8.1 

24.1 

5 

ft 

It 

II 


566-570 

8.7 

18.9 

7 1/2 

If 

II 

0.4 

tl 

571-575 

6.7 

29.9 

7 

II 

0.7 

0.0 

If 

586-590 

9.6 

25.9 

8 1/2 

II 

II 

0.2 

II 

421-415 

8.6 

20.5 

5 1/2 

If 

II 

If 

A 

426-430 

9.1 

20.9 

6 

II 

It 

II 

B 

581-585 

9.5 

17.7 

8 

II 

• 1 

0.4 

If 

576-580 

7.2 

16.7 

7 
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PILOTED SIMULATION RUN LISTING Cont 


SEA STATE 6 Cont. 

HUD3 

VELOCITY. COMMAND SYSTEM 


T/W 

max 

T 

eng 

Z 

w 

Pilot 

Run No. a 

T.D. . 
vel 

Flight 

Time 

P.R. 

1.05 

0.3 

0.0 

B 

621-625 

6.2 

30.0 

6 

II 

II 

0.2 

II 

626-630 

7.8 

32.3 

6 

II 

II 

0.4 

tl 

631-635 

7.0 

34.3 

5 3/2 

II 

0.7 

0.0 

It 

636-640 

6.4 

22.7 

6 

II 

II 

0.2, 

If 

641-645 

9.2 

39.8 

7 


If 

0.4 

II 

646-650 

7.2 

24.1 

7 

1.07 

0.3 

0.0 

A 

501-505 

8.2 

27.9 

6 

II 

II 

II 

B 

661-665 

10.0 

16.8 

7 

II 

II 

0.2 

A 

496-500 

8.2 

21.2 

5 1/4 

It 

It 

II 

B 

656-660 

4.5 

28.7 

3 1/2 

II 

II 

0.4 

A 

491-495 

5.5 

32.5 

4 3/4 

II 

It 

II 

B 

651-655 

5.3 

22.5 

3 1/2 

It 

0.7 

0.0 

A 

476-480 

6.2 

20.5 

5 1/2 

II 

II 

II 

B 

666-670 

8.1 

13.1 

7 

II 

• 1 

0.2 

A 

481-485 

7.2 

17.0 

5 1/2 


If 

II 

B 

671-675 

8.1 

20.4 

6 

l» 

II 

0.4 

A 

486-490 

5.3 

23.0 

5 

If 

II 

II 

B 

676-680 

7.9 

34.5 

6 

1.10 

0.3 

0.0 

B 

696-700 

7.6 

19.1 

7 

It 

If 

0.2 

If 

701-705 

7.0 

35.4 

5 1/2 

II 

It 

0.4 

II 

706-710 

5.7 

38.1 

4 

If 

0.7 

0.0 

II 

691-695 

10.2 

22.3 

7 1/2 

II 

If 

0.2 

II 

686-690 

7.1 

43.4 

6 1/2 

ft 

II 

0.4 

II 

681-685 

7.3 

43.8 

5 

SEA STATE 

6 







HUD3 








VELOCITY 

COMMAND 

SYSTEM 






1.01 

0.3 

0.0 

A 

156-160 

6.1 

33.3 

5 

If 

If 

II 

B 

186-190 

6.5 

24.2 

6 1/2 

II 

If 

0.2 

A 

151-155 

5.7 

34.1 

4 1/2 

It 

II 

II 

B 

181-185 

8.5 

34.9 

6 1/2 

II 

II 

0.4 

A 

146-150 

5.6 

31.2 

4 1/4 

II 

II 

II 

B 

176-180 

9.4 

18.3 

5 1/2 

II 

0.7 

0.0 

A 

171-175 

6.7 

28.5 

5 

II 

II 

II 

B 

201-205 

6.4 

35.4 

6 1/2 

If 

II 

0.2 

A 

166-170 

6.5 

26.1 

5 


II 

II 

B 

196-200 

7.9 

29.6 

7 

It 

If 

0.4 

A 

161-165 

7.7 

34.2 

5 1/2 


II 

II 

3 

191-195 

6.7 

38.6 

5 1/2 


B1 '3 



PILOTED SIMULATION RUN LISTING Cant. 


SEA STATE 6 Cont. 

MJD3 

VELOCITY COMMAND SYSTEM 


T/W 

IuaX 

T 

eng 

z 

w 

Pilot 

Run No.s 

T.D. . 
vel 

Flight 

Time 

P.R. 

1.03 

0.3 

0.0 

A 

55-1 00 

7.0 

37.7 

4 1/2 

It 

II 

II 

II 

131-135 

7.1 

28.7 

5 

It 

It 

II 

II 

271-275 

6.7 

22.0 

5 

M 

II 

II 

B 

291-295 

5.5 

21.2 

5 1/2 

II 

It 

II 

II 

556-560 

6.0 

23.8 

5 

II 

It 

0.2 

A 

91-95 

6.6 

32.3 

4 1/4 

II 

tl 

II 

II 

136-140 

4.2 

32.2 

4 1/2 

II 

II 

II 

B 

286-290 

6.5 

26.6 

5 

II 

• 1 

tt 

A 

401-405 

6.5 

33.1 

5 1/4 

»» 

II 

II 

B 

406-410 

5.3 

44.0 

5 

II 

It 

tl 

It 

551-555 

5.1 

19.1 

6 

II 

It 

0.4 

A 

86-90 

7.4 

38.3 

4 1/2 

It 

It 

tt 

tl 

141-145 

6.5 

24.5 

4 1/4 

II 

tl 

II 

tt 

266-270 

7.1 

42.4 

5 1/4 

II 

It 

It 

B 

281-285 

5.4 

26.8 

4 1/2 

It 

It 

It 

II 

546-550 

7.5 

29.8 

6 1/2 

II 

0.7 

0.0 

A 

111-115 

6.5 

38.9 

4 3/4 

It 

II 

It 

It 

116-120 

6.3 

37.1 

4 3/4 

It 

II 

II 

B 

306-310 

7.2 

17.1 

6 1/2 

It 

It 

tt 

II 

531-535 

7.6 

23.3 

7 

II 

II 

0.2 

A 

106-110 

8.0 

36.7 

4 3/4 

It 

II 

It 

It 

121-125 

6.3 

25.6 

4 1/2 

II 

It 

II 

tl 

276-280 

7.7 

22.6 

5 1/2 

II 

II 

II 

B 

301-305 

5.2 

19.9 

5 1/2 

It 

It 

• 1 

II 

536-540 

4.7 

24.5 

5 1/2 

II 

tl 

It 

A 

101-105 

8.1 

35.8 

4 1/2 

It 

tl 

II 

II 

126-130 

7.1 

33.7 

4 1/2 

II 

II 

II 

B 

296-300 

5.8 

24.6 

5 1/2 

tl 

It 

II 

II 

541-545 

4.9 

25.7 

5 

1.05 

0.3 

0.0 

A 

216-220 

8.9 

24.2 

5 1/2 

II 

II 

It 

B 

251-255 

7.0 

30.7 

5 

II 

tl 

0.2 

A 

211-215 

4.4 

26.6 

4 1/2 

II 

It 

II 

B 

256-260 

5.5 

33.6 

5 

It 

tl 

0.4 

A 

206-210 

6.6 

27.7 

5 

It 

• 1 

II 

B 

261-265 

7.1 

42.4 

4 1/2 

II 

0.7 

0.0 

A 

221-225 

7.1 

31.5 

5 1/2 

II 

It 

It 

B 

246-250 

5.0 

44.5 

6 1/2 

II 

II 

0.2 

A 

226-230 

7.0 

38.3 

5 

II 

II 

II 

B 

241-245 

8.0 

37.2 

6 1/2 

II 

II 

0.4 

A 

231-235 

6.5 

30.9 

5 

II 

II 

II 

B 

236-240 

9.8 

39.9 

6 1/2 




PILOTED SIMULATION RUN LISTING Cont. 

SEA STATE 6 Cont. 

HUD3 

VELOCITY COMMAND SYSTEM 


l- 4 

l * 

T/W 

max 

T 

eng 

z 

w 

Pilot 

Run No.s 

T.D. , 
vel 

Flight 

Time 

P.R. 

V 

1.05 

0.3 

0.0 

A 

216-220 

8.9 

24.2 

5 1/2 


M 

tl 

II 

B 

251-255 

7.0 

30.7 

5 


It 

II 

0.2 

A 

211-215 

4.4 

26.6 

4 1/2 

: 

II 

II 

II 

B 

256-260 

5.5 

33.6 

5 

. 

It 

II 

0.4 

A 

206-210 

6.6 

27.7 

5 


II 

It 

tl 

B 

261-265 

7.1 

42.4 

4 1/2 

| 

• 1 

0.7 

0.0 

A 

221-225 

7.0 

31.5 

5 1/2 

\ 


II 

II 

B 

246-250 

5.0 

44.5 

6 1/2 


It 

It 

0.2 

A 

226-230 

7.0 

38.3 

5 

( 

II 

It 

If 

B 

241-245 

8.0 

37.2 

6 1/2 

| 

II 

tl 

0.4 

A 

231-235 

6.5 

30.9 

5 

1 

II 

II 

II 

B 

236-240 

9.2 

39.9 

6 1/2 

P 

1.07 

0.3 

0.0 

B 

386-390 

5.0 

24.7 

4 1/2 

L 

II 

It 

II 

<i 

436-440 

5.8 

29.5 

4 1/2 



II 

II 

A 

446-450 

5.9 

25.2 

5 7/4 

) ' 

It 

If 

0.2 

B 

391-395 

7.8 

22.2 

4 

| 

It 

II 

ti 

•t 

441-445 

6.1 

26.7 

3 1/2 



ft 

tt 

A 

451-455 

5.9 

31.3 

5 

r 

II 

tl 

0.4 

B 

396-400 

5.5 

26.2 

3 1/2 


II 

tt 

ti 

A 

456-460 

6.0 

30.. 3 

4 1/2 

1 

It 

0.7 

0.0 

B 

381-385 

7.4 

23.5 

5 1/2 



If 

tt 

A 

471-475 

5.4 

39.9 

6 

r 

II 

If 

0.2 

B 

376-380 

7.0 

27.7 

4 1/2 

1 

II 

II 

tt 

A 

466-470 

7.4 

35.6 

5 1/2 


tl 

tl 

0.4 

B 

371-375 

7.1 

24.7 

4 1/2 

\ * 

II 

If 

tt 

A 

461-465 

5.7 

35.4 

4 3/4 

I 

1.10 

0.3 

0.0 

B 

321-325 

6.7 

33.2 

3 1/2 


If 

If 

It 

A 

356-360 

5.3 

32.5 

3 3/4 

J 

tl 

If 

tt 

B 

721-725 

5.2 

26.9 

4 

L 


II 

0.2 

<i 

316-320 

5.3 

25.8 

3 


II 

It 

tt 

A 

361-365 

5.9 

29.3 

4 

j 

II 

It 

tt 

B 

716-720 

4.5 

52.1 

3 1/2 

1 


• 1 

0.4 

ti 

311-315 

8.2 

40.2 

3 



tl 

it 

A 

366-370 

6.2 

26.2 

4 

i . 

II 

If 

it 

B 

711-715 

4.7 

53.9 

2 1/2 

i 

i 

ft 

0.7 

0.0 

it 

336-340 

6.9 

32.6 

5 

/ 

ft 

II 

tl 

A 

351-355 

6.9 

34.4 

4 1/2 


If 

II 

It 

B 

726-730 

5.6 

39.0 

4 1/2 

1 

II 

ft 

0.2 

ii 

331-335 

4.4 

29.2 

3 

p 

It 

tt 

it 

A 

346-350 

6.1 

35.0 

4 


ft 

If 

•i 

B 

731-735 

5.4 

22.3. 

4 

i 

II 

ft 

0.4 

it 

326-330 

6.4 

40.1 

3 

» 

i 


It 

ti 

A 

341-345 

5.2 

48.9 

4 

1 .:. 

ft 

II 

tt 

B 

736-740 

5.0 

33.7 

5 1/2 
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PILOTED SIMULATION RUN LISTING Cont. 

SEA STATE 6 
HUD3 

ATTITUDE COMMAND SYSTEM 


T/W max T eng Z w Pllot *»“ N °* s T * D * ve l F1± B ht P - R - 

Time 


1.01 

0.3 

0.0 

B 

1051-1055 

8.7 

36.9 

8 

II 

II 

0.2 

ti 

1046-1050 

5.7 

30.0 

7 

ft 

II 

0.4 

it 

1041-1045 

6.6 

29.6 

6 

It 

0.7 

0.0 

it 

1056-1060 

6.3 

25.3 

8 

II 

If 

0.2 

ii 

1061-1065 

5.9 

35.3 

6 1/2 

II 

II 

0.4 

it 

1066-1070 

6.9 

39 .'6 

6 

1.05 

0.3 

0.0 

ti 

1026-1030 

9.5 

24.8 

8 

II 

II 

0,2 

it 

1031-1035 

7.5 

29.6 

6 1/2 

II 

It 

0.4 

ii 

1036-1040 

4.8 

35.9 

5 1/2 

II 

0.7 

0.0 

ii 

1021-1025 

8.7 

27.9 

8 

II 

II 

0.2 

it 

1016-1020 

9.9 

22.2 

7 

II 

II 

0.4 

ii 

1011-1015 

6.3 

27.3 

6 

1.10 

0.3 

0.0 

ti 

1081-1085 

5.0 

33.4 

5 

II 

If 

0.2 

it 

1076-1080 

6.8 

34.8 

5 

It 

II 

0.0 


1071-1075 

5.4 

49.2 

4 1/2 

II 

0.7 

0.0 

ti 

996-1000 

7.3 

23.0 

7 1/2 

II 

II 

0.2 

it 

1001-1005 

6.7 

31.6 

6 1/2 

If 

tl 

0.4 

it 

1006-1010 

6.7 

29.0 

6 
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ABSTRACT 

The problem of determining the vertical axi3 control 
requirements for landing a VTOL aircraft on a moving ship deck in 
various sea states i3 examined. Both a fixed-base piloted 
simulation and a non-rpiloted simulation were used to determine the 
landing performance as influenced by thrust-to-weight ratio, 
vertical damping, and engine lags. 

The piloted simulation was run using a fixed-base simulator at 
N.A.S.A. Ames Research Center, simplified versions of an existing 
AV-8A Harrier model and an existing head-up display format were 
used. The ship model used wan that of a DD963 class destroyer. 

Simplified linear models of the pilot, aircraft, ship notion, 
and ship air-wake turbulence were developed for the non-piloted 
simulation. A unique aspect of the non-piloted simulation was the 
development of a model of the piloting strategy used for shipboard 
landing. This model was refined during the piloted simulation until 
it provided a reasonably good representation of observed pilot 
behavior. Further refinement could lead to a model suitable for 
prediction of landing performance of VTOL aircraft on ships and as 
the basis of control logic for automatic landing. 

A surprising result of this simulation was that, with a good 
station keeping control system and with statistical ship motion 
displayed on the head-up display, pilot3 could consistently perform 
safe landings in sea state 6, with handling qualities that were 
adequate at thru3t-to-weight ratios greater than 1.03 and even . 
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marginally adequate down to thrust-to-weight ratios of 1.01. These 
results should hold quite generally provided that a thrust- to-weight 
ratio' of 1 + A is interpreted a3 meaning that the pilot always has 
the capability of accelerating the aircraft at Ag upward even in the 
presence of ground effect and hot gas reingestion. 
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